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Abstract 



The aim of this thesis is twofold. First, to study methods to calculate one-loop 
corrections in the context of perturbative theories. Second, to apply those methods 
to calculate the leading electroweak (EW) corrections to the important process of 
Higgs production associated with two bottom quarks at the CERN Large Hadron 
Collider (LHC). Our study is restricted to the Standard Model (SM). 

The first aim is of theoretical importance. Though the general method to cal- 
culate one-loop corrections in the SM is, in principle, well understood by means of 
renormalisation, it presents a number of technical difficulties. They are all related 
to loop integrals. The analytical method making use of various techniques to reduce 
all the tensorial integrals in terms of a basis of scalar integrals is most widely used 
nowadays. A problem with this method is that for processes with more than 4 exter- 
nal particles the amplitude expressions arc extremely cumbersome and very difficult 
to handle even with powerful computers. In this thesis, we have studied this problem 
and realised that the whole calculation can be easily optimised if one uses the helicity 
amplitude method. Another general problem is related to the analytic properties of 
the scalar loop integrals. An important part of this thesis is devoted to studying this 
by using Landau equations. We found significant effects due to Landau singularities 
in the process of Higgs production associated with two bottom quarks at the LHC. 

The second aim is of practical (experimental) importance. Higgs production asso- 
ciated with bottom quarks at the LHC is a very important process to understand the 
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bottom-Higgs Yukawa coupling. If this coupling is strongly enhanced as predicted by 
the Minimal Supersymmetric Standard Model (MSSM) then this process can have 
a very large cross section. In this thesis, based on the theoretical study mentioned 
above, we have calculated the leading EW corrections to this process. The result 
is the following. If the Higgs mass is about 120GeV then the next-to-leading order 
(NLO) correction is small, about —4%. If the Higgs mass is about 160GeV then the 
EW correction is strongly enhanced by the Landau singularities, leading to a signif- 
icant correction of about 50%. This important phenomenon is carefully studied in 
this thesis. 
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Introduction 



In the realm of high energy physics, the Standard Model (SM) of particle physics 
[HElEllllElEllTj is the highest achievement to date. Almost all its predictions have 
been verified by various experiments [HI [9]. The only prediction of the SM which has 
not been confirmed by any experiment is the existence of a scalar fundamental particle 
called the Higgs boson. The fact that we have never observed any fundamental scalar 
particle in nature so far makes this the truly greatest challenge faced by physicists 
today. For this greatest challenge we have the world largest particle accelerator to 
date, the CERN Large Hadron Collider (LHC) [TU]. The LHC collides two proton 
beams with a center-of-mass energy up to 14TeV and is expected to start this year. 
It is our belief that the Higgs boson will be found within a few years. 

The prominent feature of the Higgs boson is that it couples mainly to heavy 
particles with large couplings. This makes the theoretical calculations of the Higgs 
production rates as perturbative expansions in those large couplings complicated. The 
convergence rate of the perturbative expansion is slow and one cannot rely merely 
on the leading order (LO) result. Loop calculations are therefore mandatory. The 
most famous example is the Higgs production mechanism via gluon fusion, the Higgs 
discovery channel. The LO contribution in this example is already at one-loop level. 
The two-loop contribution, mainly due to the gluon radiation in the initial state and 
the QCD virtual corrections, increases the total cross section by about 60% for a Higgs 
mass about lOOGeV at the LHC [TT]. Indeed, loop calculations are required in order 
to understand the structure of perturbative field theory and the uncertainties of the 
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theoretical predictions. The only way to reduce the error of a theoretical prediction 
so that it can be comparable to the small error (say 10%) of precision measurements 
nowadays is to pick up higher order terms, i.e. loop corrections. 

There are two methods to calculate loop integrals: analytical and numerical meth- 
ods. The traditional analytical method decomposes each Feynman diagram's numer- 
ator into a sum of scalar and tensorial Passarino-Veltman functions. The advantage is 
that the whole calculation of cross sections involving the numerical integration over 
phase space is faster. The disadvantage is that the numerator decomposition usu- 
ally results in huge algebraic expressions with various spurious singularities, among 
them the inverse of the Gram determinant (defined as det(G) = det(2pj.pj) with pi 
are external momenta) which can vanish in some region of phase space. Recently, 
Denner and Dittmaier have developed a numerically stable method for reducing one- 
loop tensor integrals [121 US], which has been used in various electroweak processes 
including the e~^e~ — > 4 fermions process [131 For the numerical method, the 
loop integration should be performed along with the integrations over the momenta 
of final state particles. In this method one should not decompose the various nu- 
merators but rather combine various terms in one common denominator. Thus the 
algebraic expression of the integrand is much simpler this way and no spurious sin- 
gularities appear. The disadvantage is that the number of integration variables is 
large resulting in large integration errors. In both methods, the ultra-violet (UV)-, 
infrared (IR)- and collinear- divergences have to be subtracted before performing the 
numerical integration. 

Recently, there has appeared on-shell methods to calculate one-loop multi-leg QCD 
processes (see [TB] for a review). These methods are analytical but very different from 
traditional methods based on Passarino-Veltman reduction technique. On-shell meth- 
ods have already led to a host of new results at one loop, including the computation of 
non-trivial amplitudes in QCD with an arbitrary number of external legs [T71 [THi [T^ . 
These methods work as follows. A generic one-loop amplitude can be expressed in 
terms of a set of scalar master integrals multiplied by various rational coefficients. 
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along with the additional purely rational terms. The relevant master integrals con- 
sist of box, triangle, bubble and (for massive particles) tadpole integrals. All these 
basic integrals are known analytically. The purely rational terms have their origin 
in the difference between D = A — 2e and four dimensions when using dimensional 
regularization. One way to calculate the rational terms is to use on-shell recursion 
[201 [21] to construct the rational remainder from the loop amplitudes' factorization 
poles [221 dZl dS]- The various rational coefficients are determined by using gener- 
alized unitarity cuts [23l [2^ . The evaluation is carried out in the context of the 
spinor formalism. Like the traditional analytical method, spurious singularities occur 
in intermediate steps. However, it is claimed in [16] that they can be under control. 
More detailed studies on this important issue are necessary to confirm this statement 
though. On-shell methods can also deal with massive internal/external particles [23] 
and hence can be used for electroweak processes. It is not clear for us whether these 
on-shell methods can be extended to include the case of internal unstable particles. 

Although the on-shell methods differ from the traditional analytical methods in 
many respects, they have a common feature that one-loop amplitudes are expressed 
in terms of a set of basic scalar loop integrals. One may wonder if there is a method 
to express a one-loop amplitude in terms of tree-level amplitudes? The answer was 
known 45 years ago by Feynman [26l [27] . Feynman has proved that any diagram 
with closed loops can be expressed in terms of sums (actually phase-space integrals) 
of tree diagrams. This is called the Feynman Tree Theorem (FTT) whose very simple 
proof can be found in ^27j. This theorem can be used in several ways. The simplest 
application is to calculate scalar loop integrals needed by other analytical methods 
described above. The best application is to calculate loop corrections for physical 
processes. Feynman has shown that this important application can be realized for 
many processes. Let us explain this a little bit more. After making use of the FTT, 
one has a lot of tree diagrams obtained by cutting a A^-point one-loop diagram with 
multiple cuts (single-cut, double-cut, . . ., A^-cut). One can re-organize this result as a 
sum of sets of tree diagrams, each set representing the complete set of tree diagrams 
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expected for some given physical process. In this way, one obtains relations among 
the diagrams for various processes. Surprisingly, no one has applied this FTT to 
calculate QCD/EW one-loop corrections to important processes at colliders, to the 
best of our knowledge. However, there is ongoing effort in this direction by Catani, 
Gleisberg, Krauss, Rodrigo and Winter. They have very recently proposed a method 
to numerically compute multi-leg one-loop cross sections in perturbative field theories 
[2H]- The method relies on the so-called duality relation between one-loop integrals 
and phase space integrals. This duality relation is very similar to the FTT. The 
main difference is that the duality relation involves only single cuts of the one-loop 
diagrams. Interestingly, the duality relation can be applied to one-loop diagrams with 
internal complex masses [28] . 

In general, Higgs production processes involve unstable internal particles. If these 
unstable particles can be on-shell then the width effect can be relevant and therefore 
must be taken into account. In particular, scalar box integrals with unstable inter- 
nal particles can develop a Landau singularity (to be discussed below) which is not 
integrable at one-loop amplitude square level. In this case, the internal widths are 
regulators as they move the singularity outside the physical region. Thus, a good 
method to calculate one-loop corrections must be able to handle internal complex 
masses. 

Independent of calculation methods, the analytic structure of S-matrix is intrinsic 
and is related to fundamental properties like unitarity and causality [29]. Analytic 
properties of S-matrix can be studied by using Landau equations [201 [21] applied 
to an individual Feynman diagram. Landau equations are necessary and sufficient 
conditions for the appearance of a pinch singularity of Feynman loop integrals [31j . 
Solutions of Landau equations are singularities of the loop integral as a function of 
internal masses and external momenta, called Landau singularities. These singular- 
ities occur when internal particles are on-shell. They can be finite like the famous 
normal threshold in the case of one-loop two-point function. The normal thresholds 
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are branch points [22] • Landau singularities can be divergent like in the case of three- 
point and four-point functions. The former is integrable but the latter is not at the 
level of one-loop amplitude squared. This four-point Landau divergence can be due 
to the presence of internal unstable particles and hence must be regularized by taking 
into account their widths. A detailed account on this topic is given in chapters H] and 

El 

The main calculation of this thesis is to compute the leading electroweak one-loop 
correction to Higgs production associated with two bottom quarks at the LHC in 
the SM. Our calculation involves 8 tree-level diagrams and 115 one-loop diagrams 
with 8 pentagons. The loop integrals include 2-point, 3-point, 4-point and 5-point 
functions which contain internal unstable particles, namely the top-quark and the 
W gauge boson. Interestingly, Landau singularities occur in all those functions. We 
follow the traditional analytical method of Veltman and Passarino [32j to calculate 
the one-loop corrections. For the 5-point function part, we have adapted the new re- 
duction method of Denner and Dittmaier [12], which replaces the inverse of vanishing 
Gram determinant with the inverse of the Landau determinant and hence replaces 
the spurious Gram singularities with the true Landau singularities of loop integrals. 
In our opinion, this is one of the best ways to deal with those spurious Gram singular- 
ities. However, as will be explained in chapter HI the condition of vanishing Landau 
determinant is necessary but not sufficient for a Landau singularity to actually occur 
in the physical region. Thus, spurious singularities can still be encountered but very 
rarely. This new reduction method for 5-point functions has been implemented in the 
library LoopTools [33l [3l] based on the library FF [35] . Our calculation has proved 
the efficiency of this method. The reason for us to choose this traditional method is 
that our calculation involves massive internal particles. Furthermore, in order to deal 
with Landau singularities, our calculation must include also complex masses. 
Although the calculation method is well understood, the difficulty is that we have to 
handle very huge algebraic expressions since we have to expand the numerator of each 
Feynman diagram. Thus, we cannot use the traditional amplitude squared method 
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as it will result in extremely enormous algebraic expressions of the total amplitude 
squared. Fortunately, there is a very efficient way to organize the calculation based 
on the helicity amplitude method (HAM) (SB]. Using this HAM, one just needs to 
calculate all the independent helicity amplitudes which are complex numbers. This 
way of calculating makes it very easy to divide the whole complicated computation 
into independent blocks therefore factorizes out terms that occur several times in the 
calculation. 

Our calculation consists of two parts. In the ffist part, we calculate the NLO 
corrections, i.e. the interference terms between tree-level and one-loop amplitudes. 
Although Landau singularities do appear in many one-loop diagrams, they are inte- 
grable hence do not cause any problem of numerical instability. The bottom-quark 
mass is kept in this calculation. In the second part, we calculate the one- loop cor- 
rection in the limit of massless bottom-quark therefore the bottom-Higgs Yuakawa 
coupling vanishes. The process is loop induced and we have to calculate one-loop 
amplitude squared. In this calculation, the Landau singularity of a scalar four-point 
function is not integrable and causes a severe problem of numerical instability if 
Mh > 2Mw- This problem is solved by introducing a width for the top-quark and W 
gauge boson in the loop diagrams. It turns out that the width effect is large if Mh is 
around 2Mw- 

Although the main calculation of this thesis is for a very specific process, we 
have gained several insights that can be equally used for other practical calculations. 
First of all, the method to optimise complicated loop calculations using the HAM is 
general. Second, the method to check the final/intermediate results by using QCD 
gauge invariance in the framework of the HAM can be used for any process with at 
least one gluon in the external states. Third, some general results related to Landau 
singularities are new and can be used for practical purposes. They are equations 
f l4.27p and (14.491) . Finally, we have applied the loop calculation method of 't Hooft 
and Veltman to write down explicitly two formulae to calculate scalar box integrals 
with complex internal masses. They are equations flE.lSp and flE.40p . The restriction 
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is that at least two external momenta are lightlike. We have implemented those two 
formulae into the library LoopTools. 

The outline of this thesis is as follows. First, a short review of the SM including 
QCD is presented in chapter [H We pay special attention to the one-loop renormali- 
sation of the EW part and the Higgs sector. We also give a short introduction to the 
Minimal Supersymmetric Standard Model and discuss its Higgs sector in the same 
chapter. In chapter [2] we discuss the dominant mechanisms for SM Higgs production 
at the LHC and Higgs signatures at the colliders. In chapter [3] we present the main 
calculation of this thesis, one-loop Yukawa corrections to the SM process pp —>■ bbH 
at the LHC, for the case Mh < ISOGeV. There are two reasons to start with small 
values of the Higgs mass: it is preferred by the latest EW data and there is no problem 
of numerical instability related to Landau singularities. The framework of a one-loop 
calculation based on the helicity amplitude method is also given in this chapter. In 
chapter H] we explain in detail the Landau singularities of a general one-loop Feyn- 
man diagram. We emphasize the conditions to have a Landau singularity and its 
nature. In chapter El we complete the study of chapter [S] for larger values of Mh, 
up to 250GeV. We show that the one-loop process gg —>■ bbH is an ideal example for 
understanding Landau singularities. It contains several types of Landau singularities 
related to two-point, three-point and four-point functions. The conclusions are given 
in chapter O 

This thesis includes several appendices. In appendix El we explain the helicity am- 
plitude method and how to check the correctness of the result by using QCD gauge 
invariance. In appendix[B]we show how to optimise the calculation of various one-loop 
helicity amplitudes and how that can be easily achieved by using FORM. Appendix 
O concerns the phase space integral of 2 3 process. We explain how to use the 
Fortran routine BASES [37] to do numerical integration. Appendix [D] gives useful 
mathematical formulae related to loop integrals. In appendix [E] we explain the analyt- 
ical calculation of scalar one-loop four-point integrals with complex internal masses. 
The restriction is that at least two external momenta are lightlike. 



Chapter 1 

The Standard Model and beyond 



The Glashow-Salam- Weinberg (GSW) model of the electroweak interaction was pro- 
posed by Glashow [T], Weinberg [2] and Salam [3] for leptons and extended to the 
hadronic degrees of freedom by Glashow, Iliopoulos and Maiani [38]. The GSW model 
is a Yang-Mills theory [32] based on the symmetry group SU {2)l x f/(l)y. It describes 
the electromagnetic and weak interactions of the known 6 leptons and 6 quarks. The 
electromagnetic interaction is mediated by a massless gauge boson, the photon (7). 
The short-range weak interaction is carried by 2 massive gauge bosons, Z and W. 
The strong interaction, mediated by the massless gluon, is also a Yang-Mills the- 
ory based on the gauge group SU{3)c- This is known as Quantum chromodynamics 
(abbreviated as QCD) [U El El [7]. The Standard Model of particle physics is just 
a trivial combination of GSW model and QCD. The particle content of the SM is 
listed in Table. II. 1[ There is an additional scalar field called the Higgs boson (H), 
the only remnant of the spontaneous symmetry breaking (SSB) mechanism invented 
by Brout, Englert, Guralnik, Hagen, Higgs and Kibble (iQl iH SSI iSl E] . The SSB 
mechanism is responsible for explaining the mass spectrum of the SM. 

To date, almost all experimental tests of the three forces described by the Standard 
Model agree with its predictions [HI El US]- The measurements of Mw and Mz together 
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Table 1.1: Particle content of the standard model 





Particles 


Spin 


Electric charge 


Leptons 


(e,/i,r) 


1/2 


-1 






1/2 





Quarks 


{u,c,t) 


1/2 


2/3 




{d,s,b) 


1/2 


-1/3 


Gauge bosons 


gluon (g) 


1 







il,Z) 


1 









1 


±1 


Higgs 


H 









with the fact that their relation = M|c^ (with ^ 0.77 defined in Eq. 
( ll.lOp ) has been experimentally proven imply two things. First, the existence of 
massive gauge bosons means that the local gauge symmetry is broken. Second, the 
mass relation indicates that the effective Higgs (be it fundamental or composite) is 
isospin doublet [l5]. Experiments have also confirmed that couplings that are mass- 
independent like the ones of quarks and leptons to the and Z gauge bosons or 
triple couplings among electroweak gauge bosons agree with those described by the 
gauge symmetry [15]. It means that the only sector which remains untested is the 
mass couplings or in other words the nature of SSB mechanism. 

The primary goal of the LHC is to find the scalar Higgs boson and to understand 
its properties. The main drawback here is that we do not know the value of the 
Higgs mass which uniquely defines the Higgs profile. The LEP direct searches for the 
Higgs and precision electroweak measurements lead to the conclusion that 114GeV < 
Mh < 190GeV [9]. The most prominent property of the Higgs is that it couples mainly 
to heavy particles at tree level. This has two consequences at the LHC: the Higgs 
production cross section is small and the Higgs decay product is very complicated and 
usually suffers from huge QCD background. Thus, it is completely understandable 
that searching for the Higgs is not an easy task, even at the LHC. 



1.1. QCD 
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1.1 QCD 

The classical QCD Lagrangian reads 



jCqcd = i>{ilp-m),p - ^ Tr F^,F''^ (1.1) 



where 



with a = 1,...,8; iJj is a fermion field belonging to the triplet representation of 
SU{3)c' group; A the gauge boson field and Qs is the strong coupling; Ta are Gell- 
Mann generators. The corresponding Feynman rules in the 't Hooft-Feynman gauge 
read: 

i k j -Sjj 

^ ^ — m + ie 



a k b 
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+ 9mf''{9ap9,s-9aS9p,) 
+ 9'r''f'''{9af39,s-9.,9(ss) 



We have adopted the Feynman rules of [l6l HTj (derived by using C) which differ 
from the normal Feynman rules (derived by using iC) by a factor i. One can use 
those Feynman rules to calculate tree-level QCD processes or QED-like processes by 
keeping in mind that the gluon has only two transverse polarisation components. 
However, in a general situation where a loop calculation is involved one needs to 
quantize the classical Lagrangian (11.11) . The covariant quantization following the 
Faddeev-Popov method [48j introduces unphysical scalar Faddeev-Popov ghosts with 
additional Feynman rules: 



Jab 



The main difference between QCD and QED is that the gluon couples to itself while 
the photon does not. In QED, only the transverse photon can couple to the electron 
hence the unphysical components (longitudinal and scalar polarisations) decouple 
from the theory and the Faddeev-Popov ghosts do not appear. The same thing hap- 
pens for the gluon-quark coupling. However, an external transverse gluon can couple 
to its unphysical states via its triple and quartic self couplings. Those unphysical 
states, in some situation, can propagate as internal particles without coupling to any 
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quarks and give an unphysical contribution to the final result. In that situation, one 
has to take into account also the ghost contribution for compensation. 

Indeed, there is another way to calculate QCD processes by taking into account 
only the physical contribution, i. e. only the transverse gluon components involve and 
no ghosts appear. This is called the axial (non-covariant) gauge [IH]. The main differ- 
ence compared to the above covariant gauge is with the form of the gluonpropagator. 
The covariant propagator includes the unphysical polarisation states visL. 

2 

i=l 

where are two unphysical polarisation states and with z = 1,2 are the two 
transverse polarisation states. In the axial gauge, the gluon propagator takes the 
form 

p _ Sab \^ 



/c^ + ie 

i=l 



6, 



ab 



k"^ + ie 



-9fiu + - 



n.k 



;i-4) 



with n"^ = and n.k ^ 0, which includes only the transverse polarisation states. The 
main drawback of this axial gauge is that the propagator's numerator becomes very 
complicated. 

The main calculation of this thesis is to compute the one-loop electroweak cor- 
rections to the process gg — » bbH. Though the triple gluon coupling does appear 
in various Feynman diagrams, it always couples to a fermion line hence the virtu- 
ally unphysical polarisation states cannot contribute and the ghosts do not show up. 
We will therefore use the covariant Feynman rules and take into account only the 
contribution of the transverse polarisation states of the initial gluonj^. 



iSee p.511 of [KD]. 

^If one follows the traditional amplitude squared method and wants to use the polarisation sum 
identity ^ e^e,y = —9^11 then one has to consider the Feynman diagrams with two ghosts in the 
initial state. 
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1.2 The Glashow-Salam- Weinberg Model 

The classical Lagrangian of the GSW model is composed of a gauge, a Higgs, a fermion 
and a Yukawa part If 

Cc = 'Cg + Ch + Cp + Cy- (1.5) 
Each of them is separately gauge invariant and specified as follows: 

1.2.1 Gauge sector 

The Lagrangian of the gauge part of the group SU{2)l x U{1)y reads 

Cg = -lid.W,^ - d,W^ + ge^'^WlWlf - \{d^B, - d,B,)\ (1.6) 

where a,b,c E {1, 2, 3}, W^^ are the 3 gauge fields of the SU (2) group, 5^ is the U (1) 
gauge field, the SU{2) gauge coupling g, the U{1) gauge coupling g' and e"'^'^ are the 
totally antisymmetric structure constants of SU{2). The covariant derivative is given 
by 

D^ = d,-igT-W;-ig'YB,, (1.7) 

where T'^ = cr"/2 with a" are the usual Pauli matrices, the hypercharge according to 
the Gell-Mann Nishijima relation 

Q = T^ + Y. (1.8) 
The physical fields W^, Z, A relate to the W"" and B fields as 

= cwW^^ - swWf, (1.9) 
= swWl + cwW% 

■^For more technical details of the GSW model, its one-loop renormalisation prescription and 
Feynman rules, we refer to [STl l46l l47]. 
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with 



V ^ + f V f + 9 



the electromagnetic couphng e 

99' e , e 

e = / , , =, 9 = —, 9=—- (1-11) 
V 5" + 9 

1.2.2 Fermionic gauge sector 

Left-handed fermions L of each generation belong to SU{2)l doublets while right- 
handed fermions it! are in SU (2) l singlets. The fermionic gauge Lagrangian is just 

U^iY, L^D^L + I R^D^R, (1-12) 

where the sum is assumed over all doublets and singlets of the three generations. 
Note that in the covariant derivative D/^ acting on right-handed fermions the term 
involving g is absent since they are SU (2) ^ singlets. Neutrinos are left-handed in the 
SM. Fermionic mass terms are forbidden by gauge invariance. They are introduced 
through the interaction with the scalar Higgs doublet. 



1.2.3 Higgs sector 

Mass terms for both the gauge bosons and fermions are generated in a gauge invariant 
way through the Higgs mechanism. To that effect one introduces minimally a complex 
scalar SU{2) doublet field with hypercharge F = 1/2 

$=('^^1 = 1 ) , (0 I $ I 0) = v/V2, (1.13) 

where the electrically neutral component has been given a non-zero vacuum expec- 
tation value V to break spontaneously the gauge symmetry SU{2)l x U{1)y down to 
U{1)q. The scalar Lagrangian writes 

jCh = {Di.^yiD^'^) - V{^), V{^) = + A($^*)l (1.14) 
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After rewriting in terms of , Xs > H and imposing the minimum condition on the 
potential V{^) one sees that ^ind X3 ^i-re massless while the Higgs boson obtains 
a mass 



M 



H 



I? = Xv^. 



;i.i5) 



X , Xs are called the Nambu-Goldstone bosons. They are unphysical degrees of 



freedom and get absorbed by the and Z to give the latter masses given by 

ev ^ ^ ev 



w 



2s 



w 



(1.16) 



1.2.4 Fermionic scalar sector 



Fermion masses require the introduction of Yukawa interactions of fermions and the 
scalar Higgs doublet 

fuD^ 

down 



up 



;i.i7) 



where fjjj^ with i,j G {1,2,3} the generation indices are Yukawa couplings, $ = 
zcr2$*. Neutrinos, which are only right-handed, do not couple to the Higgs boson and 
thus are massless in the SM. The diagonalization of the fermion mass matrices 
introduces a matrix into the quark- W-boson couplings, the unitary quark mixing 
matrix [H] 



^ Vud Vus Vuh \ ( 0.97383 0.2272 0.00396 \ 



V 



Vcd Ks Kfe 

V Vtd Vts Vtb I 



1.18) 



0.2271 0.97296 0.04221 
y 0.00814 0.04161 0.9991 / 

which is well-known as Cabibbo-Kobayashi-Maskawa (CKM) matrix. There is no 
corresponding matrix in the lepton sector as the neutrinos are massless in the SM. 

For later reference, we define = ^J^m^jv where mj is the physical mass of a 
fermion. 
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1.2.5 Quantisation: Gauge-fixing and Ghost Lagrangian 

The classical Lagrangian jCq has gauge freedom. A Lorentz invariant quantisation 
requires a gauge fixing (otherwise the propagators of gauge fields are not well-defined) . 
The 't Hooft hnear gauge fixing terms read 

This leads to a gauge fixing Lagrangian 

^fix = —[{Fy + {F^r + 2F+F-]. (1.20) 

Cfix involves the unphysical components of the gauge fields, i.e. field components 
with negative norm, which lead to a serious problem that the theory is not gauge 
invariant and violates unitarity. In order to compensate their effects one introduces 
Faddeev Popov ghosts u°'{x), u°'{x) {a — A, Z, W^) with the Lagrangian 

C,nost = u^{x)-^^u\x), (1.21) 

where -^j^ is the variation of the gauge fixing operators under infinitesimal 
gauge transformation parameter 9^{x). An element of the SU(2)l x f/(l)y group 
has a typical form G = Q-^9T°'^oL{x)-ig'YeY{x) ^ Faddeev Popov ghosts are scalar fields 
following anticommutation rules and belonging to the adjoint representation of the 
gauge group. 

In a practical calculation, the final result does not depend on gauge parameters. 
Thus one can choose for these parameters some special values to make the calculation 
simpler. For tree-level calculations, one can think of the unitary gauge = S,'^ — 
oo where the Nambu-Goldstone bosons and ghosts do not appear and the number 
of Feynman diagrams is minimized. For general one-loop calculations, it is more 
convenient to use the 't Hooft Feynman gauge — — — 1 where the numerator 
structure is simplest. 



18 



Chapter 1. The Standard Model and beyond 



It is worth knowing that the 't Hooft hnear gauge fixing terms defined in Eq. 
fll.lQp can be generahsed to include non-hnear terms as follows [521 HZ! 



d^Z, + Mz^'zX3 + ^ezeHx3 



ii^X-i)x ; (1-22) 

with the gauge fixing term for the photon remains unchanged. It is simplest 
to choose ^'zw ~ ^z,w- Those non-linear fixing terms involve five extra arbitrary 
parameters C = (a, /3, 5, e). The advantage of this non-linear gauge is twofold. First, 
in an automatic calculation involving a lot of Feynman diagrams one can perform 
the gauge-parameter independence checks to find bugs. Second, for some specific 
calculations involving gauge and scalar fields one can kill some triple and quartic 
vertices by judiciously choosing some of those gauge parameters and thus reduce 
the number of Feynman diagrams. This is based on the fact that the new gauge 
parameters modify some vertices involving the gauge, scalar and ghost sector and 
at the same time introduce new quartic vertices 07]. In the most general case, the 
Feynman rules with non-linear gauge are much more complicated than those with 't 
Hooft linear gauge, however. 

With Cfix and C ghost the complete renormalisable Lagrangian of the GSW model 
reads 

C-GSW = ^^C + ^^fix + ^ghost- (1.23) 



1.2.6 One-loop renormalisation 

Given the full Lagrangian Cgsw above, one proceeds to calculate the cross sec- 
tion of some physical process. In the framework of perturbative theory this can 
be done order by order. At tree level, the cross section is a function of a set of 
input parameters which appear in Cgsw- These parameters can be chosen to be 
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O = {e, Miv, Mz, Mh, ^} which have to be determined experimentally. There 
are direct relations between these parameters and physical observables at tree level. 
However, these direct relations are destroyed when one considers loop corrections. 
Let us look at the case of Mw as an example. The tree- level W mass is directly 
related to the Fermi constant through 

slM^^ = (1.24) 
When one takes into account higher order corrections, this becomes [531 EU [55] 

where Ar containing all loop effect is a complicated function of Myy and other input 
parameters. A question arises naturally, how to calculate Ar or some cross section 
at one-loop level? The answer is the following. If we just use the Lagrangian given 
in Eq. f ll.23p . follow the corresponding Feynman rules to calculate all the relevant 
one-loop Feynman diagrams then we will end up with something infinite. This is 
because there are a lot of one- loop diagrams being UV-divergent. This problem 
can be solved if Cqsw is renormalisable. The renormalisability of nonabelian gauge 
theories with spontaneous symmetry breaking and thus the GSW model was proven 
by 't Hooft [Sni EZ|- The idea of renormalisation is that we have to get rid of all UV- 
divergence terms originating from one-loop diagrams by redefining a finite number of 
fundamental input parameters O in the original Lagrangian Cqsw- This is done as 
follows 

e ^ {l + 6Y)e, 
M M + 6M, 

ip {l + 5Z^/'^)%l). (1.26) 

The latter is called wave function renormalisation. The renormalisation constants 5Y , 
5M and 5Z^/'^ are fixed by using renormalisation conditions to be discussed later. The 
one-loop renormalised Lagrangian writes 

^^jsw' = ^Gsw + ^C,GSw- (1-27) 
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The parameters O in jC^T^^^ are now called the renormalised parameters determined 
from experiments. From this renormalised Lagrangian one can write down the cor- 
responding Feynman rules and use them to calculate Ar or any cross section at one 
loop. The results are guaranteed to be finite by 't Hooft. 

We now discuss the renormalisation conditions which define a renormalisation 
scheme. In this thesis, we stick with the on-shell scheme where all renormalisation 
conditions are formulated on mass shell external fields. To fix 6Y, one imposes a 
condition on the e~^e~A vertex as in QED. The condition reads 

{e^e~A one-loop term + e'''e~y4 counterterm) |g=Qp|=.^2= 0, (1-28) 

where q is the photon momentum, p± are the momenta of respectively. All 5Ms 
are fixed by the requirement that the corresponding renormalised mass parameter is 
equal to the physical mass which is the single pole of the two-point Green function. 
This translates into the condition that the real part of the inverse of the corresponding 
propagator is zero. 5Z^^'^s are found by requiring that the residue of the propagator 
at the pole is 1. To be explicit we look at the cases of Higgs boson, fermions and 
gauge bosons, which will be useful for our main calculation of pp bbH. The Higgs 
one-particle irreducible two-point function is 11^ (g^) with q the Higgs momentum. 
One calculates this function by using Eq. (11.271) 

where the counterterm contribution is denoted by a caret, the full contribution is 
denoted by a tilde. The two renormalisation conditions read 

Ren^(M^)=0, -^Ren^(g^) ^ ^=0. (1.30) 



This gives 



5Z;/' = _l^Ren^(g2) . (1.31) 
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For a fermion with ijj = ijj ^ + ijj r {ip^R = Pl,r'^ with Pl^r = i(l ip 75), respectively), 
the one-particle irreducible two-point function takes the form 



S(g^) + S(g^), 

Ki + K^i + i^57^75, 



(1.32) 



with 



k^ 



K, 



57 



-{6Z 



bZf) 



6m 



a/2 



(1,33) 



The two renormalisation conditions become 



■mfReK^{mj) + ReKi{m'j) = and Re i^'5^(mj) = 



Re 



f]=mf 



This gives 



5m/ = Re{mfK^{m^j) + Ki{m'^j)^, 



5Z]I^ 

JL 



5Zl^ 



~ Re (^Kr^^{mj) + K^{m 



ruf 



(1.34) 



rrif 



dq^ 



mfReK^{q^) + ReKi{q^) 



mfReK^{q^) + ReKi{q^) 



q^=mf 



For gauge bosons, the one-particle irreducible two-point functions write 



(1.35) 



KM') 

KM") 



)K{<f) + ^nr(g^) 



qfj.qu 2 
qf^qu. 



q^ 



, ^/^^-nV(g2)^ 



n 



V 



6M^ + 2(M^ - q'^)6Z, 



a/2 



flY = 5Ml + 2MUZ. 



1/2 



(1.36) 



*For massless gauge bosons like the photon, the longitudinal part vanishes. 
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where V = W, Z. We do not touch the photoio since it is irrelevant to the calculations 
in this thesis, which are only related to the Yukawa sector. It is sufficient to impose the 
two renormalisation conditions (for the pole-position and residue) on the transverse 
part nj^(g^) to determine 6 My and SZy"^. The longitudinal part is automatically 
renormalised when the transverse part is, if the theory is renormalisable. The two 
conditions write 



Ren^(M2) = 0, ^Ren^(g2) 



0, 



;i.37) 



which give 



5M^ 



Ren^(M^), SZ. 



V 



1 d 
2d^ 



Ren^(g^) 



;i.38) 



In practical calculations, one has to calculate n^(g^), i^i(g^), K^{q^), K^^{q^) and 
n^(g^) as sums of various two-point functions. The full results in the SM can be 
found in |l^lll71EI]• 



1.3 Higgs Feynman Rules 

In order to understand the phenomenology of Higgs production, it is important to 
write down the relevant Feynman rules. 



The Feynman rules listed here are taken from [47]. Their Feynman rules derived 
from Casw differs from the normal Feynman rules derived by using ijCqs^y by a 
factor cl- A particle at the endpoint enters the vertex. For instance, if a line is 
denoted as W~^, then the line shows either the incoming or the outgoing W~ . 
The momentum assigned to a particle is defined as inward. The following Feynman 
rules are for the linear gauge. 



^In a general case, one should keep in mind that there is mixing between the photon and the Z 
boson. 

^In the QCD section [Ol we have adapted the same rules of this section. 



1.3. Higgs Feynman Rules 
Propagators 



z 

f 









1 




-M| 




-1 


^- 


-ruf 


-1 



IW 



k/jki^ 



k^ - iwM^ 



H 



± 



X 



X3 



k^-Mjj 
-1 

k-" - iwM^ 



w 



-1 



k^-izMl 



Vector- Vector-Scalar 
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Scaleir-Scaleir- Vector 



Pi Pi P2 P3 il^) 



\ 

\ 

\ 



9\f\r>^f^y\r\j Phi'' 
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/ 

/ 



P2 
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Pi P2 P3 



H H H -e K^Ml 



s 

\ 



/ 
/ 

/ 

/ 



y 



s 

s 

P3 ^ Ml 

H X- X+ -e ^ 5. 
/ 2swMw 



Pi A/r2 

H X3 X3 -e 



2swMw 



1.4. Problems of the Standard Model 
Fermion-Fermion-Scalar 
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Pi 



Pi P2 P3 


f f H - 


1 rrif 

e — 

2sw 




U/D U/D X3 (- 


-/+he 7^ 

2swMw 


U D 
■ 1 1 w 


- mu) + 




D U X- 
■ 1 1 u 


- mo) + 


[mu + ^0)75] 



We would like to make some connections between the underlying Feynman rules 
of the SM and the main calculation of this thesis, one-loop Yukawa corrections to the 
process gg bbH. The relevant vertices will be "scalar-scalar-scalar" and "fermion- 
fermion-scalar" . Of these, the vertex (bbxs) will be excluded as it will result in 
Feynman diagrams proportional to Xlf^^^ which are neglected in our calculation. 



1.4 Problems of the Standard Model 



In spite of its great experimenta 



success, the SM suffers from a conceptual problem 
known as the hierarchy problem |3. This problem is related to the quantum corrections 
to the Higgs mass. In the calculation of one-loop corrections to the Higgs mass, we 
see that quadratic divergences appear. Of course, these UV-divergences have to be 
canceled by the corresponding counter terms. The leading correction is proportional 
to the largest mass squared, assumed to be mf. Since the value of nit ~ 174GeV 



'''Indeed, there are other conceptual as well as phenomenological problems of the SM such as those 
related to gravity and dark matter. These discussions can be found in the recent review of Altarelli 
[45j and references therein. The discussion on the hierarchy problem can be found also in [581 159j. 
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is not so large, this correction is well under control in the SM. However, the SM 
is just an effective theory of a more general theory with heavy particles at some 
high energy scale, say the GUT scale Aqut ~ 2 x lO^^GeV where the three gauge 
couphng constants unify. The masses of those heavy particles are at the order of 
h-GUT- Those heavy particles must couple to the SM Higgs boson and hence give 
enormous corrections to Mh- The fact the Mh/ ^gut ~ lO^"*^^ means that an extreme 
cancelation occurs among those huge corrections. This is known as the naturalness or 
fine-tuning problem. A related question, called the hierarchy problem, is why Agut 3> 
Mz- These problems can be solved if there is a symmetry to explain that cancelation. 
There are a few options for such a symmetry, among them supersymmetry is the most 
promising candidate. 



1.5 Minimal Super symmetric Standard Model 

The Minimal Supersymmetric Standard Model is a theory describing the interactions 
of all SM fundamental particles, their superpartners and some additional Higgs parti- 
cles. None of these superpartners and new Higgs bosons has been seen in experiment. 
The fundamental superpartners arise as a consequence of the so-called supersymmetry 
(SUSY) imposed on the Lagrangian of the theory. The SUSY generator Q transforms 
a fermion into a boson and vice versa: 

Qjfermion) = |boson), Q|boson) = |fermion). (1.39) 

It means that each SM particle has a corresponding superpartner. The superpartners 
of a fermion, a vector gauge boson, a scalar Higgs boson are called a sfermion, a 
gaugino, a higgsino respectively. SUSY requires that a superpartner has the same 
quantum numbers as its corresponding particle except for the spin. Sfermions are 
scalar while gauginos and higgsinos have spin 1/2. One notices immediately that 
some mixings among gauginos and higgsinos are allowed. The MSSM Lagrangian has 
three symmetries: Lorentz symmetry, SM gauge symmetry and SUSY. 
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The fact that we have never observed a fundamental scalar selectron with the 
same mass as the electron means that SUSY is broken. To date there is no com- 
pletely satisfactory dynamical way to break SUSY. In the MSSM, SUSY is broken by 
introducing extra terms that explicitly break SUSY into the Lagrangian [60]. They 
are called soft-SUSY-breaking terms, all contained in Csoft- The purpose of Csoft is 
to give (quite heavy) masses to superpartners [591 EOl EH [62] . 

The SM particles obtain masses by the Higgs mechanism. In the SM, we just 
need one Higgs doublet $ (and $ = ia2^*) to generate masses for down quarks (up 
quarks). However, the same trick cannot be used for the MSSM since it will break 
SUSY. Thus one needs two complex Higgs doublets with opposite hypercharges 



H^={ M with Yh, = -1/2; H, = { M with Yh, = 1/2, (1.40) 



to give masses for down fermions and up fermions respectively. Before symmetry 
breaking, these two Higgs doublets have 8 independent real fields. After symmetry 
breaking, 3 vector gauge bosons Z, get masses by "eating" 3 Goldstone bosons, 
so five real fields remain. The MSSM therefore predicts the existence of 3 neutral 
Higgs bosons denoted H, h, A and 2 charged Higgs bosons denoted H^. 

In the unconstrained MSSM, Csoft introduces a huge number (105) of unknown 
parameters {e.g. intergenerational mixing, complex phases), in addition to 19 parame- 
ters of the SM [SniES]- This makes the phenomenology study of the MSSM extremely 
difficult if not impossible. There exists however the so-called contrained MSSMs with 
only a handful of parameters. Among them, mSUGRA is most well-known with the 
5 following parameters [60, EU ES [Ml [65] 



This is achieved by imposing some conditions on the soft-SUSY-breaking parameters. 
These parameters are required to be real and satisfy a set of boundary conditions at 
the GUT scale {Aqut ~ 2 x lO^^GeV) where the three gauge coupling constants unify. 





tan/5, mi/2, tuq, Aq, sign{fi). 



(1.41) 
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These boundary conditions say that: all gauginos have the same masses {mi/2), all 
sfermions and Higgs bosons have the same mass (mo) and all trilinear couplings in 
Csoft are equal at the GUT scale. 



1.5.1 The Higgs sector of the MSSM 

The scalar Higgs potential Vh comes from three different sources [591 EHl EH] : 



Vh = Vd + Vf + Vs 



softi 



J2 



vf = ^^'m\' + \H2\'), 

Vsoft = ml^HlH, + ml^HlH2 + Bij{H2.H, + h.c.), (1.42) 

where g, g' are the usual two couplings of the groups SU{2) and U{1) respectively; 
/i and are bilinear couplings; = + l-f^fP ^"^^ same definition for 

|if2p- The first two terms of Vh are the so-called D- and F- terms. The last term 
Vsoft is just a part of Cgoft discussed above. The MSSM Higgs potential contains the 
gauge couplings while the SM one given in Eq. fll.l4p does not. 

The neutral components of the two Higgs fields develop vacuum expectations 
values 

{Hi) = ^ , {H',) = ^. (1.43) 



One defines 



Comparing to the SM we have 



tan/? = — . (1.44) 

Vl 



v{ + vi = v\ (1.45) 



We now develop the two doublet complex scalar fields Hi and H2 around the vacuum, 
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into real and imaginary parts 

Hi = (H^, Hi 
H2 = 



V2 
1 



{vi + H^,+tP^ , H{) 



) = ^ (^2+ , V2 + Hi + iP^) 



(1.46) 



where the real parts correspond to the CP-even Higgs bosons and the imaginary parts 
corresponds to the CP-odd Higgs and the Goldstone bosons. By looking at the Eq. 
I1.42[ we see that those fields mix. After diagonalizing the mass matrices, one gets 

cos (3 sin (3 

— sin [3 cos [3 

cos [3 sin [3 

— sin (3 cos (3 




cos a sm a 
— sin a cos a 



) 










( Ht\ 


) 






ii) 


) 





1.47) 



with the mixing angle a given by 



cos 2a 



cos 2/3 



Ml 



Ml 



, sin 2a 



M 



H 



Mr 



1.48) 



where xs? ^^^6 massless Goldstone bosons to be eaten by the Z, respectively; 
A, H^, H and h are five physically massive Higgs bosons; the tree-level masses are 
given by 

25/i 



Ml 
MIh 



-5/i(tan/5 + cot/3) 

Ml + Ml,, 
1 
2 



sin 2/3' 



Ml + M| T ^(Ml + M|)2 - 4MlM| cos2 2/3 . 
We remark that < Mz at tree level. From Eqs. fll.48p and (11.491) we get 

cos^(/3-a)-^^^(^^-^^^) 



1.49) 



Ml{Ml - Ml 



2\ ' 



1.50) 



which will be useful later. 
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1.5.2 Higgs couplings to gauge bosons and heavy quarks 

Like in the SM, the Higgs boson couphngs to the gauge bosons are obtained from the 
kinetic terms with covariant derivatives of the Higgs fields Hi and H2. The Yukawa 
Higgs boson couphngs to the fermions are obtained from the Yukawa Lagrangian. We 
hst here some relevant couplings needed in this thesis. For a full account of Higgs 
couplings in the MSSM, we refer to [591 EH EH]- With Xwwh, ^zzh, ^bbH, ^uh are 
the SM couplings and using the same Feynman rules as in section 11.31 we have: 
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b 
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— XbbH 


[cos(/3 — 


a] 


) — tan /? sin(/? - 


■a)] 
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h 


— XttH\ 


sin(/3 — 


a) 


+ cot (3 cos(/5 — 


a)] 


b 


b 


h 


— XbbH 


[sin(/? - 


a) 


— tan 13 cos(/? - 


-a)] 



We remark that the bb{tt) coupling of either the H or h boson is enhanced (suppressed) 
by a factor tan (3 with the enhancement (suppression) magnitude depending on the 
value of sin(/3 — a) or cos(/3 — a). Thus one can have very large value of bottom-Higgs 
Yukawa coupling, leading to a large cross section if tan (3 is large. In the decoupling 
limit where Ma — > 00, i.e. cos(/3 — a) (see Eq. f ll.SOp ). the h is SM-like (the same 
couplings) while the H Yukawa coupling to bb{tt) is exactly enhanced(suppressed) by 
a factor tan (3. 
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Standard Model Higgs production 
at the LHC 



2.1 The Large Hadron Collider 

The Large Hadron Collider (LHC) is the world largest particle accelerator to date 
[10] . It collides two proton beams with the center-of-mass energy up to 14TeV. It 
is expected to start this year. It has four main experiments: ATLAS, CMS, LHCb 
and ALICE. ATLAS and CMS are general-purpose detectors. Their goals are to find 
the Higgs boson and discover new physics expected to be Supersymmetry. LHCb 
is for B-physics and CP violation. ALICE aim is to study the physics of strongly 
interacting matter at extreme energy densities, where the formation of a new phase 
of matter, the quark-gluon plasma, is expected. The number of events per second 
generated in the LHC collisions given by 

event -^'^eiient; (^■-^) 

where advent is the cross section for the event under study and L the machine lumi- 
nosity. The machine luminosity depends only on the beam parameters and can be 
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Table 2.1: LHC beam parameters relevant for the peak luminosity 



Number of particles per bunch (A*";,) 


1.15 X 10^^ 


Number of bunches per beam (n^) 


2808 


Revolution frequency (frev) 


11245Hz 


Relativistic gamma (7^) 


7461 {E = 7TeV) 


Normalized transverse emittance (e^) 


3.75 X 10-^cm 


Full crossing angle at the IP (6^) for ATLAS/CMS 


285/irad = 0.0163° 


RMS bunch length (cr^) 


7.55cm 


Transverse RMS beam size (a*) at ATLAS/CMS 


16.7/im 


Geometric luminosity reduction factor (F) at ATLAS / CMS 


0.84 


Optical beta function at ATLAS/CMS {(3*) 


55cm 



written for a Gaussian beam distribution as 

N^Ubfrevl: 



(2.2) 



where Nh is the number of particles per bunch, rii, number of bunches per beam, 
frev the revolution frequency, 7^ the relativistic gamma factor, the normalized 
transverse beam emittancfl, /3* the optical beta function at the collision point|^ and F 
the geometric luminosity reduction factor due to the crossing angle at the interaction 
point (IP): 



I/W1 + 



2a* 



(2.3) 



where 6c is the full crossing angle at the IP, the root-mean-square (RMS) bunch 
length and a* the transverse RMS beam size at the IP. Note that F < 1 since the angle 
between two beams at the collision point is greater than zero. The above expressions 
assumes equal beam parameters for both circulating beams. In order to calculate the 



^The beam emittance (e), units of length, is the extent occupied by the particles of the beam in 
phase space. A low emittance particle beam is a beam where the particles are confined to a small 
distance and have nearly the same momentum. The normalized beam emittance e„ = ^rPr^- 

^The optical beta function /3(s) appears in the amplitude of the solution of the equation of a 
harmonic ossilator x" [s) + Kx{s) — describing the transverse beam dynamics (linear force in x 
and y, s is the beam direction). 
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peak luminosity, we need the LHC beam parameters given in Table I2.1[ We then get 
the value of the peak luminosity at ATLAS/CMS L = 10^^cm~^s~^. Note that in 
luminosity formula (12.21) there are only two beam parameters (3* and F which depend 
on the collision point (experiment position), other beam parameters are the same for 
LHCb and ALICE. Other experiments LHCb and ALICE with Pb — Pb collision have 
lower luminosity, lO^^cm^^s^^ and 10^''cm"^s~^ respectively. If the Higgs production 
cross section is Ipb then one has 10~^ events per second at ATLAS/CMS. 

The luminosity in the LHC is not constant over a physics run but decays due 
to the degradation of intensities and emittances of the circulating beams. The LHC 
integrated luminosity {Lint = Jq L{t)dt) expected per year is between 80fb~^ and 
120fb~^ [69]. 

Despite of a number of enormous advantages: the high values of luminosity and 
center-of-mass energy, the dedicated detectors ATLAS and CMS designed to discover 
the Higgs boson and Supersymmetry particles; to get something out of a huge amount 
of data produced by the LHC is not easy. The main challenge for LHC physics is at 
the way one analyses data. Let us make it clear. The total cross section for inelastic 
scatterings which can be seen by LHC detectors is about 60mb [70]. It is translated 
to 6 X 10^ events per second (event rate). To be optimistic and forget about the 
background, we say that the typical total cross section for Higgs production at the 
LHC is 60pb which corresponds to 6 x 10^^ event rate. Thus we have to know how to 
find out 6 x 10~^ events in 6 x 10^ events per second. The task becomes much more 
complicated when backgrounds are taken into account. The only way to distinguish 
signal from backgrounds is to study various distributions. 



2.2 SM Higgs production at the LHC 

The SM Higgs profile is: 
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Figure 2.1: Four mechanisms to produce the Higgs at the LHC. The right figure shows 
the cross sections as functions of Mh, which include the full NLO QCD corrections 



(1) Electric charg 


e: neutral 


(5) Mass: 114GeV < Mh < 190GeV ^ 


(2) Color charge: 
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massive gauge boson, / a fermion (see 


(4) CP: even 




section ll.3|) . 



In this section, we assume that the above Higgs profile is correct. The last property 
means that the Higgs boson couples mainly to heavy particles: W and Z gauge bosons, 
the top quark and, to lesser extent, the bottom quark. The four main mechanisms 
for single Higgs boson production are [TT]: 



1. gluon-gluon fusion: gg H 

2. associated production with heavy quarks: gg — > QQ + H with Q = b,t 

3. vector boson fusion: qq —>■ V*V* qq + H with g is a light quark 

4. associated production with W{Z): qq{qq) — > W{Z) + H 
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We will use the item number to refer to the corresponding process. The Feynman 
diagrams and total cross sections are shown in Fig. 12.11 It is important to know that 
in the second mechanism, there is also a contribution from qq in the initial state. 
However, this contribution is very small compared to the gg contribution and can 
be neglected. There are two reasons for this. First, the quark density in the proton 
is very small compared to the gluon density. Second, qq contribution contains only 
S-channel diagrams strongly suppressed at high energy while gg contribution contains 
T and U channels. 

From the plot of total cross sections as functions of Mh in Fig. 12.11 we see that 
for Mh G [114, 190] GeV the total cross section for gg fusion process is largest and 
almost 10 times bigger than the second one from vector boson fusion process. For 
Mh = 120GeV, the total cross sections for the four processes are about: 



H ttH qqH WH{ZH) 
40pb 0.7pb 5pb 2(l)pb 



(2.4) 



We can conclude that: for Higgs production, the LHC is just the gg collider. The 
Higgs discovery channel is clearly the gg fusion process. 



The first important signal for a Higgs boson is a peak in the invariant mass 
distribution of its leptonic decay products [f. However, if we observe a new bump 
at around 120GeV in some production channel it is not necessarily the SM Higgs 
boson. All we can say is that it is new physics, a new particle with mass around 
120GeV. Now we look back at the Higgs profile and assume that all the properties 
are confirmed except for the last one (the coupling property). We still cannot say 
that it is the SM Higgs boson. Remember that the Higgs originates from the Higgs 
potential defined by its self coupling (A) and interacts with vector gauge bosons and 
heavy fermions in some special ways. The best way to check the SSB mechanism 
in the SM is to determine the three independent Higgs couplings: Xvvh^ ^ffH and 
^HHH- If all these three couplings are consistent with the SM prediction then we can 



^The Higgs hadronic decay products suffer from extreme QCD jet backgrounds. 
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definitely say that it is the SM Higgs boson. However, measuring accurately those 
three Higgs couplings is an extremely difficult task if not impossible at the LHC. The 
Higgs self-coupling is the most difficult one. For this, one can think of gg — > HH but 
the total cross section is small and background is large. The two other couplings can 
be accessed at the LHC [TT], [7T] . 



2.3 Experimental signatures of the SM Higgs 

The SM Higgs boson is a heavy particle and decays. Its decay branching ratios 
(Fj/r//) and total width {Th) are shown in Fig. 12. 2[ We ffist look at the plot for the 
Higgs decay branching ratios (BR) and impose the following conditions: BR > 10^^, 
leptonic (photonic) decay product or heavy quarks in the decay product. These are 
the conditions for observing the Higgs in experiment. With the Higgs profile in hand 
{Mh e [114, 182]GeV) we are left with 7 branching ratios: bb, W+W', ZZ, t+t', 
cc, 77, Z'j. Among these, cc and Z7 can be discarded. For cc there are two reasons: 
BR{cc) is 10 times smaller than BR{bb) and it is more difficult to tag a charm-quark 
in experiment than a bottom-quark. For Z7 there are also two reasons: BR{Z'~f) is 
at the same order as BR{'~f'~f) and the Z is heavy and decays dominantly into hadrons 
and neutrinos and thus the combined branching ratio for this channel is very small. 
Thus we now have 5 potential Higgs signatures 

bb 77 T+T- W+W- ZZ 

0.68 2.16 X 10-3 6.78 x 10"^ 0.13 1.49 x 10"^ 

where the second row are branching ratios for Mh = 120GeV. Of these bb and 77 
can be observed "directly" in experiments. The three other decay modes cannot be 
directly seen in experiments and more branching ratios must be taken into account. 
77 is the most beautiful signal in experiment but its branching ratio is smallest. This 
is unlikely to be a discovery channel. The gg H ^ bb suffers from huge QCD 
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100 130 160 200 300 500 700 1000 

Mh [GeV] 

Figure 2.2: Upper: The SM Higgs boson decay branching ratios as a function of Mh- 
Lower: The SM Higgs boson total decay width as a function of Mh- Ref- fTlf - 
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backgrounds □ [72] and cannot be realised in LHC experiments [TT]. For bb signal, we 
have to use the mechanism of Higgs production associated with massive gauge bosons 
or heavy quarks. If we stick with the gg fusion mechanism then the Higgs discovery 
signal can be r+r", W^W^ or ZZ depending on the value of Mh- t^t^ signal is 
important for small Higgs mass. ZZ can give a beautiful 4-lepton signal but only for 
very large Higgs mass {Mh > 2Mz)- For the large range of Mh left, W^W~ decay 
mode is the best signal to pursue. 

The total Higgs decay width is shown in the Fig. 12.21 We first notice that Th < 
0.63GeV for Mh < 180GeV. The SM Higgs boson is a heavy long-lived particle. This 
makes it easy to determine the Higgs mass accurately by using if 77 for small 
Mh and H ^ ZZ ^ l+l'l+l' for large Mh ^WM- However, it will be difficult to 
measure Vh precisely. 



2.4 Summary and outlook 

If there exists only one Higgs boson in nature as predicted by the SM, it should be 
found at the LHC via the gluon-gluon fusion channel. However, if nature prefers 
having more than a unique Higgs boson, say five Higgs bosons as anticipated by the 
MSSM, then the situation of the SM-like Higgs (with the same Higgs profile given 
in section 12.21 except for the last coupling property) production at the LHC will be 
very different. From the results of subsection 11.5.21 we know that the SM-like Higgs 
coupling to the bottom quark can be strongly enhanced for large value of tan (3 while 
the same coupling with the top quark is much suppressed by the same factor. Thus 
the SM-like Higgs production associated with two bottom quarks can be the dominant 
mechanism at the LHC [59]. This channel with H W^W~ -signaXnie is also very 
valuable for the determination of the bottom-Higgs Yukawa coupling. In order to 
identify those new physics characters or just to confirm the SM Higgs properties, the 

''in a small window about the Higgs mass, the QCD background pp — > bb is still very large. 
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study of SM Higgs production associated with two bottom quarks at the LHC is very 
important. Indeed, it is the topic of the next chapter. 



Chapter 3 



Standard Model bbH production at 
the LHC 



The content of this chapter is based on our pubhcations [THl [77] . 

3.1 Motivation 

The study of the Higgs properties such as its self-couphngs and couphngs to the other 
particles of the standard model (SM) will be crucial in order to establish the nature of 
the scalar component of the model. In this respect most prominent couplings, in the 
SM, are the Higgs (Xhhh), the top (Xuh), and to a much lesser degree the bottom 
(Afcfe/f), Yukawa couplings. The top Yukawa coupling is after all of the order of the 
strong QCD coupling and plays a crucial role in a variety of Higgs related issues. 
In our main calculation of EW corrections to bbH, the leading contribution includes 
terms with largest powers of Xuh- 

The next-to-leading order (NLO) QCD correction to pp bbH has been cal- 
culated by different groups relying on different formalisms. In a nut-shell, in the 
five-flavour scheme (5FNS)[78l [79], use is made of the bottom distribution function 
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so that the process is approximated (at leading order, LO) by the fusion bb H. 
This gives an approximation to the inclusive cross section dominated by the untagged 
low pt outgoing b jets. If only one final b is tagged, the cross section is approximated 
by gb bH . The four flavour scheme (4FNS) has no b parton initiated process but 
is induced by gluon fusion gg bbH, with a very small contribution from the light 
quark initiated process qq — > 66-fil|. Here again the largest contribution is due to low 
Pt outgoing 6's which can be accounted for by gluon splitting into bb. The latter needs 
to be resummed and hence one recovers most of the 5FNS calculation while retaining 
the full kinematics of the reaction. QCD NLO corrections have been performed in 
both schemes [79l l80l ISTl [82] and one has now reached a quite good agreement [83]. 
The 5FNS approach, which at leading order is a two-to-one process has allowed 
the computation of the NNLO QCD correction [84^ [53] and very recently the elec- 
troweak/SUSY (supersymmetry) correction [86] to 66 — » 0, any of the neutral 
Higgs boson in the MSSM. SUSY QCD corrections have also been performed for 
gg —>■ 66/i[871 [HE] where h is the lightest Higgs in the MSSM as well as to gb bd\89\. 

In order to exploit this production mechanism to study the Higgs couplings to 
6's, one must identify the process and therefore one needs to tag both 6's, requiring 
somewhat large pr b. This reduces the cross section but gives much better signal over 
background ratio. For large pt outgoing quarks one needs to rely on the 4FNS to 
properly reproduce the hight pt b quarks. The aim of this and next chapters is to 
report on the calculation of the leading electroweak corrections to the exclusive bbH 
final state, meaning two 6's are detected. These leading electroweak corrections are 
triggered by top-charged Goldstone loops whereby, in effect, an external b quark turns 
into a top. This transition has a specific chiral structure whose dominant part is given 
by the top mass or, in terms of couplings, to the top Yukawa coupling. Considering 
that the latter is of the order of the QCD coupling constant, the corrections might 
be large. In fact, as we shall see, such type of transitions can trigger gg bbH even 

^In fact qq —^ bbH is dominated by qq HZ* bbH and does not vanish for vanishing bottom 
Yukawa couphng. However this contribution should be counted as ZH production and can be 
exchided by imposing an appropriate cut on the invariant mass of the bb pair. 
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with vanishing X^hH in which case the process is generated solely at one-loop. We will 
quantify the effect of such contributions. 

3.2 General considerations 

Before discussing the details of the calculation it is educative to expose some key 
features that appear when one considers the electroweak corrections at one-loop 
compared to the structure we have at tree-level or even the structure that emerges 
from QCD loop calculations. In particular the helicity structure is quite telling. So 
let us set our definition first. The process we consider is g{pi,Xi) + g{p2, ^2) 
b{p3, A3) + 6(p4, A4) + H{p^). Aj = ± with z = 1, 2, 3, 4 are the hehcities of the gluons, 
the bottom and anti-bottom while pi are the momenta of particles B The correspond- 
ing helicity amplitude will be denoted by A{\i, A2; A3, A4). 

3.2.1 Leading order considerations 



Figure 3.1: All the eight Feynman diagrams can he obtained by inserting the Higgs 
line to all possible positions in the bottom line. 

At tree-level, see Fig. 13.11 for the contributing diagrams, the Higgs can only at- 
tach to the 6-quark and therefore each diagram, and hence the total amplitude, is 

^The cross section for — > g* — > bbH (with input parameters given in section [5751 and q = u,d, s) 
is about 0.7% of the a{gg — > bbH) and thus can be totally neglected. There are two reasons for 
this. First the density of gluon at the LHC is much bigger the the density of quarks. Second the 
qq bbH which is S-channel is rather suppressed at high energy while the gg — > bbH contains also 
T and U channels which give the dominant contribution. 



g(Pi.'H) 




b(P3, h) 
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proportional to the Higgs coupling to bb, XbbH- Compared to the gluon coupling this 
scalar coupling breaks chirality. These features remain unchanged when we consider 
QCD corrections. Moreover the QCD coupling and the Higgs coupling are parity 
conserving which allows to relate the state with helicities (Ai, A2; A3, A4) to the one 
with (— Ai, — A2; — A3, — A4) therefore cutting by half the number of helicity ampli- 
tudes to calculate. With our conventions for the definition of the helicity states, see 
Appendix |Xl parity conservation for the tree-level helicity amplitude gives 

v4o("~Ai, — A2; —A3, — A4) = A3A4^o(Ai, A2; A3, A4)*. (3.1) 

This can be generalised at higher order in QCD with due care of possible absorptive 
parts in taking complex conjugation. 

The number of contributing helicity amplitudes can be reduced even further at 
the leading order, in fact halved again, in the limit where one neglects the mass of 
the 6-quark that originates from the 6-quark spinors and therefore from the b quark 
propagators. We should in this case consider the XbbH as an independent coupling, 
intimately related to the model of symmetry breaking. In this case chirality and 
helicity arguments are the same, the b and b must have opposite helicities for the 
leading order amplitudes and hence only ^o(Ai, A2; A, —A) remain non zero. In this 
limit, this means that only a string containing an even number of Dirac 7 matrices, 
which we will label in general as F'^^^" as opposed to r°'^'^ for a string with an odd 
number of 7's, can contribute. 

In the general case and reinstating the b mass, we may write the helicity amplitudes 

as 

^(Ai,A2;A3,A4) = u{\s){TZ:i + TrX)^i>^^) 

= 6x„-x, (^™'^" + mU°"") + 6x„x, (^°"" + m,^^^'^") .(3.2) 

The label in and are the contributions of T^^^"" to the amplitude 

and likewise for This way of writing shows that nib originates from the mass 
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insertion coming from the massive spinors and are responsible for chirality flip. In 
the limit rrib — > 0, T^'^'^^^ and r^^'^;^^ contribute to different independent helicity am- 
plitudes. In general T'^'"'^"' and r°'^'^ differ by a (fermion) mass insertion. In fact r°'^'^ 
is proportional to a fermion mass insertion from a propagator. At leading order the 
mass insertion is naturally nif,, such that r°'^'^ is 0{mf,). This shows that at leading 
order, corrections from = to the total cross section are of order 0{ml). Of 
course there might be some enhancement of the 0{ml) terms if one remembers that 
the cross section can bring about terms of order ml/{p^)'^. However, in our calcula- 
tion where we require the 6's to be observed hence requiring a cut, the effect will 
be minimal. With rrib = 4.62GeV, the effect of neglecting nib is that the cross section 
is increased by 3.7% for |p^*| > 20GeV and 1.1% for |p^^| > 50GeV. At one-loop, the 
chiral structure of the weak interaction and the contribution of the top change many 
of the characteristics that we have just discussed for the tree-level. 

3.2.2 Electroweak Yukawa-type contributions, novel charac- 
teristics 



t H 9 b 




Figure 3.2: Sample of one-loop diagrams related to the Yukawa interaction in the SM. 
Xw represents the charged Goldstone boson. 

Indeed, look at the two contributions arising from the one loop electroweak correc- 
tions given in Fig. 13.21 Now the Higgs can attach to the top or to the W . Therefore 
these contributions do not vanish in the limit XbbH = 0. The mass insertion in what 
we called F"*^"^ is proportional to the top mass and is not negligible. In fact the 
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diagrams in Fig. 13.21 show the charged Goldstone boson in the loop. The latter trig- 
gers a. t —>■ bxw transition whose dominant coupling is proportional to the Yukawa 
coupling of the top. We will in fact be working in the approximation of keeping 
only the Yukawa couplings. This reduces the number of diagrams and if working in 
the Feynman gauge as we do in this computation, only the Goldstone contributions 
survive. The neutral Goldstone bosons can only contribute corrections of order 
(see section for the Feynman rules). We will neglect these 0{Xl) contributions 
at the amplitude level. However the order 0{Xb) corrections will be kept. All the 
corrections are then triggered hj t ^ bxw and apart from the QCD vertices, only 
the Yukawa vertices shown in Fig. 13.31 below are needed to build up the full set of 
electroweak corrections. Note that in the MSSM, the Higgs coupling to the fermion 




Figure 3.3: Relevant vertices appearing at one loop, eu = A^/At and A is the Higgs 
self- coupling, related to the Higgs mass in the Standard Model. The relations to the 
gauge couplings can be obtained by comparing to the SM Feynman rules given in 
sectioning 



/, XffH, can involve other parameters beside the corresponding Yukawa coupling A/, 
as shown in subsection 11.5.21 The Higgs coupling to the charged Goldstone involves 
the Higgs self-coupling or Yukawa coupling of the Higgs, A = Mjj/2v'^ proportional to 
the square of the Higgs mass. The latter can be large for large Higgs masses. These 
considerations allow to classify the contributions into three gauge invariant classes. 
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Figure 3.4: All the diagrams in each group can be obtained by inserting the two gluon 
lines or one triple gluon vertex (not shown) to all possible positions in the generic 
bottom line, which is the first diagram on the left. We have checked the number of 
diagrams through Grace- loop fj^ . 



3.2.3 Three classes of diagrams and the chiral structure at 
one-loop 

All the one-loop diagrams are classified into three gauge invariant groups as displayed 
in Fig. 13. 4[ The Higgs couples to the bottom quark in the first group (Fig. 13.4b .). to 
the top quark in the second group (Fig. 13.4b ) and to the charged Goldstone boson 
in the third group (Fig. 13.4b ). As shown in Fig. 13.41 each class can be efficiently 
reconstructed from the one-loop vertex bbH, depending on which leg one attaches the 
Higgs, by then grafting the gluons in all possible ways. We have also checked explicitly 
that each class with its counterterms, see below, constitutes a QCD gauge invariant 
subset. Note that these three contributions depend on different combinations of 
independent couplings and therefore constitute independent sets. 
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The chiral structure t — > bxw impacts directly on the structure of the hehcity 
amphtudes at one-loop. The split of each contribution according to T'^'"'^"' and r°'^'^, 
see Eq. (13.21) will turn out to be useful and will indicate which helicity amplitude can 
be enhanced by which Yukawa coupling at one-loop. We show only one example in 
class (6) of Fig. 13. 4[ It is straight forward to carry the same analysis for all other 
diagrams. We choose the first diagram in group (b) in Fig. 13.41 For clarity we will 
here take = 0, we have already shown how rrify insertions are taken into account, 
see Eq. (13.21) . Leaving aside the colour part which can always be factorised out (see 
Appendix [B]) and the strong coupling constant, we write explicitly the contribution 
of this diagram as 

Ai(Ai, As; A3, A4) = AitHAiM(A3,p3)/(Ai,pi)^Cfei%^(A2,P2)M(A4,P4)- (3.3) 

Pl3 P24 

Cfji is the Yukawa vertex correction. In D-dimension, with q the integration variable, 
the momenta as defined in Fig. 13.11 with pij = Pi + pj and pij = pj — pi we have 

d^q {Pr - ebtPL){mt + <i+pi3){mt + 4 - MjPL - SuPr) , . 
{2nri ^M'^ - q2)[ml - [q + p,^f][ml - {q - p2,y] ' 

where Su = Afe/A* as defined in Fig. 13.31 The numerator of the integrand of Eq. (13. 4p . 
neglecting terms of (9(A^), can be re-arranged such as 

Ai(Ai,A2; A3,A4) - su {ml + + - ^2a)) 

" V ' 

Y'cven 

+ mtP/j (2^ + ^13 - ^24) • (3.5) 
^ V ' 

This shows explicitly that F"*^*^ structures with a specific chirality, Pr, can indeed be 
generated. They do not vanish as Xhbh 0. The even one- loop structures on the 
other hand are 0{\)- The structure in class (c), Higgs radiation off the charged 
Goldstones, is the same. For class (a), radiation off the 6-quark, the structure of the 
correction is different, the odd part is suppressed and receives an 0{\b) correction. 
To summarise, with mb = 0, making explicit the Yukawa couplings and the chiral 
structure if any, for example Pr, that characterise each class and comparing to the 
leading order, one has 
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Despite the existence of the simple relation Xj/h = —Xf/V2 in the SM, we have 
kept XffH and A/ separate to distinguish their different origins (A/ comes from the 
Goldstone couplings). As discussed in subsection ll.5.2[ in the MSSM X^bH is enhanced 
by tan j3 but not Xb- We clearly see that all one-loop r^''^" contributions vanish in the 
limit Xb = and XbbH = 0. On the other hand this is not the case for the one-loop 
■podd contribution belonging to class (b) and (c). However for these contributions to 
interfere with the tree- level LO contribution requires a chirality flip through a m?, 
insertion. Therefore in the SM for example, the NLO cross section is necessarily of 
order ml, like the LO, with corrections proportional to the top Yukawa coupling for 
example. On the other hand, in the limit of XbbH = 0, the tree level vanishes but 
gg — » bbH still goes with an amplitude of order g'^XfXuH or g^X^X^^n- For XbbH 7^ 
these contributions should be considered as part of the NNLO "corrections" however 
they do not vanish in the limit m;, — > (or XbbH = 0) while the tree level does. These 
contributions can be important and we will therefore study their effects. For these 
contributions at the "NNLO" we can set rrib = 0. 

The classification in terms of structures as we have done makes clear also that 
the novel one-loop induced F"'^'^ contributions must be ultraviolet finite. This is 
not necessarily the case of the T'^^^"' structures where counterterms to the tree-level 
structures are needed through renormalisation to which we now turn. 
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3.3 Renormalisation 



We use an on-shell (OS) renormalisation scheme exactly along the lines described in 
subsection 11.2.61 Ultraviolet divergences are regularised through dimensional regular- 
isation. In our approximation we only need to renormalise the vertices bbg and bbH 
as well as the bottom mass, m^. For the bbg vertex, its counterterm reads 

5^,, = 2g,Y{SZli'PL + 5Zl';^Pn). (3.6) 

1 /2 

5Zj^ ^ are calculated by using Eq. fll.35p . For this one needs to know the coefficients 
-^1,7,57 which are very simple in our approximation: 




We get 



K,{q') = -K,,{q') = ^,(Cuv-2F,{t,W)) 

with Cuv = - - + In 47r , = 4 - 2e, 
e 

Fn{A,B) = J (ixx"ln ^(1 - a;)m^ + xml - a;(l - a;)g^ 



(3.7) 



The reason we get K-y{q^) = — i^5^(g^) is due to the particular chiral structure 
of the t bxw loop insertion. In particular for rrib = 0, one recovers that these 

1/2 1 /2 

corrections only contribute to and not 5^^^ • 

For the bbH vertex with the bare Lagrangian term CbbH = one needs 

to renormalise nib, v, '^h- With v — > v{l + 6v) and the rules given in 

Eq. f ll.26p we get the counterterm 

6nih 



^bbH — \bH 



nib 



I'b , cryl/2 , cryl/2 , /j;ryl/2 r X 



{3.1 
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Again SZl^^^^ and Srrib are calculated by using Eqs. (11.351) and (13.71) . SZ]^'^ is calcu- 
lated through Eq. (11.311) where 11^ (g^) comes from the diagrams with heavy particles 
in the loop. As we are interested in the Yukawa corrections, these particles are the 
top quark, the Higgs boson, the W-Goldstone bosons and the Z-Goldstone boson. 
Indeed, those corrections include all the leading Higgs couplings: Xuh and Xhhh- We 
get 



^Cuv - \f,{H, H) - F,{W, W) - ^Fo{Z, Z) - 3(1 - InM^) 



+ 



3m 



* {q'iCuv - Fo{t,t)] + 4[1 - Inm? + Fo{t,t)]} 



(3.9) 



Prom this and Eq. (I1.3ip we get 



SZ 



1/2 
H 



Re 



3A? 



Cuv - Fo{t, t) - M^Goit, t) + 4m2Go(t, t) 



9Go{H, H) + 2Go(TV, W) + G'o(^, Z 



Gn{A,B) 



q^—Fn{A,B) = q^ I dx- 



dq^ 



x"x(l — x) 



1 — x)m\ + xm\ — x{l — x)q- 



:.(3.10) 



For 5v, we use the relation 



V = , Sw 



Ml 



(3.11) 



to write 5v in terms of (5M^, (5M| 



5v 



2M| 



(3.12) 



We do not need to renormalise e since we are interested only in the Yukawa sector 
hence do not touch the photon. We now use Eq. (11.381) to calculate and (5M|. 

Like in the case of 11^, 11^ and Ilf include all the leading contributions related to 
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the Goldstone bosons: 



U^{M^ 



WJ 



n|(M| 



+ 



ML 



ML 

^^2^2 



[F^iH, W)-Fo{H,W)] 



+ 



InM^ 



[3Cc/y-6Fi(6,t)] 



+ 



mr 



^^2^2 



[Fr{H,Z)-F,{H,Z)] 



[3Cf/v-3Fo(t,t)] 



InM^ 



(3.13) 



Then we get 

5v 



Re 



3A? 



Cc/y-2Fi(6,t) 



- A 



Fo{H, W)-F,{H, W)--lnMl 



3A? 



r2 



'VI/ 



- + A( - Fo{H, W) + ly)) 

Q \2 

^Fo(t, t) + A(Fo(i/, Z) - F,{H, Z)) 



(3.14) 



From Eqs. f l3.10p and f l3.14p . one clearly sees that {SZjj ~Sv) is UV-finite. Therefore, 
by looking at Eqs. (13.61) and (13. Sp . we conclude that to make all the contributions 



of diagrams in Fig. 13.41 UV-finite, it is sufficient to renormalise the mass and wave 

1 /2 

function of bottom-quark as done above. On the other hand {SZjj — 6v) can be 
seen as a universal correction to Higgs production processes. We will include this 
correction as it has potentially large contributions scaling like A^ and A which fall 
into the category of the corrections we are seeking. 

In the actual calculation, the counter term 5^^^ belongs to class (a) in the clas- 
sification of Fig 13.41 This makes class (a) finite. The counterterm we associate to 
class (6) is the part of 6bbH from the t — * bxw loops and therefore does not include 
what we termed the universal Higgs correction, i. e does not include the contribution 
{SZ^^'^ — 6v) . This is sufficient to make class (6) finite. In our approach (c) is finite 

-1 /o 

without the addition of a counterterm. We will keep the {SZ^I — 6v) contribution 
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separate from the contributions in classes (a), (6), (c). We will of course include it in 
the final result. 

3.4 Calculation details 

We have written two independent codes. In the first one we set = in all 
propagators and other spinors that emerge from the helicity formalism we follow. In 
this limit, the helicity formalism is very much simplified and the expression quite 
compact. This code is in fact subdivided in two separate sub-codes. One sub-code is 
generated for the "even" part (constituted by the r*''^^" contributions, see Eq. fl3.2l) ) 
and the other by the "odd" part. We also generate a completely independent code 
for the case ^ where in particular we use the helicity formalism with massive 
fermions. Details of the helicity formalism that we use are given in Appendix El 

The steps that go into writing these codes are the following. In the first stage, 
we use FORM|90] to generate expressions for the tree level and one loop helicity 
amplitudes. Each helicity amplitude is written in terms of Lorentz invariants, scalar 
spinor functions {A,B,C)\^Xj defined in Appendix IX] and the Passarino-Veltman[32j 
tensor functions for a tensor of rank M for A^-point function. We have also 
sought to write the contribution of each amplitude as a product of different structures 
or blocks that reappear for different graphs and contributions. For example colour 
factorisation is implemented, this further allows to rearrange the amplitude into an 
Abelian part and a non-Abelian part which will not interfere with each other at 
the matrix element squared level. The helicity information is contained in a set of 
basic blocks for further optimisation. Another set of blocks pertains to the loop 
integrals and other elements. The factorisation of the full amplitude in terms of 
independent building blocks is easily processed within FORM. These building blocks 
can still consist of long algebraic expressions which can be efficiently abbreviated 
into compact variables with the help of a Perl script which also allows to convert 
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the output of FORM into the Fortran code ready for a numerical evaluation. More 
details on the FORM code as well as the optimisation we implemented can be found 
in Appendix [Bi 

3.4.1 Loop integrals, Gram determinants and phase space 
integrals 

The highest rank M of the Passarino-Veltman tensor functions with M < N that 
we encounter in our calculation is M = 4 and is associated to a pentagon graph, 
= 5. We use the library Loop Tools [HSl ES] to calculate all the tensorial one loop 
integrals as well as the scalar integrals, this means that we leave it completely to 
LoopTools to perform the reduction of the tensor integrals to the basis of the scalar 
integrals. In order to obtain the cross section one needs to perform the phase-space 
integration and convolution over the gluon distribution function (GDF), g{x, Q) with 
Q representing the factorisation scale. We have 



The integration over the three body phase space and momentum fractions of the 
two initial gluons is done by using two "integrators": BASES [37| and DADMUL|91j. 
BASES is a Monte Carlo that uses the importance sampling technique while DAD- 
MUL is based on the adaptive quadrature algorithm. The use of two different phase 
space integration routines helps control the accuracy of the results and helps detect 
possible instabilities. In fact some numerical instabilities in the phase space inte- 
gration do occur when we use DADMUL but not when we use BASES which gives 




\A{gg bbH)\^6^{pi + P2 - - P4 - Ps) , 

(3.15) 
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very stable results with small integration error, typically 0.08% for 10^ Monte Carlo 
points per iteration (see section IC. 21 for more details). For the range of Higgs masses 
we are studying in this chapter, the instabilities that are detected with DADMUL 
were identified as spurious singularities having to do with vanishing Gram determi- 
nants for the three and four point tensorial functions calculated in LoopTools by using 
the Passarino-Veltman reduction methocj^. Because this problem always happens at 
the boundary of phase space, we can avoid it by imposing appropriate kinematic cuts 
in the final state. In our calculation, almost all zero Gram determinants disappear 
when we apply the cuts on the transverse momenta of the bottom quarks relevant for 
our situation, see section [3.5.1l for the choice of cuts. The remaining zero Gram deter- 
minants occur when the two bottom quarks or one bottom quark and the Higgs are 
produced in the same direction. Our solution, once identified as spurious, was to dis- 
card these points by imposing some tiny cuts on the polar, 6, and relative azimuthal 
angles, of the outgoing 6-quarks, the value of the cuts is Ol'^^ = \ sin0^|c„i = 10^^. 
DADMUL then produces the same result as BASES within the integration error. 



3.4.2 Checks on the results 

i) Ultraviolet finiteness: 

The final results must be ultraviolet (UV) finite. It means that they should be 
independent of the parameter Cuv defined in Eq. (13.71) . In our code this parameter 
is treated as a variable. The cancellation of Cuv has been carefully checked in our 
code. Upon varying the value of the parameter Cuv from Cuv = to Cuv = 10^, the 
results is stable within more than 9 digits using double precision. This check makes 
sure that the divergent part of the calculation is correct. The correctness of the finite 
part is also well checked in our code by confirming that each helicity configuration is 
QCD gauge invariant. 

■^The reduction of the five point function using the method of Denner and Dittmaier [12 [31] 
which avoids the Gram determinant at this stage as implemented in LoopTools gives very stable 
results. 
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ii) QCD gauge invariance: 

In the physical gauge we use, the QCD gauge invariance reflects the fact that the 
gluon is massless and has only two transverse polarisation components. In the helicity 
formalism that we use, the polarisation vector of the gluon of momentum p and helicity 
A is constructed with the help of a reference vector q, see Appendix |X] for details. 
The polarisation vector is then labelled as e^(p, A;g). A change of reference vector 
from q to q' amounts essentially to a gauge transformation (up to a phase) 



QCD gauge invariance in our case amounts to independence of the cross section in 
the choice of the reference vector, q. We have carefully checked that the numeri- 
cal result for the norm of each helicity amplitude at various points in phase space 
is independent of the reference vectors say qi^2 for gluon 1 and 2, up to 12 digits 
using double precision. By default, our numerical evaluation is based on the use of 
9i,2 = {P2,Pi)- For the checks in the case of massive b quarks the result with the 
default choice gi^2 = {P2,Pi) is compared with a random choice of qi^2, keeping away 
from vectors with excessively too small or too large components, see Appendix El for 
more details. 

iii) As stated earlier, the result based on the use of the massive quark helicity ampli- 
tude are checked against those with the independent code using the massless helicity 
amplitude by setting the mass of the b quark to zero. This is though just a consistency 
check. 

iv) At the level of integration over phase space and density functions we have used 
two integration routines and made sure that we obtain the same result once we have 
properly dealt with the spurious Gram determinant as we explained in section 13.4.11 

v) Moreover, our tree level results have been successfully checked against the results 




(3.16) 



of CalcHEP[92]. 
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3.5 Results: Mh < 2Mw 

3.5.1 Input parameters and kinematical cuts 

Our input parameters are a(0) = 1/137.03599911, Mw = 80.3766GeV, Mz = 
91.1876GeV, a,(Mz) = 0.118, = 4.62GeV, = 174.0GeV with cw = Mw/Mz- 
The CKM parameter Vtb is set to be 1. We consider the case at the LHC where 
the center of mass energy of the two initial protons is ^/s = 14TeV. Neglecting the 
small light quark initiated contribution, we use CTEQ6Lp3l [9ll [95l [96] for the gluon 
distribution function (GDF) in the proton. The factorisation scale for the GDF and 
energy scale for the strong coupling constant are chosen to be Q = Mz for simplicity. 

As has been done in previous analyses [STl [97] , for the exclusive bbH final state, 
we require the outgoing b and b to have high transverse momenta |p^^| > 20GeV 
and pseudo-rapidity \ri^'^\ < 2.5. These kinematical cuts reduce the total rate of the 
signal but also greatly reduce the QCD background. As pointed in [50] these cuts 
also stabilise the scale dependence of the QCD NLO corrections compared to the case 
where no cut is applied. In the following, these kinematical cuts are always applied 
unless otherwise stated. 

Talking of the NLO QCD scale uncertainty and before presenting our results, 
let us remind the reader of the size of the QCD corrections. Taking a renormali- 
sation/factorisation scale as we take here at Mz, the QCD corrections in a scheme 
where the bottom Yukawa coupling is taken on-shell amount to ~ —22% for a Higgs 
mass of 120GeV. 

3.5.2 NLO EW correction with XbbH 7^ 

The cross sections with two high-p^ bottom quarks at LO and NLO at the LHC 
are displayed in Fig. 13.51 as a function of the Higgs mass. The NLO EW correction 
reduces the cross section by about 4% to 5% as the Higgs mass is varied from llOGeV 
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Figure 3.5: Left: the LO and NLO cross sections as functions of Mh- Right: the 
relative NLO EW correction normalized to tree level ctlo- {o); {b), (c) correspond 
to the three classes of diagrams as displayed in Fig. \3.4\ to which counterterms are 



1 /2 

added (see section [Ql) . {5Z^ — 6v) is the correction due to the universal correc- 
tion contained in the renormalisation of the bbH vertex. "Total" refers to the total 
electroweak correction, of Yukawa type, at one-loop. 



to 150GeV. The first conclusion to draw is that this correction is small if we compare 
it to the QCD correction or even to the QCD scale uncertainty. Considering that we 
have pointed to the fact that the contributions could be grouped into three gauge 
invariant classes that reflect the strengths of the Higgs coupling to the b, the t or its 
self-coupling, one can ask whether this is the result of some cancellation. It turns 
out not to be the case. All contributions are below 3%, see Fig. 13.51 Class (a) with 
a Higgs radiated from the bottom line is totally negligible ranging from —0.09% to 
—0.06%. We have failed in finding a good reason for the smallness of this contribution 
compared to the others. Those due to the Higgs self-coupling are below 1%. Radiation 
from the top contributes about —2% and is of the same order as the contribution of 
the universal correction. We had argued that the Yukawa corrections brought about 
by the top might be large. It seems that the mass of the top introduces also a large 
scale which can not be neglected compared to the effective energy of the hard process 
even for LHC energies. 
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Figure 3.6: Effect of the NLO electroweak corrections on the pseudo-rapidity and 
transverse momentum distributions of the Higgs for Mh = 120, 150 GeK The relative 
corrections daNLo/dcTLO — ^ is also shown. 



The NLO corrections are spread rather uniformly on all the distributions we have 
looked at. We have chosen to show in Fig. 13.61 the effect on pseudo-rapidity and 
transverse momentum distributions of the Higgs for two cases = 120GeV and 
Mh = 150GeV. As Fig. 13.61 shows the relative change in these two distributions is 
sensibly constant especially for Mh = 120GeV. For Mh = 150GeV, the corrections 
are largest for p^ around 140GeV, however this is where the cross section is very 
small. A similar pattern, i.e. a constant change in the distributions, is observed for 
the bottom variables. 
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3.5.3 EW correction in the limit of vanishing A^^ij 




MHlGeV] MufGeV] 



Figure 3.7: The one-loop induced cross section as a function of Mh in the limit 
of vanishing bottom-Higgs Yukawa coupling. The right panel shows the percentage 
contribution of this contribution relative to the tree level cross section calculated with 

^bbH 7^ 0. 

The cross section for X^bH = can be induced at one-loop through the top loop. 
This "NNLO" contribution rises rather quickly as the Higgs mass increases even in 
the narrow range Mh = 110 — ISOGeV as can be seen in Fig. 15. Ill Indeed relative 
to the tree level, the cross section with Mh = 120GeV amounts to 3% while for 
Mh = ISOGeV it has increased to as much as 17%. Going past Mh > 2M\y we 
encounter a Landau singularity [30] (a pinch singularity in the loop integral) from 
diagrams like the one depicted in Fig. 13.21 (right) with the Higgs being attached to the 
W^s or their Goldstone counterpart. It corresponds to a situation where all particles 
in the loop are resonating and can be interpreted as the production and decay of the 
tops into (longitudinal) W^s with the later fusing to produce the Higgs. This leading 
Landau singularity is not integrable, at the level of the loop amplitude squared and 
must be regulated by the introduction of a width for the unstable particles. This 
issue together with a general discussion of Landau singularities will be considered in 
the next chapters. Fig. 13.81 shows the pseudo-rapidity and transverse momentum 



3.5. Results: Mh < 2Mw 



61 




-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



> 


35 


CD 




eg 






30 


n. 




H H 
Q. 
T3 


25 


"B 










20 




15 




10 




5 








xlQ- 



pp->bbH 
Vs=14TeV 



\ M„=150GeV 



MH=120GeV 




20 40 60 80 100 120 140 160 180 200 

pSGeV] 



> 

C5 0.05 



0.04 



0.03 



0.02 



0.01 



, Mu=150GeV 




20 40 60 80 100 120 140 160 180 200 

p^IGeV] 



-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 




20 40 60 80 100 120 140 160 180 200 

p'JIGeV] 



60 
50 7 
40 
30 

20 b 
10 




n MH=150GeV 



pp- 


^bbH 




14TeV 



MH=120GeV 




20 40 60 80 100 120 140 160 180 200 

p^[GeV] 



Figure 3.8: The pseudo-rapidity of the Higgs and transverse momentum distributions 
of the Higgs and the bottom for Mh = 120, ISOGeF arising from the purely one- 
loop contribution in the limit of vanishing LO (XtbH = Oj. Its relative percentage 
contribution da{XbbH = 0)/d<^LO is also shown 
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distributions of the Higgs as well as the the pt of the bottom for two cases Mh = 
120GeV and Mh = 150GeV in the limit of vanishing bottom-Higgs Yukawa coupling. 
These distributions are significantly different from the ones we observed at tree-level 
(and with the electroweak NLO corrections), see Fig. 13.61 The Higgs prefers being 
produced at high value of transverse momentum, about 130GeV. In the case of a 
Higgs with Mh = ISOGeV this contribution can significantly distort the shape of 
the pff distribution for hight pf^ with a "correction" of more than 70% over a rather 
large range. The distribution in the pt of the bottom is also very telling. The new 
contributions do not produce the bottom preferentially with low p^ as the case of the 
LO contribution. 



3.6 Summary 

We have calculated the EW radiative corrections triggered by the Yukawa coupling 
of the top to the process pp bbH at the LHC through gluon fusion in the Standard 
Model. This process is triggered through Higgs radiation of the bottom quark with a 
small coupling proportional to the mass of the bottom. Yet in order to analyse this 
coupling, precision calculations that include both the QCD and electroweak correc- 
tions are needed. In this perspective, to identify the process one needs to tag both 
6-jets. Our calculation is therefore conducted in this kinematical configuration. 

Inserting a top quark loop with a Yukawa transition of the type t bxw, Xw 
is the charged Goldstone, allows now the Higgs to be radiated from the top or from 
the Goldstone boson. The latter coupling represents the Higgs self-coupling and 
increases with the Higgs mass. The former, the top Yukawa coupling, is also large. 
As a consequence, the one-loop amplitude gg — > bbH no longer vanishes as the Higgs 
coupling to 6's does, like what occurs at leading order. 

We find that in the limit of vanishing XbbH, the one-loop induced electroweak pro- 
cess should be taken into account for Higgs masses larger than 140GeV or so. Indeed, 
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though this contribution is quite modest for a Higgs mass of llOGeV it increases quite 
rapidly as the Higgs mass increases, reaching about 17% of the leading order value, 
calculated with mf, = 4.62GeV, for Mh = 150GeV. For these new corrections to in- 
terfere with the leading order requires helicity flip. Therefore at next-to-leading order 
in the Yukawa electroweak corrections, all corrections involve either a bottom mass 
insertion or a bottom Yukawa coupling. At the end the total Yukawa electroweak 
NLO contribution brings in a correction which is within the range —4% to —5% for 
Higgs masses in the range llOGeV < Mh < ISOGeV. They are therefore negligible 
compared to the NLO QCD correction and even the remaining QCD scale uncer- 
tainty. This modest effect translates also as an uniform rescaling of the distributions 
in the most interesting kinematical variables we have looked at (pseudo-rapidities 
and pt of both 6-quarks and the Higgs). This is not the case of the one-loop induced 
contributions which survive in the limit of (and XbbH —* 0). Here the distri- 

butions for the Higgs masses where the corrections for the total cross section is large 
are drastically different from the LO distributions. A summary for the corrections 
including the NLO with X^hu ^ and the part of the NNLO counted as loop induced 
in the hmit XbbH — > is shown in Fig. 13.91 
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Figure 3.9: 6ew = ^nlo + 
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as a function of Mh ■ 
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The analysis we have performed in this chapter does not cover Higgs masses over 
150GeV and rests within the range of Higgs masses preferred by indirect precision 
measurements. In fact as the threshold for H —>■ WW opens up, important phenom- 
ena take place. Foremost Landau singularities, or pinch singularities in some loop 
integrals, develop. This important issue will be addressed in the next two chapters. 

There is another contribution which does not vanish for vanishing XbbH and which 
contributes to gg — > bbH through a closed top quark loop. This contribution rep- 
resents gg Hg* Hbb. We have not included this contribution in the present 
work as we do not consider it to be a genuine bbH final state. This correction can be 
counted as belonging to the inclusive gg ^ H process. The same line of reasoning 
has been argued in [83]. Nonetheless from the experimental point of view it would be 
interesting to include all these effects together with the NLO QCD corrections and 
the electroweak corrections that we have studied here. 
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Landau singularities 



Let us start this chapter by quoting a paragraph in p. 30 of the important book The 
Analytic S-matrix written by Eden, Landshoff, Ohve and Polkinghorne (1966) |29j : 

For an individual Feynman integral, the singularities, and hence the an- 
alytic structure of the integral, arise from singularities in the integrand. 
Those that arise from ultraviolet divergences are removed by renormali- 
sation and will not concern us. The singularities that do concern us come 
from zeros of the denominator factors. These denominator factors are the 
same for scalar particles as they are for particles with spin, namely of the 
general form 

p^-m^ + ze (4.1) 
for mass and four-momentum p. 

The issue of loop-integral singularities is also discussed in other text books [981 ESI 
llUUj . see also the recent lecture notes of Denner for a practical review [lUlj . 

In this chapter we will restrict ourselves to one-loop Feynman integrals. To set 
the mathematical background, we first consider the case of some general complex 
integrals. 
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4.1 Singularities of complex integrals 

We would like to consider two simple examples taken from |29j . 

(i) With w complex 

f{w)= / = \n{h-w)-\n{a-w). (4.2) 

Ja z-w 

(ii) With a complex and real x 



Q {x — w){x — a) 
[ln(l - w) - ln(l - a)] - [In(-w) - In(-a)] 



w — a 



(4.3) 



In the first example we see that the integrand l/{z — w) is singular at the end-points 
a, b if w = a, b, corresponding to the singularity of the logarithm at these points. We 
see clearly that fiw) is always singular if w = a,b whatever the integration contour 
is. w = a,b are called end-point singularities. We notice also that if w lies somewhere 
on the integration contour but not at the end points then we can always use Cauchy's 
integral theorem to deform the contour to avoid the pole. 

The second example is more interesting. Again we see that w = 0,1 are end-point 
singularities. One notices also that there may be a problem when w = a. In order to 
see the point, we consider the case where a = a — ie with e is infinitesimal positive 
and w = w — ip with p is also infinitesimal but can be positive or negative, a and w 
now should be considered as real numbers. We distinguish the three following cases: 

1. For < a < 1 and w a we have ln(— a + ie) = ln(a) + iir, ln(— w + ip) = 
ln{w) + iTT sign(p) and 

f(w—*a) = i7r[sign(p) — 1] lim . (4.4) 

1 — a a w — a 
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2. For a < then all logarithmic arguments in Eq. 04.31] are strictly positive. We 
then can forget about e and p to get 

f{w a) = --^ - -. (4.5) 
I — a a 

3. For a > 1 then all logarithmic arguments in Eq. fl4.3p are negative. We have 

f{w — > a) = — h {i7r[sign(p) — 1] — 'i7r[sign(p) — 1]} lim 



1 — a a w — a 

7^ - -■ (4.6) 

1 — a a 




Figure 4.1: (a) Two singularities in the complex plane coming together and pinch- 
ing the contour of integration, [0,1] along the real axis, (b) and (c) Examples of 
coincidence of two singularities that do not pinch the integration contour. 



We remark immediately that w = a is a singular point if and only if two conditions 
are satisfied: 

< a < 1 and sign(p) = -1, (4.7) 

which simply mean that the real integration contour [0, 1] is trapped when w ^ a 
(see Fig. 14.1b .). This is a pinch singularity. We learn also from the first case that if 
sign(p) = 1, i.e. both w and a approach the real axis from the same side, then the 
contour is not trapped and there is no singularity (see Fig. 14.1b ). In the second and 
third cases there is no singularity because w and a come together nowhere near the 
contour (see Fig. 14.1b ). In this simple example, it is easy to find out the situation 
where the integration contour is trapped. However, in a general case with more 
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than one integration variables finding whether the contour is actually trapped is 
very difficult. The "e-prescription" and multivalued function play crucial roles in 
the second example. These two elements will appear again when one considers loop 
integrals E. 

Those two examples will guide us all the way. For now they lead us to two 
important statements for simple complex integrals. For this we ffist define the context 
Let g(w,z) be an analytic function of two complex variables and let C be some 
finite contour in the complex z-plane. Define a function f{w) by 

/H= / g{w,z)dz. (4.8) 
Jc 

It is supposed that the singularities of the integrand g are known and that their 
locations in the 2;-plane are 

z = Zr{w), r = l,2,.... (4.9) 

f{w) is analytic as long as C does not intersect with any Zr{w) or can be deformed 
accordingly. From the above examples, we see that f{w) can be singular at Wi for 
one of two reasons: 

(i) End-point singularities: One of the singularities Zr reaches one of the end-points 

of the contour C when w ^ wi. No deformation of C can avoid them and wi 
is a singularity of f{w). 

(ii) pinch singularities: If two (or more) singularities approach the contour from 

opposite sides and coincide, the contour C will be trapped between them and 
no deformation can avoid them, wi is a singularity of fiw). 



^In fact, we will see later in this thesis that one way to deal with pinch singularities of loop 
integrals, named as Landau singularities, is to kill the e-prescription by introducing the widths of 
internal unstable particles. The width moves the singularities away from the real axis, so they do 
not occur in the physical region (the real axis). 
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In order to deal with loop integrals, we have to generalise the above consideration 
to the case of multiple integrals. Consider 

/H= / l[dz,g{w,z,), (4.10) 

where if is a hyper contour in the multi-dimensional complex Zj-space. The singular- 
ities of the integrand g{w,Zi) are given by various equations 

Sriw,Zi) = 0, r = l,2,... . (4.11) 

The boundary of H is specified by 

Ss{w,Zi) = 0, s = l,2,... . (4.12) 

Singularities occur when a surface of singularity meets a boundary surface or when 
the hypercontour H is pinched between two or more surfaces of singularities. More 
precisely, it may be shown that a necessary condition of singularity is that there exists 
a set of complex parameters a^, not all equal to zero such that [29l [98] 

arSr = 0, for each r, 

(so that either or Sr = 0), 

asSs = 0, for each s, 
and, for each integration variable Zi, 

= 0. (4.13) 

The last condition expresses that the hypersurfaces are tangent. This is only a nec- 
essary condition and does not guarantee that the hypercontour H is pinched. This is 
easy to understand if we look back at the example two. If the hypersurfaces come to 
be tangent from one side of the hypercontour or they are tangent at nowhere on the 
hypercontour then they are harmless. 

To find necessary and sufficient conditions for a pinch singularity of multiple in- 
tegrals is very difficult and requires homology theory. This is obviously beyond the 
scope of this thesis and we refer to [1021 1103j for more study. 
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4.2 Landau equations for one- loop integrals 




We would like to clarify the important terminology used in this and the next chap- 
ters. The terminology "phase space" means the kinematical allowed region obtained 
by using the energy-momentum conservation on external momenta. The terminol- 
ogy "physical region" means the physical integration contour of the one-loop integral 
defined in Eq. (14.201) . It means {xi = x*,Xi > 0,qi = q*} where are Feynman pa- 
rameters and Qi are loop momenta. The real condition on Xi and qi accords with the 
ie prescription. 

We are now in the position to apply the previous analysis for one-loop Feynman 
integrals. Consider the one-loop process -Fi(pi) + -^2(^2) + • • • -^iv(pAr) — > 0, where Fj 
stands for either a scalar, fermion or vector field with momentum pi as in the figure 
opposite. The internal momentum for each propagator is qi with i = 1, . . . N. The 
scalar one-loop integral reads 
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In momentum space 

rpN 

° - / (27r)^.nf=iA' 

i 

Di = Qi - mf + ie, Qi = q + Vi, = ^Pj, (4.14) 

where mj are the masses of internal particles, all the external particle are out- 
going, Pi = gj+i — qi are the external momenta. 



Introducing Feynman parameters and integrating over q one gets 

_ (-i)^r(iv-D/2) ^(ELx.-i) 

° ~ (4^0^72 dXi---dXN y^-^^> 

with 

1 ^ 

A = - ^ XiXjQij - ie, Qij = ml + m^j - - rj f. (4.16) 

In representation (14.141) the physical hypercontour is [— oo,oo]'^ along the real axes, 
by definition. Each factor in the denominator is protected from zero by ie with e > 0. 
Thus the integral can be singular only when e — > 0. However, even in the case when 
e ^ and the singularities of the integrand reach the real hypercontour one can 
deform the real hypercontour to be complex to avoid the singularities. Indeed, the 
generalised definition of Feynman integral is Eq. (14.141) and its analytical continu- 
ation As already pointed out in the second example of the previous subsection, 
a singularity arises when the real hypercontour is pinched and therefore cannot be 
deformed. The pinch singularity of a Feynman integral is called Landau singularity. 
Necessary condition for this singularity is that equations (14.131) have solution with 
real q. The boundaries are at infinity so we can set all = 0. With Si = qf — mf 



^Based on this definition, a method called contour deformation technique to calculate Feynman 
integrals numerically has been realised in practice [1041 1105j . This method works as long as the 
contour can be deformed to avoid singularities, i.e. when the singularities approach the contour 
from the same side like in Fig. 14.1b . 
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equations fl4.13l) become 

These are the Landau equations corresponding to representation (I4.14p . If equations 
(14.1 7p have a solution with some complex Oi ^ and real qi then one can say that 
integral (14.141) may have a Landau singularity. 

The drawback of this representation is that there are not any constraints on the 
range of a^. By introducing Feynman parameters and requiring them to be real (this 
is the purpose of ie prescription), one imposes more constraint on the physical region. 
It will be clear later on that for Landau singularities to be in that physical region, at 
must be real and positive. 

For representation (14.151) . the physical hypercontour is [0, oo]^ along the real axes 
by definition [fj. Thus the boundaries are Si = Xi = for all i. The hypersurface of 
singularities of the integrand is 5 = A = 0. One might think that the hypercontour 
cannot be pinched now because there is only one hypersurface of singularities and 
the singularities are just end-point singularities if they occur. Indeed, a hypersurface 
can be very complicated and contains several sub-hypersurfaces which can trap the 
hypercontour. This is very difficult to imagine but the condition for that to happen 
is very simple, just like the extreme condition {29]: 

Other way to get the condition is to use equations (14.131) . We can set a = 1 since 
there is only one singular hypersurface. One gets 

^ = 0, , , 

4.19 

\fi = 0. 

' axi 

The first condition comes from the second equation of (14.131) . The second condition 
comes from the first and third equations of (14.131) . If one picks up the solution Xj = 



■^One might think that the hypercontour should be [0, 1]^. However, it can be expanded to 
infinity because of the Dirac delta function in the integrand. 
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from the second condition of fl4.19l) then that is for end-point singularities. Otherwise, 
one has the condition for pinch singularities, which is the same as Eq. fl4.18p . So 
everything is consistent. Since A is a homogeneous function of Xj, the first equation 
in (14.191) is automatically satisfied when the second is. Necessary condition for integral 
f l4.15p to have Landau singularities is that equations fl4.19p have solution with some 
Xi > 0, Xi are real. The drawback of this representation is that equations fl4.19p do 
not tell us anything about g^. 

There exists a representation which contains all the advantages of Eqs. f l4.17p and 
f l4.19p . That is the mixed representation of Feynman integrals in the space of real 
momentum and real Feynman parameters B 

= T{N) / dx, ■ ■■dx^SiY.x. - 1) / J ^ (4.20) 

The physical hypercontours are real [—00,00]^ for q and real [0, 00]^ for Xj. The 
boundaries are Si = Xi = for all i. The hypersurface of singularities of the integrand 
is S* = J2iLi ^iili — rnj) = 0. From Eq. (I4.13P one gets 



Vz Xi{qf - mf 



(4-21) 



together with 




(4.22) 



The first condition of (I4.2ip comes from the second equation and third equation (with 
Zi = Xi) of (I4.13p . The second condition of fl4.2ip comes from the third equation of 
(14. 13p with Zi = q. One notices immediately that Eq. (I4.2ip can be obtained from 
Eq. (I4.17P by replacing complex ctj with real Xj. Conditions (I4.22p come from the 
definition of the physical hypercontours. The Landau singularities may occur in 
the physical region if equations (14.2 ip and (I4.22p are satisfied. These are necessary 



*Landau has used this representation to devise the condition for singularities [30j 
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conditions since we cannot be sure that the real hypercontours are pinched when 
e — * (remember the second example of the previous subsection). 

If all strictly positive then we have a leading Landau singularity (LLS). 

Otherwise one has conditions for sub-leading Landau singularities (sub-LLS). 



4.3 Necessary and sufficient conditions for Landau 
singularities 

It may be shown that if Landau matrix Qij (defined in Eq. (14.161) . see also Eq. (14.251) ) 
has only one zero eigenvalue then the necessary and sufficient conditions for the 
appearance of a singularity in the physical region are equations (I4.2ip and (I4.22p . 

This important conclusion has been pointed out in the paper of Coleman and 
Norton [31] □. The proof is very simple and will be given in the next subsection (after 
equation (14.431) ). It is based on the underlying fact that the Landau matrix Qij is 
real and symmetric hence can be diagonalized by a real orthogonal co-ordinate trans- 
formation. It means that for unstable particles with complex masses the argument 
fails and the conditions based on Landau equations are no longer sufficient. 

We seek conditions for Eq. (14.211) to have a solution Xi = for i = M + 1, . . . , N 
with 1 < M < and Xi > for every i e {1, . . . , M}. The Eq. ( K^ becomes 

Xi = for i = M + 1,...,N, 

qf = m\ for i = 1, . . . , M, (4.23) 

EM ri 
i=i Xiqi = 0. 

For M = N one has a leading singularity, otherwise if M < this is a sub-leading 
singularity. Multiplying the third equation in (14.231) by qj leads to a system of M 



^See aslo Itzykson and Zuber ,98J in p. 306. 
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equations 



QuXi + Q12X2 H QimXm = 0, 

Q21X1 + Q22X2 + ■ ■ ■ Q2MXM = 0, 



QmiXi + Qm2X2 + ■ ■ ■ QmmXm — 0, 

where the Q matrix is defined as 



(4.24) 



Qij = 2qi.qj = ml + m] - {qi - qjf = m] ^ m] - (r^ - r^f, i,j = 1,..., M,(4.25) 



which agrees with Eq. f l4.16p . The necessary and sufficient conditions for the appear- 
ance of a singularity in the physical region now become 



det(Q) = 
Xi> 

2 2 

qf = m\ 
Qi = Qi 



(4.26) 



for 2 = 1, . . . , M. 

The condition det(Q) = defines a singular surface or a Landau curve. 

If some internal (external) particles are massless like in the case of six photon scat- 
tering, then some Qtj are zero, the above conditions can be easily checked. However, 
if the internal particles are massive then it is difficult to check the second condition 
explicitly, especially if M is large. In this case, we can rewrite the second condition 
as following 



det(4A/)/det(Q 



MM, 



>0, j = l. 



M- 1, 



(4.27) 



where Qij is obtained from Q by discarding row i and column j from Q and det((5jAf) = 
d[det{Q)]/{2dQjM), det(QMM) = d[det{Q)]/dQMM- The proof for Eq. (1^271) is the 
following. It is obvious that when the condition det{Q) = is satisfied one can set 
XAf = 1 and discard the last equation in (14.241) . After moving the M-column from the 
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left hand side to the right hand side, one obtains a system of M — 1 equations with 
M — 1 variables. The solution of this is clearly equation (14.271) . If det (Qmm) = 
then condition (14.271) becomes det{QjM) = with j = 1, . . . , M — 1. 

Conditions (I4.2ip and (14.221) admit a beautiful physical interpretation. This was 
discovered by Coleman and Norton [31]. Consider the case where all Xi are strictly 
positive. All the internal loop particles are therefore on-shell and have real momenta. 
An internal particle has a real four-momentum: qi = niiUi (for each i) with Ui is a 
four- velocity. Each vertex can be regarded as a real space-time point. The space-time 
separation between two vertices reads 

dXi = dTiUi = — -Qi, for each i, (4.28) 

TTli 

where dTi is the proper time. Following a closed loop, one has dXi = 0. Comparing 
this to the second equation of (I4.2ip we get the correspondence: dTi = rriiXi > for 
each i. It means that the loop particle is moving forward in time. dXi can be positive 
or negative depending on the sign of g,. If one chooses a reference frame where vertices 
are ordered in time, i.e. dX^ > 0, then > in that frame. This information can 
be very useful in practice. Let us illustrate this point. Consider two important 
Feynman diagrams in Fig. 14.31 We choose a reference frame where the arrows of the 




Figure 4.3: Typical triangle and box Feynman diagrams. 



internal lines follow the time direction. We look at vertex 2 (connected to P2) of the 

^One can regard different particles running in a loop as different states of one particle. Each 
vertex is associated with an external " force" . 
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triangle diagram and choose a co-ordinate system such that qi = (mi, 0,0,0), q2 = 
(e2, q2x, 0, 0). With pf = Mf, from the energy-momentum conservation P2 = Qi — <i2 
we get 

m? + ma — Mo 
62 = - 



2mi 

2 \{Mlmlml) [M| - (mi - m2)^][M| - (mi + m^f] 

- M " • ^ 

From the conditions q2 = q2 and 62 > we get 

M| < (mi - m2)^ (4.30) 

Similarly, for other vertices of the triangle diagram we have 

Ml > (mi + m4)^ M| > (m2 + mg)^ (4.31) 

By using the same trick one can easily see that necessary conditions to have a leading 
Landau singularity in the box diagram are 

Mf > (mi + m4)2, M| > (m2 + mg)^ 

(4.32) 

M| < (mi - m2)^ M| < (mg - m^f. 

Similarly, with t = {p2 + p-if and u = {p^ + p^Y and using the energy-momentum 
conservation, we get the constraint 

t > (mi + m^Y and u > (m2 + m4)^. (4.33) 

Thus a necessary condition for any diagram to have a leading Landau singularity is 
that it has at least two cuts which can produce physical on-shell particles 0. Other 
external particle which does not correspond to those cuts must have mass smaller 
than the difference of the two internal masses associated with it. 



^In the case of two-point function, the two cuts coincide due to energy-momentum conservation. 
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4.4 Nature of Landau singularities 
4.4.1 Nature of leading Landau singularities 

Our purpose is to extract the LLS by using Feynman parameter representation fl4.15p . 
The matrix Q which appears in the denominator is real and symmetric hence can be 
diagonahzed by a real orthogonal co-ordinate transformation. In general, Q has N 
real eigenvalues called Ai, . . . , A^v- The characteristic equation of Q is given by 

/(A) = A^ + (-l)a^_iA^'i + (-l)V-2A^-'-...(-l)'^~VA +(-!)% 

= (A-Ai)(A-A2)...(A-A„) = 0. (4.34) 

For the case = 4 we have 

Co = A1A2A3A4 = det(Q4), 

0-1 = A1A2A3 + A1A2A4 + A1A3A4 + A2A3A4, 

02 = A1A2 + A1A3 + A1A4 + A2A3 + A2A4 + A3A4 = ^[Tr(g4)'-Tr(g2)], 

03 = Ai + A2 + A3 + A4 = Tr(g4), (4.35) 

Consider the case where Q has only one very small eigenvalue A^r -C 1. Then, to 
leading order 

Atv = — , ai = A1A2 . . . Aa^-i 7^ 0. (4.36) 
Let V = {x?, X2, . . . , x%} be the eigenvector corresponding to Aat. is normalised to 

N 

E^° = l- (4.37) 

i=l 

For later use, we define 



v'^ = V.V. (4.38) 
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The expansion of A around V reads 

N N 



i,j=l i=l 
1 ^ 1 

where i/i = Xj — For the leading part of the singularity it is sufficient to neglect 
the linear terms. The Q-matrix can be diagonalised by rotating the y-vector 

N 

y, = Y,Aj^j, (4-40) 
i=i 

where A is an orthogonal matrix whose columns are the normalised eigenvectors of 
Q. Thus we have 

j = l 

N-l 
i=l 

Note that the term Xnz% in the rhs has been neglected as this term would give a 
contribution of the order 0{X%) to the final result. Equation (14.151) can now be 
re-written in the form 

Although the original integration contour is some segment around the singular point 
Zi = with i = 1, . . . ,N, the singular part will not be changed if we extend the 
integration contour to infinity, provided the power (A^ — D/2) of the denominator in 
Eq. (I4.42P is sufficiently large. Integrating over zn gives 

{—1)^T{N — D/2)v f^°^ 1 

where the factor v comes from the ^-function. One sees clearly that each integration 
contour is pinched when e ^ if all Aj 7^ with z = l,...,A^ — 1. 
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Asumming that Aj > for i = 1, . . . , K and < for j = K + 1, . . . , N — 1 with 
< i^' < — 1. We then change the integration variables as follows 

U = y/Xizi for i = 1, . . .K, 

(4.44) 

tj = ^/^jZj ioT j = K + 1, . . . N - 1. 
This makes sure that all are real. We get 

_ {-irm-D/2)v 

J- n 



■0 



7r^/223^/2-^^(-l)^-^-iai 

+00 f + OO 



dti-'-dtx I dtx+i - ■ ■dtN_i jTj— j ,(4.45) 

7-00 ^ i-ElK+ltl + b')''-''/' 

where 

K 

= J2 ^1 + ^nv^ - ie, Re{b^) > 0. (4.46) 

i=l 

Changing to spherical coordinates by using formulae ( ID. 181) we get 



T 



^ {-1)^T{N - D/2)v 27r(^-^-i)/2 



^D/223D/2-7V^(_l)iV-A'-l^^ T{{N -K- l)/2) 

_ rfti ■ ■ ■ dtK I drj^—^^^^^. (4.47) 
Note that (6^ - r'^)N-D/2 ^ ^-i.{N -d 12) _ ^2)JV-d/2 ^j^g f^^^^ ^j^^^^ ^ > 0. 

Then by using formula (ID. 190 we have 

(_\\Ni-K{N-K-\)l2 

= ^ ^ ^ ^iN~K~l)/2r((M-n+K+U/9.) 

7rD/223D/2-7V^(_l)JV-/r-la^ 

r>+CO -j^ 

dti---dtK T? . (4.48) 

(E.=l ^? + _ (iV-D+X+l)/2 ^ ' 

Repeat the above steps to get 

^ (-l)^e-(^-^-^/^t; (47r)(^-^+^)/^r((Ar - D + l)/2) 

2(3+iV)/2^^(_l)iV-i.-l^^ _ ^^ (iV-D+l)/2 • I • 



This result holds provided 

ai ^ and - D + 1 > 0. (4.50) 
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In the case where N — D + 1 < one can write D = 4 — 2e {e > 0) and do the 
expansion in e. Apart from a divergent term of the form 1/e related to the artificial 
infinite boundary, the other terms give the nature of singularities. 
A similar result for the nature of the singularity has been derived in |106j in the 
general case of a multi-loop diagram including the behaviour of the sub-LLS. The 
extraction of the overall, regular, factor which is the i^'-dependent part in Eq. (14.491) 
(see also Eq. (14.721) for the sub-LLS) is more transparent in our derivation. 
For = 4, D = 4 

(To^).» = e-(='-^)/2 ^— 



= e'-(^-^)/^^= ^ (4.51) 
4v/(-l)3-^ao-2e 

For A^ = 3, D = 4 - 2e, we use T{e) = (1/e) - -fE to get 

ini2-K)/2 

87rv/(-l)2-^AiA2 
For AT = 2, D = 4 - 2e, we use r(-l/2) = -20F to get 

{nU = --^{\2v'-^eY'\ (4.53) 
Sttv Ai 

For A^ = 1, D = 4 - 2e, we have 



_ -r(-l+e)(47r)% 
Ai fl 



IGtt^ V e 



-7^ + ln(47r)-ln(Ai) , (4.54) 



with Ai = m 



2 



Remarks: The leading Landau singularity appears when Aat ^ 0. The nature of 
the leading singularities for the scalar one-, two-, three-, four- functions are 1, 1/2, 
log, —1/2 respectively. One remarks that in the case A^ = 4, 3 the LLS is divergent, 
i.e. becomes infinite. The LLS is finite but singular, i.e. the derivative is divergent 
at the singular point, in the case N = 2 and is regular in the case A^ = 1. The 
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scalar three-point function and its square are integrable at the LLS point. The scalar 
four-point function is also integrable at the LLS point but its square is not. 

One may wonder if we can use the general result in Eq. (14.49 1) for the case > 5. 
The answer is YES as long as ai 7^ 0. As will be proved in the next subsection, ai is 
proportional to the Gram determinant det(G') at the singular point. Since det(G) = 
for > 6 in four dimensional space, we conclude that Eq. (14.491) cannot be used for 
the case A^ > 6. However, the LLS can occur for N = 5. If this happens, we will have 
five on-shell equations qf = mf with i = 1, . . . , 5 to solve for q^. We just need four 
equations to find g^, the rest is a 5-function to give some constraint on the internal 
masses and external momenta. Thus the nature of 5-point function LLS is a pole 
[29] . Indeed, it is highly nontrivial to find a physical process which contains a 5-point 
function LLS. 



4.4.2 Nature of sub-LLS 

In order to understand the nature of sub-leading Landau singularities, one should 
integrate over xn from Eq. (14.151) . This gives 

_ {-irnN - D/2) vii- 



= ^ \ ' ' / dxi--- dxN-i-^ — , (4.55) 

where rj is Heaviside step function and 

A(xi, . . . ,Xiv^i) = A(a;i, . . . ,a;Ar_i, 1 - xi - . . . -X7v_i) 

^ N-1 N-1 ^ 

i,j=l i=l 

where 

Gij = Qij - QiN - QjN + Qnn = '^ri.Tj, Pi = Qnn - QiN = m% - + r'j. (4.57) 
Thus det(G) is just the Gram determinant. From Eq. (I4.57P we get 

N-1 

det(Q) = Qnn det(G) - ^ PiP.Gi,. (4.58) 
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The Landau equations for representation (I4.55P arc 



A = 0, 

X, = 0, i = l,2,...u, (4.59) 
11 = 0, ^ = z. + l,...iV-l. 



The third equation of fl4.59p gives 



N-l 

GijXj = f3i, i = u+l,...N -1. (4.60) 



If det(G') 7^ then the solution reads 



N-l , N-l 



= E = E (^Av t = u + l,...N-l. (4.61) 

Thus the solution of the second and third equations of fl4.59p is Xi = Xi with 

X = (0, . . . , 0, Xy+i, . . . , XAr_i). (4.62) 



The first equation of fl4.59p gives the equation of the surface of singularity [107] 

I iv-i 



(1 - 

Ax = 



2 2 ^ 

where we have used Eqs. f l4.6ip and (14.580 . Not surprisingly, one obtains again 
det((5) = 0. In the case det(G') = 0, the condition for the second and third equations 
of fl4.59p to have solution is Yl!j=v+i(^fiij ~ ^ (see Eq. fl4.6ip ). This together with 
Eq. fH38l) give det(Q) = 0. 



^It is important to notice that when one performs the x^r-integration the boundaries become: 
Xi — for i = 1, . . . — 1 and 1 — X^i^i^ ~ 0. In the mean time A becomes inhomogeneous, so 
that the first equation of (|4.59p is not automatically satisfied when the others are. 
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In the neighbourhood of a point x that hes on the surface of singularity, we expand 
A, keeping only the lowest terms: 



i=l 



dA 

dxj 



N-l 



i,j=u-{-l 



dxidxj 



N-l 

+ ^ yiVjCi 



(4.64) 



with 



, ^ dA 



i=l 



— ze, 



Xi — Xi, i = u + 1, . . . , N — 1. 



(4.65) 



Integral (14.551) becomes 



-^0 



-iyT{N-D/2) 

^D/223D/2~N 
+00 N-l 



1 



n 



i=u+l 



(4.66) 



where, similar to the calculation in the previous subsection, we have let each yi- 
integration contour run from — oo to +oo, provided the power {N — D/2) of the 
denominator is sufficiently large. To understand the difference between the LLS and 
sub-LLS we should compare Eq. (I4.39P to Eq. (I4.64p . One remarks that the linear 
terms only appear in the case of sub-LLS. The ?/j-integration is exactly the same 
for the two cases. Gij is a real symmetric matrix hence can be diagonalized by a 
real orthogonal co-ordinate transformation. Using the same method described in the 
previous subsection we integrate over yi to get 



N 



(N-D-v-l)/2 



23D/2-N ^(_l)^-.-l-Xdet(G) 

T{{N-D + iy + l)/2) / r[da;,[C(a;0] 

^0 i=l 



-{N-D+u+l)/2 



(4.67) 
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where K is the number of positive eigenvalues of Gram matrix Gij. Of course, one 



can recover Eq. fl4.49p by setting u = 0. Asumming that h = §^ 
leading part of the singularity occurred when Xi 0"*" we have 



7^ 0, for the 



/ l\dx,[C{x,)f = 2"/" f\dx,[b,x, + A{x)] 
Jo ^0 

r(a + l) 



where a = -{N -D + u+l)/2. With 7 = -a - z/ = (iV - z/ - D + l)/2 we get 



{Tq )div 



zTX 



{N-D-u-l)/2 



X 



23D/2^N^^_l^N-u~l-K det(G') 

2-- (_i)-r(7 + z/)r(-7-z/ + i) 



UUb, r(-7 + i) 

We then make use of the following identity 



[2A(x) - ie] 



(4.69) 



r(i-z)r(z) 



TT 



SmLZTT 



(4.70) 



to get 



The final result reads ^ 



-l)"r(7 + z/)r(-7 -;/+!) 

r(-7 + i) 



r(7). 



(47r)^r(7) 



det(Q) 
2det(G) 



(4.71) 



~7 

(4.72) 



2(3+^-)/%v^(-l)^--^-i det(G') Wi=ibi 
which holds provided 

dA 

det(G)^0, bi = -rf- ^0 and 7 = (A^ - z/ - D + l)/2 > 0. (4.73) 
If some bi = then expansion fl4.64p must be modified to take into account the higher 



order terms xAx^ — xA-J^-^ 



with i 



V. If 7 < one can write D = 4: — 2e 



{e > 0) and do the expansion in e to get the divergent term independent of e as in 



^Similar result has been obtained by Polkinghorne and Screaton |106j . 
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the case of LLS. 

Comparing Eq. (14.491) to Eq. f l4.72p we see that, apart from the finite factor related 
to bi, the nature of the singlarities given by the latter can be obtained from the 
former by simply replacing N hj N — u which are the number of on-shell internal 
particles. We remark also that there must be some relation between ai and det(G') 
at the singular point, namely det(G') oc ai when det{Q) =0. As a consequence, if 
ai = det(Q) = (Landau matrix Q has at least two zero-eigenvalues) then det(G') = 
(Gram matrix G has at least one zero eigenvalue). There is a beautiful way to prove 
this mathematically ^. Let Va = {x^""^} with a = 1,2 and i = 1, N he two linearly 
independent vectors of the degenerate zero-eigenvalue lj. We can always normalise 
Va such that 

N 

J]a;S'^) = l, a =1,2. (4.74) 



One has 



i=l 



N N-l N-l 



j=i j=i j=i 

N-l 

= ^^{Qij ~ QiN)x^f' + QiN 

i=i 

N-l N-l N-l 

= ^^(Qii ~ QiN — QjN + QNN)x^j'^ + QnjX^j"^ — QnN X^j"^ + QiN 

j=l j=l j=l 

N-l 

= Y,Gijxf-f3., (4.75) 
where relations (I4.57P have been used. Thus one gets 

N-l 

J2Gi,xf=Pi, a = 1,2. (4.76) 



-"^^We have learnt this trick from Eric Pilon, many thanks! 

^^In the case of the real symmetric matrix Q, the degree of degeneracy of one zero-eigenvalue is 
equal to the number of zero-eigenvalues. 
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With Vo = Vi - 1/2 7^ this gives 

G.Vo = 0. (4.77) 
This means that the Gram matrix has at least one zero-eigenvalue hence det(G') = 0. 

4.5 Conditions for leading Landau singularities to 
terminate 



(b) 





(a) 





• 









(c) 




• 






• 





Figure 4.4: Mechanism for termination of LLS in Xi-plane. 



This section concerns the termination of LLS as we vary the external parameters 
denoted by Mj (without any loss of generality, we assume that the internal masses 
are fixed for the sake of simplicity). 

It is obvious that the position of LLS and its properties depend on the values 
of Mj. If we vary Mi continuously, while maintaining the pinch conditions, the only 
mechanism for the termination of a LLS is the followinj^. The LLS moves to the end 
of the integration contour (xj = 0). Thus the LLS will coincide with a sub-LLS and 
move off the physical region afterwards |108[ 11091 1110] . This is illustrated in Fig. 14.41 

The following remark will be useful later. The 4-point LLS terminates when it 
coincides with a 3-point sub-LLS which in turn will terminate when coinciding with 



^^In the case of more than one loop, there is another mechanism. The LLS can terminate if two 
pinches meet on the integration contour. If this happens, the singularity may somehow leave the 



physical region [I08lfT09]. 
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a normal threshold. The normal threshold coincides with itself, i.e. it occurs at one 
point. 

A good question to ask is " How does a LLS terminate?" This is a mathematically 
complicated question and we do not attempt to give a complete answer. What we 
understand is as follows. When moving from Fig. 14.4b . to Fig. 14.4b . the leading Landau 
curve (det((5Ar) = 0) changes continuously until it makes an effective intersection with 
a sub-leading Landau curve (det(Q7v-i) = 0) |1U81 IllUi I1U9] . At the point of contact, 
both curves have the same slope and both corresponding Landau equations have the 
same solution of Xj \WS\ I109j . At effective intersections the nature of the Landau 
singularity may change [1101 1109] . In Fig. 14.4b . we have a A^-point LLS whose nature 
is given by 

Ta{Mi) = Aa{M,)+Ba{M,)fa{detQN), (4.78) 

where the functions A and B are analytic in a neighborhood of the singular point Mi, 
the function /a(det Qn) is singular at this point. fa{z) can be z^^'^, Inz or z~^/'^ if A^ 
is 2, 3 or 4 respectively. For Fig. 14.4b we have a similar equation 

T,(M,) = A,{M,) + 5,(M,)/,(detQ^_i), (4.79) 

where we have assumed that Fig. 14.4b has a (A^ — l)-point Landau singularity. The 
nature of the coincident singularity in Fig. 14.4b is a product of two factors which are 
similar to Tq and Tc |110j . Thus, we have 

n{M,) = A + Bfaidet Qn) + Cf,{det Qjv-i) + Dfa{det QN)fc{det Qn^i). (4.80) 

li D ^ then the leading singularity is given by the last term which means that the 
Landau singularity can be enhanced at termination point. This kind of enhancement 
can be somehow understood if we look at some formulae in this thesis: from Eq. (14.721) 
we see that if a leading A^-point singularity coincides with a sub-LLS then bi = 
1^ =0 leading to an enhancement from the prefactor; from Eq. (14.1011) we 

' x=x 

observe a product of two singularities (a leading Landau singularity and a collinear 
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divergence); from Eq. (lE.lSp we see the possibility that a product of two singularities 
can occur if the integrals related to 3-point functions develop a pinch singularity at 
the same position where the prefactor related to the leading Landau determinant 
(1/ det((54)) vanishes. 

If a four-point LLS coincides with a three-point sub-LLS, the nature of this sin- 
gularity can be z~^^'^\nz which is integrable but its square is not. 



4.6 Special solutions of Landau equations 
4.6.1 Infrared and coUinear divergences 

In this section, we will show that infrared and collinear singularities are solutions 
of Landau equations. However, in order for them to become divergent additional 
conditions must be satisfied. As one might anticipate from Eq. fl4.72l) sub-leading 
Landau singularities can be enhanced by various factors. 

Infrared divergence 

We consider the case of a sub-LLS where xi = . . . = xjv-i = and xat > 0. The 
Landau equations become 

q% = m\ and = 0. (4-81) 

We get m^v = 0. As remarked in subsection 14. 4[ for the case = 1,2 the Landau 
singularities are finite hence there is no infrared divergence in those cases. We thus 
consider the case A^ = 3. With z/ = 2, equation fl4.66p becomes 

f dxidx2-^^, (4.82) 

with 

2 2 

C{xi) = 2ml + 2 X] ^if^i + XI '^ij^i^j ~ (4-83) 

1=1 i,j=^ 
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If Pi = mg — + rf ^ then one can neglect the quadratic terms in C{xi) to get 



(r„^)..„ ~ J^^ . _,„^ ^ 0. (4.84) 

If I3i = m\ — mf + rf = Q then we get 

iT^)div ~ / dxidx2-—: — — Inlmg) oo. (4.85) 



0+ 



2m| + X]i,j=i GijXiXj 



The nature of Landau singularity is m| if /?j 7^ and is enhanced to ln(m|) if (3i = 
/?2 = 0. The condtions to have an infrared divergence for the case of three-point 
function therefore are 

m-s = 0, m\ = r\ = p\, m\ = r\= p\. (4.86) 

For the cases > 3 we can always reduce them to three-point functions hence we get 
the same conclusion. Physically, one sees that conditions fl4.86p can be satisfied only 
by the photon or gluon. For the electroweak theory, if we take the limit M\y 
then the one-loop diagrams involving the W-gauge boson as an internal particle have 
no infrared divergence since it couples to particles with different masses. 

Collinear divergence 

For M = 2, i.e. xi^2 > and X3 = . . . = xn = equations (14.231) become 

3^1,2 > 0, 

qf = ml, q2 = ml, (4.87) 
XiQi + X2q2 = 0. 

One gets Xinii = ^2^2. If mi 2 7^ then Xi{qi + ^q2) = whose solution corresponds 
to the normal threshold = (gi — ^2)^ = {fni + ^2)^. If mi^2 = one gets 

Ximl + X2{qi.q2) = 0, (4.88) 

which gives qi.q2 = or = {qi — ^2)^ = 0. This solution corresponds to a collinear 
divergence whose nature is also logarithmic [lllj. Clearly, this collinear divergence 
cannot occur if = 1. 
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It is important to remark that the solutions for infrared and colhnear divergences 
appear in the hmit of massless internal particles. These solutions require no constraint 
on relevant Feynman parameters, even the positive condition is not necessary. 

4.6.2 Double parton scattering singularity 



Figure 4.5: A typical box Feynman diagram which has a double parton scattering 
singularity. 

There exists a special case of Landau singularity called double parton scattering 
(DPS) singularity [104[ [22] which appears also in the massless limit. Unlike the 
sub-leading infrared and collinear divergences, the DPS singularity is a LLS and its 
solution requires some sort of constraint on relevant Feynman parameters (the positive 
sign condition is important). 

Let us consider the case of g{pi)g{p2) W {p3)W {p4) box diagram displayed in 
Fig. 14.51 The internal particles w-quark and (i-quark are massless. The Q-matrix is 
given by 



Pi 



g 




g 











-t 




Qa 












—u 



(4.89) 



-t 












V 








/ 
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where t = {pi — ps)^, u = {p2 — Ps)^- The Landau determinant, characteristic poly- 
nomial and Gram determinant read 

det(Q4) = {tu-M^)^ = 0, 

det(G3) = 2s{tu-M^), (4.90) 

with s = (pi+P2)^- One sees that ai = and det(Q4) = det(G'3) = at the boundary 
of phase space where the Landau matrix has two zero-eigenvalues. The phase space 
is defined as 

s > 4:M^, tu-M^>0 with t + u = 2M^ -s < -2M^. (4.91) 

In this special case where the Landau matrix has two zero-eigenvalues at the singu- 
lar point, the conditions given in subsection 14.31 (see Eq. fl4.26p ) are necessary but 
no longer sufficient. One should keep in mind that this box diagram always has a 
collinear divergence associated with the reduced two-point functions as discussed in 
the previous subsection. The necessary conditions for a LLS read 

tu = ML, 

^ (4.92) 

t < 0, ?i < 0, 

which are compatible with Eq. (I4.9ip . We now have two questions to be answered: 
whether they are sufficient conditions for a pinch singularity? If Yes then what is the 
nature of the singularity? 

For these questions, we come back to equation (14.151) and perform the following 
nonlinear change of variables [1071 1104] 

Xi = era, X2 = a(l — a), X3 = r/3, a;4 = r(l— /5). (4.93) 

The inverse solution gives us the range of new variables 

0<(T,r<oo; 0<a,/3<l. (4.94) 
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The Jacobian is simply ar. With D = 4 + 2e {e > 0) equation fl4.15l) factorizes 



1 1 



The first integral is just the Beta function B{e,e) producing a collinear divergence. 
We use 

i?(^,^) = ^^^^ = ^ + (^(^) (4.96) 



to get 



Tq = — -- / dad(3 



■° 87r% Jo + - + /3) - ?i - ze]2 

+ ^ rc^ac^/j ln[.a/3+(.-M^)(a + /3)-.] 

We are interested in the collinear divergent term. The relevant integral reads 

1 1 



(M^ - t) {Ml - u) io (« - ^) (« - a) ' ^^-^^^ 



where we have integrated over (3 and 



M2 . e 



— U 6 

This is nothing but the second example fl4.3p . From Eq. (14. 7p we have the necessary 
and sufficient conditions for a pinch singularity: 

< w = a < 1 and (M^ - t)(M^ - m) > 0, (4.100) 

which are completely equivalent to Eq. (I4.92p . From Eq. (14.41) we get the following 
result at the singular point tu = 

1 2m 



(I 



1 )div 



i 1 i 1 

^^'^''"^ ^ 47r(M^-t«)i + --- = "4vrv/d^e+--- " ^^'^^^^ 
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We conclude that in this special case where the Landau matrix has two zero-eigenvalues, 
conditions (14.921) given by the Landau equations are necessary and sufficient for the 



appearance of a LLS. The nature of this singularity is 1/ ^/deU^ij which is consistent 
with Eq. (14.491) . The LLS goes together with the collinear divergence (see also section 
. This LLS is called double parton scattering singularity ffist pointed out in |"lU4j . 



We notice that Eg. (I4.10ip disagrees with the result of Ellis and Zanderighi given 
in Eq. (4.21) of |112)^^l where they claimed that, apart from the collinear divergence. 



a finite result can be obtained by doing an expansion around the LLS point tu = M^. 

This DPS singularity is not integrable. In a practical calculation, the cross section 
must be finite. This singularity must somehow disappear order by order. For the 
process gg —>■ WW there is no tree level diagram hence there is no real radiation at 
the leading order. One infers that the numerator must vanish at the singular point. 
The actual calculations |113[ I114j have confirmed this. We have conjectured it as a 
consequence of the gauge dynamics [22] • To understand this dynamical cancellation, a 
physical investigation is necessary. The same phenomenon of DPS singularity happens 
in the case of six photon scattering whose numerical analysis around the singular point 
is given in [22] . 
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Many thanks to Denner for pointing out this reference to us. 



Chapter 5 



SM bbH production at the LHC: 

Mh > 2Mw 




The content of this chapter is based on our pubhcations |115l [77] . 

An important notation used in this chapter should be clarified. The notation of 
center-of-mass energy, y/s, can be used for the sub-process gg bbH or pp — > bbH. 
Which use will depend on the context, as you will notice very easily. 



5.1 Motivation 

The aim of this chapter is to extend the study we made in chapter [3] to higher Higgs 
masses. 

When trying to do so for the case > 2My/ we have encountered severe nu- 
merical instabilities for the cross section involving the one-loop amplitude squared. 
We have tried to understand this problem and found the following facts. At the 
level of NLO which involves the interference term between the tree-level and one-loop 
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amplitudes no instability was present. On close inspection it was found that the in- 
stabilities were only due to the box and pentagon diagrams of class (c) of Fig. 13.41 
The contribution of various triangle diagrams does not show any numerical instabil- 
ities. At the partonic gluon-gluon level it was found that there is no instability for 
y/s < 2mt and that independently of Mh and the result was completely stable 
for irit = Mw- These threshold conditions were a sign for the possible existence of a 
leading Landau singularity for the box diagrams whose square is not integrable. The 
five point functions are reduced to four point functions hence should have the same 
problem. Indeed, some triangle diagrams of class (c) of Fig. 13.41 have also LLS whose 
nature is integrable. That's why they do not show any numerical instability. 

In order to solve this problem of Landau singularities, we first have to understand 
them in detail by applying the general analysis of the previous chapter to the specific 
problem of gg bbH. 

5.2 Landau singularities in ^ bbH 

In this section, we discuss all Landau singularities that occur in the Feynman diagrams 
of class (c) of Fig. 13. 4[ It is pedagogical to start with the three point function. 

5.2.1 Three point function 

As concluded in section 14.31 a necessary condition for a three point function to have a 
LLS is that it has two cuts which can produce physical on-shell particles. The other 
external particle which does not correspond to those cuts must have mass smaller 
than the difference of the two internal masses associated with it. The diagram in 
Fig. 15.11 satisfies those conditions. The scalar three point function of this diagram 
reads 

T^{S2) = Co{s2, Ml, 0, ml M^, M^) (5.1) 
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Pi 



P2 



b, p. 



/ N 

/ Xw Xw \ 
> < 1—* 




b, P, 



Figure 5.1: Left: A triangle diagram contributing to gg —>■ bbH that can develop a 
leading Landau singularity for Mh > 'iMy/ and ^/s2 > Trtt + My/, i.e. all the three 
particles in the loop can he simultaneously on-shell. Right: the three point function of 
the left diagram where the arrows indicate momentum flow. 



with S2 = {p4 + p^Y and the bottom-quark mass has been neglected. Necessary- 
conditions to have LLS are 

(5.2) 



(5.3) 



Mh > 2Mw and ^/s^ > rrit + Mw 
The phase space is defined by 

Mh<S2< S. 

The characteristic equation writes 

— + a2\^ — ai\ + ao = 

with 



(5.4) 



ao 
ai 

02 



< 



det(g3) = -2M^s^ + 2Mfj{mi + M^)s2 - 2Mfj[Mfjmi + {mf - M'^ 

- [Ml{Ml - AM^) + {S2 - m] - M^f + mt + - Qm^.M^] 

- [mUmI - 4M^) + {ml - Mlf] < 0, 



2 n2i 



2{mi + 2M^), (5.5) 

where we have used conditions (15. 2p to prove that oi is negative. The equation 
det((53) = has two roots 
1 



MUmIt + m?) ± {ml - MI)Mh\ MI - 



(5.6) 
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The sign condition (xj > 0), Eq. fl4.27l) . is very simple. For this, we need 

det(Q33) = -Mi(Mi-4M^) <0, 
det(Qi3) = -M^ (m? + M^) + 2s2M^ < 0, 

det(g23) = -M2(m,2 + M^)+S2(M^-2M^)<0, (5.7) 
which together with Eq. fl5.2l) give 

Only the minus solution given in Eq. (15.61) satisfies this, thus 
1 



^LLS 



2M^ 



MjjiM^ + m^) - (m? - M^)Mh^ Mj^ - AM, 



(5.9) 



This result tells us many things. First, if Mh < 2Mw there is no LLS. Second, if 
Mh = 2My/ the LLS occurs at $2^^ = 2{mf + M^) which is the maximum value 
of S2^^ given by Eq. (15. 8p . If Mh increases then s^'^'^ becomes smaller and smaller. 
The maximum value of Mh is determined when Sj^"^ reaches the normal threshold 
(this condition for the termination of a LLS was discussed in section : 

4M^ < Mi < 4M^ + ^{mt - Mw)\ 

rrit 

{mt + Mwf <S2<2{ml + Ml,). (5.10) 



Numerically, we have 



160.7532GeV < Mh < 172.889GeV, 

254.3766GeV < < 27L059GeV. (5.11) 



If one keeps increasing Mh > 172.889GeV, S2 will increase from its minimum value 
254.3766GeV and become larger than the limit j^tz^^I-M^ (see Eq. (15.81) ) required 
by the sign condition (xj > 0) hence be outside the physical region. In the mean time, 
one should notice that the normal threshold is moving towards the left boundary of 
phase space as Mh increases. When Mh = rrit + Mw = 254.3766GeV, the normal 
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Figure 5.2: Left: the real part of Cq as functions of y/s2 with various values of Mh- 
Right: the same plots for the imaginary part. 



threshold is at the boundary and disappears if Mh passes that value. The function 
Tq then has no structure. All this phenomenon is shown in Fig. 15.21 

We would like to explain the behaviour of the real and imaginary parts of Cq at 
the LLS point. Since A3 0, we have A1A2 = ai < as proved in Eq. (15. 5p . Thus 
there is one positive eigenvalue, i.e. K = 1. From Eq. (14.521) we get 



OJdiv 



ln(A3t; — ie) 



{i In lAat; | + vr) 



(5.12) 



Thus, one observes a positive jump in the real part and a logarithmic singularity 
in the imaginary part as described in Fig. 15.21 This singularity is integrable at the 
NLO as well as one-loop amplitude square level hence does not cause any numerical 
instability. 



5.2.2 Four point function 

The 5-point functions have no LLS but contain a sub-leading Landau singularity which 
is exactly the same as the LLS of the box diagram in Fig. 15.31 It is thus enough to 
study the singularity structure of this box diagram. With si = {ps + ps)^, S2 = 
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Ps 

"h 



Figure 5.3: A box diagram contributing to gg bbH that can develop a Landau 
singularity for Mh > 2Mw and y/s > 2mt, i.e. all the four particles in the loop can 
be simultaneously on-shell. The arrows indicate momentum flow. Internal momenta 
qi are real, on-shell and have positive energies. 

{pi + Ps)^, and the on-shell conditions Pi = pi = 0; Pa = = = 0, pi = M]j, 
fixing s and Mh, the scalar box integral is a function of two variables si,2 

To^(si, S2) = Do{Ml, 0, s, 0, si, S2, M^, M^, ml ml). (5.13) 

The kinematically allowed region (phase space) is 

Mh<sx< s, Ml— < S2 < Mh + s - si, 

Sl 

where the latter is obtained by writing S2 = (p4 +^5)^ = + 264(65 
in the center-of-mass system of (ps + ps). 

Conditions for the opening of normal thresholds 

The condition ^ XiQi = with a;^ > of Eq. (14.231) can be re- written in the form 

X2q2 + xsqs = xiqi + Xiq4 (5.15) 

with all > as shown in Fig. 15. 3[ Indeed, there are other possibilities like x^q^ = 
Xiqi + X4q4 + X2q2 with all > but this will require mi, > (jrit + Mw) which is 




(5.14) 



- IP5I cos 6^45) 
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impossible in our case. Thus equation flS.lSp is the unique possibihty for a LLS. This 
condition of positive energy together with the real and on-shell condition (g^ = q*, 
qf = ml) give a beautiful physical picture as shown in Fig. 15.31 where internal and 
external momenta share the same physical properties. 

As already discussed in section the above picture gives the following necessary 
conditions for the appearance of a LLS: 



Mh > 2Mw and y/s>2mt, 

Si > {rrit + Mw)"^ and S2 > {rrit + Mw] 

mt > Mw, 



(5.16) 
(5.17) 
(5.18) 



where we have used the fact that the momenta of the bottom-quarks and the Higgs 
boson flow in the same positive direction to get the last constraint (if we consider 
the inverse process H bbgg where momenta of the bottom-quarks and the Higgs 
boson are in opposite directions then we get Mw > rrit which cannot be satisfied by 
experimental data.). They are conditions for the opening of 4 normal thresholds. 



Landau determinant 

The reduced matrix, S^^\ which is equivalent in this case to the Q matrix for studying 
the Landau singularity, is given by 



^4 



2M^-Mfj 



2M, 



w 
M2 - 



2Mwmt 
2Mwmt 



2M^- 



2M^ 



2Mwmt 
2Mwint 



2Mwmt 



2Mwmt 
1 



2mi- 
2rni 



2Mwmt 

2Mwmt 
2m^—s 
2mi 



\ 



(5.19) 
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The determinant reads 



det {Q4 



16M^m^det(S4) 



bs2 + c, 



a = X{si,mlM^) = [si-{mt + MwY][si-{mt-Mw)% 

b = 2{-sl{m',+M^) + s^[{m',+M'^y-{s-2mf){Ml-2M^)] 

c = s2(m2-M^)2 + 2M^sK + M^)si 



+ sM^[(s-4m2)(M^-4M2,)-4K + M, 



(5.20) 



We remark that a = corresponds to a normal threshold and defines the asymptotes 
of the Landau curve. 



The Landau determinant can be written in the following beautiful form 



det(Q4) 
A 

det Q3(M^,si) 



s'2 



b 



a(s2-S2) -— , 
det Q3(s,Si)detQ3(M|,si), 



2M2, 



det Qz{s,sx) 
2ml 



2M2, 



-M\ 



If 



4M^ 



2m\ 



2 \2 



(5.21) 



which is very useful when one knows that the LLS coincides with a three-point sub- 
LLS. If that happens the Eq. det{Q4) = has only one root S2 = S2. The fact 
that the discriminant of a Landau determinant can be written as a product of lower- 
order Landau determinants is known as the Jacobi ratio theorem for determinants 
[TTUlfTTU] . 



Discussion of singular structure 

We would like to understand the singular structure of the scalar 4-point function 
defined in Eq. f l5.13p . Namely, we will look at the behaviour of its real and imaginary 
parts as functions of Si and S2 while other parameters are fixed. We take = 
174 GeV, Mw = 80.3766 GeV, = 353 GeV and Mh = 165 GeV. 
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The behaviour of the Landau determinant, the real and imaginary parts of the 
4— point function Tq are displayed in Fig. 15.41 as function of si, S2 within the phase 
space defined by Eq. fl5.14p . We clearly see that the Landau determinant vanishes 
inside the phase space and leads to regions of severe instability in both the real 
and imaginary parts of the scalar integral. We notice that this region of Landau 
singularities is localised at the center of phase space. The boundary of singular 
region will be explained later, see Eqs. fl5.35p and fl5.36l) . 




Figure 5.4: The Landau determinant as a function of Si and S2 (upper figure). The 
real and imaginary parts of Dq as a function of Si and S2- 



To investigate the structure of the singularities in more detail let us fix = 
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^2(m2 + M^) ^ 271.06 GeV, so that the properties are studied for the single vari- 
able S2- In order to do so, we have to find all the reduced diagrams containing S2 as an 
external momentum squared. There are 3 diagrams shown in Fig. 15.51 one self-energy 
and two triangle diagrams. The plots for the real and imaginary parts are shown in 




Xj=X3=0 X3 = Xj=0 



Figure 5.5: Three reduced diagrams of the box diagram in Fig. 15.31 that contain S2 
as an external momentum squared. The self-energy diagram has a normal threshold. 
The two triangle diagrams contain anomalous thresholds. We refer to subsection \5.2.1\ 
for a detailed account of the 3-point Landau singularity. 



Fig. 15.61 This figure is very educative. We see that there are four discontinuities in 
the function representing the real part of the scalar integral in the variable a/s2: 



As 82 increases we first encounter a discontinuity at the normal threshold = 
rrit + Mw = 254.38 GeV, representing Hb — ^ Wt. This corresponds to the 
solution (for the Feynman parameters) Xi 3 = and ^2,4 > of the Landau 
equations (see Fig. 15.51 -left). 

The second discontinuity occurs at the anomalous threshold ^/s2 = 257.09 GeV 
of a reduced triangle diagram. This corresponds to the solution X3 = and 
3^1,2,4 > of the Landau equations (see Fig. 15.51 -middle). As discussed in sub- 
section 15.2.11 the singular position is given by 

= ^ (m^K + M^) - Mn^Ml-AM^im', - M^)^ (5.22) 

which gives = 257.09 GeV. 
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Figure 5.6: The imaginary, real parts of Dq and the Landau determinant as functions 
of Js^- 



• The third discontinuity corresponds to the diagram of Fig. 15.51 -right. The posi- 
tion of this Landau singularity is given by 

^2 = ^ [s{ml + M^) - ^s^s-Amliw?, - M^)) , (5.23) 

which gives ^/s2 = 259.58 GeV. 

• The last singular discontinuity is the leading Landau singularity. The condi- 
tion det(S'4) = for the box has two solutions which numerically correspond 
to ^/s2 = 263.88 GeV or = 279.18 GeV. Both values are inside the phase 
space, see Fig. 15. 6[ However after inspection of the corresponding sign condi- 
tion, only = 263.88 GeV (with xi ^ 0.53, ^ 0.75, ^ 0.77) qualifies 
as a leading Landau singularity, ^/s2 = 279.18 GeV has Xi ~ — 0.74,X2 ~ 
— 0.75,X3 ~ 1.07 and is outside the physical region. 
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The nature of the LLS in Fig. 15.61 can be extracted by using general formula fl4.5ip . 
With the input parameters given above, the Landau matrix has only one positive 
eigenvalue at the leading singular point, i.e. K = 1. The leading singularity behaves 
as 

Df" = ^ (5.24) 

WM^m^^detiSi) -ie 

When approaching the singularity from the left, det(S'4) > 0, the real part turns 
singular. When we cross the leading singularity from the right, det(S'4) < 0, the 
imaginary part of the singularity switches on, while the real part vanishes. In this 
example, both the real and imaginary parts are singular because det(S'4) changes sign 
when the leading singular point is crossed. 

At the position of the two 3-point sub-LLSs, the imaginary part shows logarithmic 
divergences while the real part has two finite negative jumps. This is similar to that 
shown in Fig. 15.21 whose explanation is given in Eq. f l5.12p . However, there is an 
important difference between Fig. 15.61 and Fig. 15.21 in the sign of the singularities. 
This is because the sign of the 3-point LLS is (—1)'^ (see Eq. (14.491) for = 3) while 
the sign of the 3-point sub-LLS is (-1)^ (see Eq. fH72|) for = 4). 



5.2.3 Conditions on external parameters to have LLS 

The conditions for the opening of normal thresholds give the lower bounds on ex- 
ternal parameters M^, s, Si^2 as given in Eqs. (15.161) and (I5.17p . However, we have 
learnt from section 14.41 that the LLS can terminate as those external parameters in- 
crease. The conditions for the termination of LLS define the upper bounds of those 
parameters, as illustrated in subsection 15.2.11 for the case of 3-point function. 

However, the situation becomes much more complicated in the case of 4-point 
function since there are 4 variables {Mh, s, si, S2) and 2 parameters {Mw, frtt) 
involved. We will show that there are two ways to find out the upper bounds by 
using numerical and analytical methods. 
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Figure 5.7: The region of leading landau singularity. 



We first explain the numerical method to find the upper bounds of Mfj and s. This 
is done by a very simple Fortran code which includes the following steps. For each Mjj, 
what is the condition on s to have a LLS in the phase space? The Landau determinant 
takes the form detlQ^) = as^ + bs2 + c as given in Eq. f l5.20p . If A = 6^ — 4ac < 
then there is no LLS. If A = 6^ — 4ac > then the Landau determinant can vanish 
at 2 points 



4 



2a 



(5.25) 



If S2 are not in the phase space defined by Eq. fl5.14p then there is no LLS. If at 
least one root is in the phase space then we have to check condition f l4.27p . If this 
condition is satisfied then there is a LLS. The result is shown in Fig. l5.7[ We conclude 
that the LLS occurs when 2Mw < Mh < 211GeV and 2mt < ^/s < 457GeV. The 
range of LLS region depends on Mw and rrit. If mt/M]^ < 1 then the first two 
conditions in Eq. f l4.26p can never be satisfied. In particular, if rrit/Mw = 1 then 
the Landau determinant can vanish but the sign condition cannot be realised. When 
Mh > 211GeV or ^/s > 457GeV the Landau determinant detlQ^) can vanish inside 



108 



Chapter 5. SM bbH production at the LHC: Mh > 2Mw 



the phase space but the sign condition > cannot be fulfilled. 
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Figure 5.8: The real part of Dq as a function of y/s^ for various values of Mh- 

4 + Ml.] 



For Mh = ^Mw we have taken si = 2{m^ + M^). For the other cases, we take 



si = 260 GeV. 



Before explaining the analytical method, we would like to show pictorially how 
the LLS moves and terminates as Mh increases. We fix ^/s = 353GeV as in Fig. 15.61 
We will increase Mh and look at the behaviour of the scalar 4-point function Dq as 
a function of S2- We will explain what value of si should be chosen. The result is 
shown in Fig. I5.8l which is just an extension of Fig. 15. 6[ The key points to understand 
this picture are as follows. At most, there are four discontinuities in the real part 
as a function of S2 as already explained (see Fig. 15. 6p . When we fix s and increase 
Mh, two of them are fixed and the other two move. The normal threshold is fixed at 
^/s^ = mt+Mw = 254.3766GeV, one 3-point sub-LLS is fixed at ^if = 259.576GeV 
as given in Eq. (15.231) . The position of the other 3-point sub-LLS depends only on 
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Mh as given in Eq. fl5.22p . The position of LLS depends on Mh and Si as given in 
Eq. (15 .25^ . As Mh increases, the Fig. 15.81 shows: 

• For Mh = 159GeV < 2Mw, only the normal threshold and the three-point 
sub-LLS ^/s^ show up. 

• For Mh = 2Mw the second three-point sub-LLS appears at \^s^ = 271.059GeV. 
One has to change Si in the range defined by Eq. (I5.14p with the condition 
Si > {irtt + My/Y make a LLS appeared. It is easy to find out that the LLS 
only occurs when ^ = ^2(m2 + M^) = 271.059GeV and the LLS position 
coincides with the position of the three-point singularity \/s^ ■ At this LLS 
point, the sign condition has the form = 0/0 for i = 1,2,3. We have the 
ordering < i/if < ^/s^ = \/sp^. 

• For Mh = 165GeV then a/s^ = 257.088GeV and we see that the LLS starts 
showing up at ^/s2 ~ 283.5GeV when ~ 260GeV (before this value there 
is no LLS) and moves to the left as ^/si increases. We have the ordering 

^/sf<^M<V4<^/W^■ 

• For Mh ~ 173GeV then \/s^ = 254.3766GeV coinciding with the normal 
threshold and we see that the LLS starts showing up at ~ 274GeV when 
y/si ~ 260GeV and moves to the left as ^/si increases. After this value of 
Mh, the y^sl^-three-point singularity disappears from the physical region and 
the LLS continues moving to the left as Mh increases. We have the ordering 

• For the special value Mh = 190.877GeV, the LLS starts showing up and coin- 
cides with the fixed three-point singularity at = -y/i| = 259.576GeV when 
y/si ~ 260GeV ~ i/if and moves off the physical region as ^/si increases. If 
Mh > 190.877GeV then the LLS disappears from the physical region. 

We conclude that for = 353GeV, the upper bound of Mh to have a LLS is 
190.877GeV which is consistent with Fig. 15.71 The above picture also leads to the 
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conclusion that the upper bound is determined when the LLS coincides with a sub- 
LLS. This important remark agrees with the explanation on the termination of LLS 
given in section 14.41 

We now are in the position to explain the equation of the upper bounds of Mh 
and ^/s, which is shown in Fig. 15.71 by using analytical method. The key points are 
as follows. The termination of LLS occurs when all LLSs coincides with a three-point 
sub-LLS. One should keep in mind that for each value of {M^, \/s) there may be a lot 
of LLSs corresponding to different values of (si, S2) which make a si,2-LLS-range. At 
the termination of LLS, this Si^2-LLS-range must become a point in order to coincide 
with a three-point sub-LLS. 

If we express det^Qi) as a quadratic polynomial of S2, there are 2 three-point 
sub-LLSs whose positions are given in Eqs. f l5.22p and f l5.23p . Without losing any 
generality, we assume that the former coincides with the LLS. Thus, at the termina- 
tion point we have 

S2 = sf and det(Q4) = 0. (5.26) 

However, since our problem is completely symmetric under exchange of Si and S2, one 
should observe the same thing when expressing detlQ^) as a quadratic polynomial of 
Si. The LLS therefore coincides with the three-point sub-LLS when 

si = sf and det(Q4) = 0. (5.27) 

We conclude that the LLS terminates when Si = S2- With this information, the Eq. 
det(Q4) = givei] 

s^ = S2 = 1[ml + Ml,) - ^(s - 4m2)(M^ - 4M2,). (5.28) 
Equating the Eq. fOHl) with Eq. (E22]), we get 

4^^-^^ = ^ (MH(m2 - M^) - (m? + M^)^(M^ - AM^,)^ . (5.29) 
^ Other roots do not satisfy the positive sign condition of the three-point sub-LLS. 
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We observe that this equation shows, in a very transparent way, that all thresholds: 

mt > Mw, Mh > 2Mw, Vs > 2mj 

need to be open simultaneously. We can invert Eq. fl5.29p to write the solution in 
terms of Mfj. To arrive at the same result, it is more judicious however to go through 
exactly the same steps but choosing instead of 82- We derive 

{Ml -AMI,) = ^ (v^(m? - Ml) - (m? + M^) ^ {s - Aml)^ . (5.30) 

The maximum value of Mh is obtained by setting y/s = 2mt {Mh = 2Mw), 

i.e. when the LLS, the two 3-point sub-LLSs and the normal threshold coincide. We 
have 



mi 



4m? < s < Am\ + ^""^ ..f^^' - (5.31) 

Numerically, it means 

348GeV < < 457.053GeV and 160.7532GeV < < 211.129GeV. (5.32) 

Of course, those analytical formulae agree with the curve obtained by numerical 
method in Fig. 15.71 

There are many other ways to derive Eq. 05.291) . Let us explain a very practical 
way. In order to obtain Eq. (15.291) . we have assumed that si,2 = S2 as in Eq. (I5.26p . 
However, we can also equally assume that S\^2 = ^2 i-e. the LLS coincides with the 
other three-point sub-LLS. Of course, this assumption does lead to the same resulc. 
But this also means that, in order to have a unique result for the equation of upper 
bounds, at the termination of LLS one must have 



(5.33) 



^This will fail if the equation of the upper bounds of s and Mfj (the maximum curve) is not 
analytic. If it happens, each assumption will give a part of the maximum curve. 
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Thus, a very practical way to quickly obtain the result Eq. fl5.29p is equating two 
equations (15.221) and (15.231) . without caring about Eq. det(Q4) = 0. 

The same argument can be repeated for the two parameters Si^2 to get their up- 
per bounds. Namely, we express det((54) as a quadratic polynomial of M^. The two 
3-point singularities are found by considering the two reduced triangle diagrams con- 
taining p5 as an external momentum. The equation of the upper bounds is obtained 
by requiring that Mfj{si) = M^(s2) (similar to Eq. fl5.33p ): 



^ S2{ml + Ml) - [ml - Ml,f - [mj - M^) ^\{s^, mf, M^) 
si(m? + M2,) - (mj - Ml,y - {m\ - M^,) ^A(si, ml ^ 



W) 



2ml 



(5.34) 



which is compatible with the sign condition < [si,2{ml+Mly)-{ml-Mly)^]/{2ml) 
(the "plus" solution does not satisfies this); the equation writes 



_ _ ml-Ml, 
' ' ml + Ml 



^A(si,m2,M2,) - ^A(s2,m?,M2,) 



(5.35) 



w 

2 



The maximum of S2 is got from this equation by setting Si = (m^ + M^y) . We get 

(m, + Mw)' < s,,2 < (m, + M^)' + ^ ' ^^ (5.36) 

mtMw 

which gives 

254.3766GeV < < 324.442GeV. (5.37) 

From Eq. 15.351 one can make a very similar plot like the one in Fig. 15.71 We remark, 
by looking at Eqs. fl5.3ip and fl5.36p . that there is no LLS if mj = Mw (because 
> Mh). 

A good question to ask is "What does the termination of the LLS mean physi- 
cally?" The answer is very simple if one uses the physical interpretation of Coleman 
and Norton discussed in section 14.31 The relation between the Feynman parameter Xi 
and the proper time dri = miXi where m^ = Mw, frtt means that at the termination 
of the LLS (when Mh or/and are large enough) at least one internal particle has 
zero proper time, i.e. it reaches the velocity of light. 
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5.3 The width as a regulator of Landau singulari- 
ties 



We will argue that the LLS which is not integrable at the level of one-loop amplitude 
squared can be tamed by introducing a width for unstable particles in the loops. 
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Figure 5.9: Effect of the width of the W , Tw and the top, Tt, on the real and imaginary 
parts of the scalar four-point function. 



It is well-known that if internal unstable particles are on their mass shell, then one, 
in general, has to introduce a width for those particles. For example in the process 
e~^e~ ff at tree level, there is a pole in a S-channel diagram associated with the 
Z boson. It is obvious that one has to include the Z-width to solve this problem. 
On the other hand, there are several processes with unstable internal particles where 
the width effect is very small hence can be safely neglected. A famous example is a 
T-channel diagram. In that case, the internal particle is forced to be far away from 
the on-shell region. 

In our case at hand, the LLS occurs when all loop internal particles (the top-quark 
and W Goldstone boson) are on-shell. Thus the width effect can be important. We 
take the simple prescription of a fixed width and make the substitution 

m^t m^t - imtVt, ^ - iMw^w (5.38) 
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Mathematically, the width effect is to move Landau singularities into the complex 
plane, so they do not occur in the physical region (the real axis). 

Applying rules fl5.38p to the case of the scalar four-point function defined by Eq. 
f l5.13p one sees in Fig. I5.9l that indeed the width regulates the LLS and gives a smooth 
result that nicely interpolates with the result at zero width away from the singularity. 
The normal threshold and the 3-point sub-leading singularity are also softened. The 
real part of the 4-point function still shows a smooth valley at the location of the 
LLS after regularisation. For the imaginary part we note that after introducing the 
width the LLS singularity is drastically reduced with a contribution of the order of 
the sub-leading singularity. 

As we will explain in the next section and in more detail in Appendix [E] the intro- 
duction of the width in a four-point function requires careful extension of the usual 
4-point function libraries. 

In our calculation of Yukawa corrections where all the relevant couplings depend 
only on the top-quark mass and the vacuum expectation value v, we will keep rrit and 
V real while applying rules (15.381) to all the loop integraljf]. 



5.4 Calculation and checks 



The one-loop calculation of this chapter is exactly the same as in the real mass case 
(see section except the fact that we now have to consider the tensorial and scalar 
loop integrals with complex masses. 

LoopTools [33] can handle the complex masses up to 3 point functions. The 5 point 
functions are reduced to 4 point functions [I2l [3l] . The tensorial 4 point functions 



•^In a general calculation of full electroweak correction and if one chooses the input parameters 
to be {a{0), Mw, Mz, mt, . . .} then rules (|5.38p should be applied everywhere in the calculation so 
that the EW gauge invariance can be preserved. This idea of doing analytical continuation on gauge 
boson masses is the philosophy of the complex-mass scheme \117\ [T5] . For a practical discussion of 
methods to deal with unstable particles, see Denner's lecture notes at [118] . 
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are reduced to the scalar 4 point function and 3 point functions. We therefore have 
to calculate only the scalar 4 point function. The analytical calculation of scalar 4- 
point function with complex masses in the most general case is practically intractable. 
The standard technique of 't Hooft and Veltman [119] has some restriction on the 
values of external momenta. In particular, the method works if at least one of the 
external momenta is lightlike. In that case, the result is written in terms of 72 Spence 
functions. In our present calculation, there are at least 2 lightlike external momenta 
in all boxes. If the positions of two lightlike momenta are opposite then we can write 
the result in terms of 32 Spence functions. If the two lightlike momenta are adjacent, 
the result contains 60 Spence functions. The detailed derivation and results are given 

le case of 



All the 



in the Appendix [El We have implemented those analytical formulae for t 
two massless external momenta into a code and added this into Loop Tools 
five point functions which have problem with Landau singularities in our calculation 
have 4 external massless momenta (two gluon and two bottom quarks). Thus they 
can always be reduced to a set of 4- point functions with at least 2 massless external 
momenta. 

Checks on the loop integrals with complex internal masses: for all the tensorial 
and scalar loop integrals (4- and 5- point functions), we have performed a consistency 
check where we have made sure that the results with complex internal masses become 
asymptotic to the ones with real internal masses in the limit widths 0^. For the 
scalar loop integrals, the results are compared to the ones calculated numerically 
in the limit of large widths, e.g. w = lOOGeV. Furthermore, for the scalar box 
integrals the results can be checked by using the segmentation technique described 
in [120] . The idea is the following. At the boundary of the phase space, the 4- point 
functions can be written as a sum of 4 three point functions. The 3-point functions 
with complex masses can be calculated by using LoopTools. In this way, we have 
verified that the results of the scalar 4-point functions are correct at the boundary of 



'*The implementation for the case of one massles external momentum is straight forward. How- 
ever, we have not done this yet since it is not necessary for our present calculation. 
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phase space. We have also carried out a comparison with a purely numerical approach 
based on an extension of the e-extrapolation technique [121] . We have found perfect 
agreement^. 

For the whole calculation, we have performed two checks at the amplitude level. 
First, by taking the limit widths 0"^ we have observed that the results approach 
to the results calculated with real internal masses. This again is just a consistency 
check. Second and it is the most important check where we have verified that the 
results calculated with complex internal masses are QCD gauge invariant. 

As the LLS is integrable at interference level, the NLO calculation with Xi,f,H 7^ 
performed in chapter [3] can be trivially extended to the region of Mh > 2Mw by using 
the same method without introducing widths for unstable internal particles. However, 
there is a small problem related to the universal correction {SzU"^ — 6v) where SZ]^"^ 
related to the derivative of the Higgs two-point function becomes singular when 
equal to 2M]y or 2Mz- This problem is separately treated by introducing the widths 
of the W and the Z. To be complete, the results are given in section [531 



5.5 Results in the limit of vanishing A^^i^ 

The input parameters and kinematical cuts are the same as given in section 13. 5[ We 
write here additional parameters related to the widths of unstable particles appearing 
in the calculation: 

Tw = 2.1GeV, Tz = 2.4952GeV, (5.39) 
the top-quark width is calculated at the tree level in the SM 



G^K- M^)^K + 2My^ 



r, = ^^^wJV'^y-^^^wJ ^ 1 5GeV (5.40) 



where the bottom-quark mass has been neglected. 



^We thank F. Yuasa for sending us the resuhs of the extrapolation technique. 
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5.5.1 Total cross section 
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Figure 5.10: Left: the cross section for the subprocess gg — > bbH as functions of Mh 
for various values of y/s including the case a/s = 2mt = 348 GeF. Right: the cross 
section for the subprocess gg — > bbH as functions of ^/s for various values of Mh 
including the case Mh = 2Mw = 160.7532 GeF. 

We start with the cross section in the case where XbbH = 0. In section 13.51 we 
reported on resuhs up to Mh = 150GeV that showed that this cross section was 
rising fast as one approached the threshold Mh = 2Mw- Beyond this threshold our 
integrated cross sections showed large instabilities. As discussed in subsection 15.2.21 
this is due to the appearance of a leading singularity which as we have advocated 
can be cured by the introduction of a width for the unstable top-quark and W gauge 
boson. Before convoluting with the gluon distribution function let us briefly look at 
the behaviour of the partonic cross section gg — > bbH paying a particular attention 
to this leading Landau singularity region shown in Fig. I5.7[ see also Eq. I5.31[ 

Figs 15.101 show that indeed the widths do regulate the cross section. Moreover it 
is within the range of LLS that the cross section is largest. The (highest) peak of the 
cross section occurs for a Higgs mass about 163GeV (slightly higher than the normal 
threshold value 2Mw = 160.7532GeV) and for ^/s ~ 351 GeV (slightly higher than the 
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normal threshold value 2mt = 348GeV). The peak positions are slightly shifted from 
the normal threshold values, this is due to the width effect. The reason for the peak 
to occur at the normal threshold position is that all Landau singularities (leading and 
sub-leading) start showing up at this point and the size of Landau singularity region 
is largest at the position of normal threshold, see Fig. 15. 7[ To see the LLS effect, we 
look at the two curves = 347GeV and ^ = 2mt = 348GeV of Fig. ICTir ieft). 
For the former, there is no LLS and the peak effect is due to sub-LLSs (the normal 
thresholds and 3-point sub-LLSs). The later includes additional LLS contribution, 
which is significant even after being regulated. We emphasize that although the 4- 
point LLS is special in the sense that it is not integrable at the level of one-loop 
amplitude squared one should not overlook this point and use it as an account for 
the bulk of the large correction around the normal threshold position. The 3-point 
(sub-) LLSs and normal thresholds have also significant contributions. 
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Figure 5.11: Left: the one-loop induced cross section as a function of Mh in the limit 
of vanishing bottom-Higgs Yukawa coupling for two cases: with and without widths. 
Right: the percentage correction of the contribution with widths relative to the tree 
level cross section calculated with XbbH 7^ 0. 



The cross section at the pp level for the 14TeV centre of mass energy (LHC) as 
a function of the Higgs mass is shown in Fig. 15.111 taking into account the width 



5.5. Results in the limit of vanishing Xb^H 



119 



of the top-quark and the W gauge boson. For comparison we also show the cross 
section without the width effect outside the leading Landau singularity range of Mh- 
The sharp rise above > ISOGeV is nicely tamed. On the other hand note that 
on leaving the leading Landau singularity region around Mh = 211GeV, the width 
effect is much smaller and the figures suggest that one could have entered this region 
from the right without having recourse to introducing a width. Indeed our numerical 
integration routine over phase space with the default Loop Tools library does not show 
any bad behaviour until we venture around values of 2Mw < Mh < 200GeV. The 
reason for this can be understood by taking a glance at Fig. 15.71 For 200GeV < Mh < 
211GeV the singularity region is considerably shrunk so that one is integrating over 
an almost zero measure. The effect of the widths outside the singularity region is to 
reduce the cross section for Mh = 120GeV, 140GeV and ISOGeV bjo respectively 
15%, 24% and 33% while for Mh = 210GeV, 230GeV and 250GeV the reduction is 
comparatively more modest with respectively 15%, 5% and 2%. 

The relative correction to the tree level cross section is shown in Fig. 15.111 (right). 
It amounts to 2.6% for Mh = 120GeV, increases to as much as 49% when Mh = 
163GeV and finally becomes almost constant at about 10% for large values of Mh- 
Large contribution is due to the effect of Landau singularities. 



5.5.2 Distributions 



^The relative difference is defined by [cr(r = 0) - cr(r ^ 0)]/cr(T ^ 0). 
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Figure 5.12: Left: The pseudo-rapidity of the Higgs and transverse momentum dis- 
tributions of the Higgs and the bottom for Mh = 150GeV arising from the purely 
one-loop contribution in the limit of vanishing LO (XbbH = ^) for two cases: with and 
without widths. Right: Its relative percentage contribution da{XbbH = ^)/daLo is also 
shown. 
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[do(r=0)-do(r^O)]/da(r*0)[%] [do(r=0)-da(r^O)]/dCT(r7iO)[%] 
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Figure 5.13: The relative difference between two cases: with and without widths, 
defined by [da(T = 0) — da(T ^ 0)]/(i(T(r 7^ 0)[%], for the transverse momentum 
distributions of the Higgs and the bottom quark. 

In order to see the width effect on distributions, we first consider the case where 
Mh = ISOGeV < 2Mw (no LIS in this case). Figures 02) and 15131 show the 
difference between the two cases: with and without widths. The relative difference 
is rather uniform, about 33%, on the Higgs pseudorapidity distribution. For the 
transverse momentum distributions, the relative difference is not uniform but has 
structure as shown in Figure I5.13[ We remark that the peaks in the transverse 
momentum distributions occur at the position where the width effect is largest, hence 
are related to the opening of sub-leading Landau singularities as discussed in the 
previous subsection. 

The largest relative corrections to the tree level distributions are shown in Fig. 
15. 141 for the case Mh = 163GeV > 2Mw (LLS does occur here) which corresponds to 
the maximum value of the cross section in the limit XbbH = as displayed in Fig. 15.111 
(right). The correction to the Higgs pseudorapidity distribution is about 60% around 
the center region. The corrections to the px distributions can be enormous in some 
region of phase space, up to 200% for the Higgs and about 170% for the bottom-quark 
case. These huge corrections to the distributions in some region of phase space are 
again due to the effect of Landau singularities. 
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Figure 5.14: Left: The pseudo-rapidity of the Higgs and transverse momentum dis- 
tributions of the Higgs and the bottom for Mh = IQSGeV arising from the purely 
one-loop contribution in the limit of vanishing LO (XtbH = 0^. Right: Its relative 
percentage contribution da{XbbH = ^)/daio is also shown 
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5.6 Results at NLO with X^^h ^ 

The purpose of this section is to complete the study of subsection 13.5.21 to cover 
higher values of Higgs mass. Moreover, we would like to know the width effect at 
NLO in the presence of Landau singularities in various one-loop diagrams. We recall 
that the LLS is integrable at NLO. 

As discussed in subsection 13.2.31 the NLO Yukawa corrections consist of 3 gauge 
invariant classes, see Fig. 13.41 The study of subsection 13.5.21 revealed that class (a) 
gives a totally negligible correction below 0.1%. We will not discuss this contribution 
any further here. Moreover, leading Landau singularities we have discussed only show 
up in class (c). We will therefore study separately the NLO correction due to class 
(c) and weigh the effect of implementing the width of the internal particles. Class (b) 
does not develop any leading Landau singularity and therefore the width effect will 
be marginal. 

5.6.1 Width effect at NLO 

The class (c) has problem with the 4-point LLS. However, this LLS is integrable 
at NLO. Thus one can calculate the NLO cross section without introducing widths 
for unstable particles in loops. However, one can still expect some effect of Landau 
singularities (leading and sub-leading) and wonder if the width effect is significant in 
this case? The answer to this question is given in Fig. I5.15[ 

The results in this plot are calculated by setting ruf, = in the kinematics (spinors 
and propagators) and loop integrals while keeping = 4.62GeV only in the XbbH 
coupling which can be regarded as an independent parameter. The first remark is 
that if Mh < 158GeV or Mh > 165GeV then the NLO width effect is below 5%. For 
Mh < 158GeV the W-Goldstone bosons in the loops can never be on-shell and thus 
the width effect is completely small. For Mh > 165GeV the W-Goldstone bosons can 
be on-shell and thus the width effect is a bit bigger than in the former case. Indeed 
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Figure 5.15: The relative difference between two cases: with and without widths, 
defined by [da(T = 0) — da(T ^ 0)]/da{T ^ 0)[%], for the NLO correction of the class 
(c) where the LLS occurs. The tree-level amplitude is calculated with massive bottom 
quark. The one-loop amplitude is calculated by keeping mi, only in the XhbH coupling. 



even at the peak Mh = 2Mw the width effect is just about 17%. From this analysis 
we conclude that the width effect at the NLO is small and can be neglected. 

5.6.2 NLO corrections with mi, ^ 

The results for the NLO corrections are shown in Fig. 15.161 (left). For classes (b) 
and (c) the widths of unstable particles are neglected. For the universal correction, 
{SZ^^'^ — 6v) where SZ^^'^ involving the derivative of the two-point function Higgs 
self-energy diverges when Mh is equal to 2Mw or 2Mz, all the widths of unstable 
particles are kept B 

The effect of Landau singularities is obvious if one compares the (c)-curve to the 
(b)-curve. The contribution from class (b) where the Higgs couples to the internal 
top decreases very slowly as the Higgs mass increases from llOGeV to 250GeV, as 
expected there is no structure as would be the case if this contribution were sensitive 

^Note that SZ^^ does not diverge when Mh = ^rrit and the top-quark width thus has a marginal 
effect. 
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to any threshold or singularity. Class (c) on the other hand does, as expected, reveal 
some structure around Mh = 2Mw where we see a fall in the relative correction. 
The correction is however, despite this fall, quite modest ranging from about — 1% 
for Mh = leOGeV to about -4% for Mh = 210GeV. 

The universal correction, {SzU"^ — 6v), contributes from —1% to —3%, where the 
highest correction is at the H — >• WW threshold {Mh = 2Mw)- 

The detailed structure of the class (c) is shown in Fig. 15.161 (right). It consists of 
two independent helicity structures where the helicities of two bottom quarks in the 
final state are the same or different. We call them S\^^\^ and 5a3,-A4 structures. In the 
massless bottom limit whose result is displayed in Fig. I5.15[ only the 5a3,-A4 structure 
survives. It is remarked that the behaviors of two different helicity structures as 
functions of Mh are very different despite the fact that they both have a common 
denominator. 




Figure 5.16: Left: The relative NLO EW corrections normalized to the tree-level cross 
section, (h) and (c) correspond to the two classes of diagrams displayed in Fig. \3.4 

1/2 ~ 

{SZjj — 6v) is the universal correction contained in the renormalization of the bbH 
vertex. "Total" refers to the sum of those 3 corrections. SzU"^ is calculated by taking 
into account the widths of W , Z and the top quark. Right: the structure of (c)- 
correction which is a sum of two independent helicity configurations. 
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5.7 Summary 

We have extended the study of section [331 to cover higher values of Mh- This study 
is important and highly nontrivial due to the appearance of Landau singularities in 
various one-loop Feynman diagrams when Mh > 2Mw- 

We have applied the general study of Landau singularities, presented in chapter 
m to the specific case of gg bbH. Namely, we have studied in detail the LLS in the 
case of one-loop 3-point and 4-point functions. The nature of those LLS is carefully 
explained by using two different methods. On one hand, we used the general formulae 
(which are for the singular part only) obtained in section 14. 4[ On the other hand, 
we used explicit results obtained by performing loop integrals using the traditional 
method of 't Hooft and Veltman. The latter results in various plots of real and 
imaginary parts of scalar loop integrals. This indeed helped us to understand much 
better the various structures of Landau singularities. 

We have performed a detailed study to understand how the LLS terminates. From 
this, we got the upper bounds of external parameters. 

We have argued that the problem of 4-point LLS, which is not integrable at the 
level of one-loop amplitude squared, is solved by introducing a width for internal 
unstable particles. In order to do so, we have applied the loop calculation method of 
't Hooft and Veltman to write down explicitly two formulae to calculate scalar box 
integrals with complex internal masses (see Appendix [E]). The restriction is that at 
least two external momenta are lightlike. We have implemented those two formulae 
into the library Loop Tools. 

We have studied the width effect in the presence of Landau singularities in various 
one-loop diagrams. At NLO, the width effect is negligible. At NNLO, the width effect 
is extremely important around the normal threshold position {Mh = 2Mw)- 

We have shown various results of one-loop corrections to the cross section as well 
as important distributions. NNLO corrections, calculated in the limit XbbH = 0, 
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can be very large in some region of phase space or parameter space where Landau 
singularities show up. For the pr distribution, those corrections can be enormous. 
NLO corrections remain small although some structure related to Landau singularities 
does show up. 



Chapter 6 



Conclusions 



In the previous chapters, we have explained how to calculate one-loop Yukawa cor- 
rections to the process of SM Higgs production associated with two high pt bottom 
quarks at the LHC and their physical content. In particular, we have studied one-loop 
Landau singularities and shown how to handle them in practice for our process. The 
properties and effect of the Landau singularities in the process pp — > bbH were also 
carefully investigated. 

The entirely dominant contribution to the process pp — > bbH is the sub-process 
gg — > bbH. The physics of one-loop Yukawa corrections to this sub-process is very 
rich due to the fact that the bottom-quark and the top-quark are in the same SU (2) l 
doublet with a large mass splitting. At tree level, the Higgs boson couples directly 
to bottom quarks with a small Yukawa coupling X^bH- At one- loop level, the Higgs 
can couple to heavy particles hke the top-quark or W Goldstone bosons. Thus, it 
is easy to see that, if we do the expansion in the small coupling XbbH-, the one-loop 
corrections can still give some contribution when this coupling vanishes. The salient 
property of the one-loop corrections is related to the smallness of and is best 
seen if we take the limit m;, — > while keeping XbbH (which can be regarded as a 
different parameter) unchanged. In this limit, the tree level contribution contains 
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only even helicity configurations ("even" means that both bottom-quarks have the 
same helicity) while the one-loop contribution includes both even and odd helicity 
configurations ("odd" means that the two bottom-quarks have opposite helicities). 
The even one-loop contribution is proportional to the small XbtH while the odd one, 
which comes from insertion in loops, is independent of XbbH and proportional to 
large couplings like Xuh and Xhhh- We remark immediately that the NLO correction 
is proportional to A^^^^^ hence must be small. However, the most interesting thing is 
that the purely one-loop contribution which can be extracted by setting XbbH = 
consists of only odd helicity configurations. Thus, this new loop induced contribution 
can be very different from the tree level or NLO contributions. Moreover, the best 
way to understand Landau singularities, which are intrinsic analytic properties of 
loop integral, is to look at purely loop corrections. 

The numerical results of chapters [3] and [5] confirm those qualitative conclusions. 
The NLO correction is small and changes from —4% for = 120GeV to —8% at 
Mh = 2Mz stabilising to around —7% for larger values of Mh up to 250GeV, despite 
the appearance of various Landau singularities. The purely one-loop contribution 
(NNLO), calculated by taking into account the widths of the top-quark and W gauge 
boson, amounts to 2.6% for Mh = 120GeV, increases to as much as 49% when 
Mh = 163GeV and finally becomes almost constant at about 10% for large values of 
Mh- The difference between the two corrections becomes clearest when one looks at 
the pt distributions. The NNLO correction to pt distributions can be enormous in 
some region of phase space for Mh about 2M]v while the NLO correction is modest. 
Large NNLO corrections are due to the effect of Landau singularities and occur around 
the normal threshold position Mh = 2M\y. 

Calculating one-loop corrections contains a lot of technical features. The ampli- 
tude squared was calculated by using the helicity amplitude method (HAM) where 
each helicity amplitude, which is a complex number, is numerically calculated. In or- 
der to do so, we have reduced all loop integrals (up to five-point function) in terms of 
Passarino-Veltman loop functions. This was easily done with the help of Loop Tools. 
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After doing so, we obtained a huge algebraic expression for each hehcity configuration 
which takes a lot of computer time to calculate. The problem was how to optimise 
the calculation. We found a systematic way to do this by using FORM based on the 
fact that each helicity amplitude can be factorised in terms of independent blocks. 
The advantage is that some comphcated blocks which appear several times in the 
computation can be put in a common sector hence just need to be calculated once. 

In our calculation with two gluons as external particles, the HAM leads to a very 
easy and convenient way to check the QCD gauge invariance which means that the 
amplitude squared is independent of the reference vectors needed to define the gluon 
polarization vectors. Indeed, this is a very powerful way to check the correctness of the 
results and can be used for other processes with at least one external gluon/photon. 

Another advantage of HAM is that since the tree-level and one-loop helicity am- 
plitudes are calculated separately, the one-loop amplitude squared is immediately 
obtained when the NLO calculation is done. 

There is a special class of one-loop gg bbH diagrams where the Higgs boson 
is produced by W gauge boson fusion. If the Higgs mass is heavy enough for this 
normal threshold to open {Mh > 2Mw), then all Landau singularities of two-, three-, 
four- point functions show up. In particular, the four-point leading Landau singularity 
leads to severe numerical instabilities when we calculate the cross section involving the 
square of a one-loop amplitude. We have solved this problem by taking into account 
the fact that W gauge boson and the top-quark are unstable particles hence have a 
width. For this, we have followed the standard technique of 't Hooft and Veltman to 
calculate scalar four-point functions with complex internal masses. The restriction is 
that at least two external momenta are lightlike. We have implemented those formulae 
into the library Loop Tools. Various checks have been performed to make sure that this 
implementation is correct. We have also observed that the same implementation can 
be done for the case of one lightlike external momentum. However, the calculation of 
the scalar four-point functions with complex internal masses in the most general case 



132 



Chapter 6. Conclusions 



with no restriction on the external momenta is not tractable if one uses the method of 
't Hooft and Veltman. Another disadvantage of this method is that the results, even 
in some special cases with massless external particle, contain many Spence functions. 
It may be better to find another way. 

Although the main calculation of this thesis is for a very specific process, it is quite 
obvious for us that some of the results discussed above can be used or generalised for 
other complicated one-loop calculations. 



Appendix A 



The helicity amplitude method 



A.l The method 



We use a combination of helicity amplitude methods as described in [361 I122 l| to 
calculate the total cross section. In the following we only want to highlight some 
key features that were most useful for our calculation, for details of the method we 
refer to[M 1122] . For our process g{pi, Ai) + g{p2, X2) b{P3, A3) + b{p4, A4) + H{p^) 
where the particles are denoted by their momentum pi and helicity Aj we write the 
corresponding helicity amplitude as ^(Ai, A2; A3, A4). 



T^"^ is a string of Dirac 7 matrices. These 7 matrices represent either interaction ver- 
tices or momenta from the fermion propagators. In our case the interaction vertices 
are the vectorial gluon vertices in which case they represent , the scalar Higgs ver- 
tex and at one-loop the pseudo-scalar Goldstone coupling. For the momenta, in our 
implementation we re-express them in terms of the independent external momenta 
Pi,P2,Ps,P4- This applies also to the loop momenta after the reduction formalism 



v4(Ai, A2; A3, A4) 
A<^"(A3,A4) 



^^^{Pl, Ai; gi)e^(p2, A2; q2)M'"'{X3, A4), 

u{ps,X3)r>'MPi,x,). 



(A.l) 
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of the tensor integrals has been performed. The first step in the idea of the hehcity 
formahsm we follow is to turn each of these 7 matrices (apart from the pseudo-scalar 
and the trivial scalar) into a combination of spinor function uu. We therefore trans- 
form our helicity amplitude into products of spinors such as the helicity amplitude 
could be written like a product u uu ...uu v with the possible insertion of 75 's in the 
string. The different u, u, v in the string we have written have of course, in general, 
different arguments. Nonetheless one can turn each spinor product of two adjacent 
uu, etc into a complex number written in terms of the momenta in our problem as 
we will see. 

In the first step, for the momentum with pf = mf we use 

jli = u{pi, -)u{pi, -) + u{pi, +)u{pi, +) - rrii. (A. 2) 

The transverse polarization vector of the initial gluon i, e^(pj, A^; g^), is also first 
expressed in terms of spinors such as 

e,{p^,\;qi) - [4(p._^^)]i/2 ' l^-^i 
where qt is an arbitrary reference vector satisfying the following conditions 

^•=0, P^.q^y^O. (A.4) 

Using the trace technique one can easily prove that definition (IA.3P indeed satisfies 
the physical polarisation sum identity given in Eq. (11. 4p . namely 

e^(p. A; q)e^{p, A; q)* = -g^^, + ^^^'^ ^''^^ . (A.5) 

Gauge invariance (transversality condition) requires that the cross sections are in- 
dependent of the choice of the reference vector as we will see later. This acts as an 
important check of the calculation, see later. It is not difficult to prove that the choice 
(1A.3P satisfies all the conditions for a transverse polarization vector. In particular, 

Pi.e{pi,\i) = 0, e{pi,Xi).e{pi,Xi) = 0, 

ef,{pi,-Xi) = ef,{pi,\iy, e{pi,Xi).e{pi,-Xi) = -1, (A.6) 
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where the reference vector is not written down exphcitly. i = 1 , 2 and no sum over i 
must be understood. Then for e^7^ one uses the so-called Chisholm identity 

u{p, Xhf,u{q, X)Y = 2[u{p, -X)u{q, -A) + u{q, X)u{p, A)], (A.7) 

where all the spinors in Eq. flA.7p are for massless states in view of the lightlike 
condition on the reference frame vector and of course the momentum of the real 
gluon. 

With U{pi, Xi) representing either u{pi, Aj) or v{pi, Aj) one uses the general formu- 
lae 

Ax,x^{pi,Pj) - MiBx,x,iPi,Pj) + MjCxiXj{Pi,Pj) 



U{pi,Xi)lbU{pj,Xj) = -Xi 



y/{pi.h){Pj-ko) 



where 

Mi = +mi if U{pi,Xi) = u{pi,Xi), 

Mi = -rrii if U{pi,Xi) = v{pi,Xi), 

Ax,x, = Sx,-x,Xi {{ko.pi){ki.pj) - {ko.pj){ki.pi) - iXie^upak'^klpIp'^) , 

Bx^x, = 5x,x,{h.pj), Cx^x, = ^Kx.ih.p^), (A.9) 

with /cq,! being auxiliary vectors such that k^ = 0, kf = —1 and kQ.ki = 0. No sum 
over repeated indices must be understood. For instance, we can choose ko = (1, 0, 1, 0) 
and ki = (0,1,0,0). With this choice, it is obvious to see that the denominator in 
Eq. flA.Sp can never vanish if the bottom mass is kept. If one would like to neglect 
rrib, that choice can bring p^.ko or p^.ko to zero in some cases. If this happens, one 
can tell the code to choose ko = (1, 0, —1, 0) instead of the above choice. In fact, that 
is what we did in our codes. 

In the case of spinors representing a massless state, the helicity formalism simplifies 
considerably. Only Ax^x is needed. Traditionally we introduce the C-numbers s{p, q) 
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and t{p, q), 

s{p,q) = u{p,+)u{q,-) = A+_{p,q), t{p, q) = u{p, -)u{q, +) = -s{p, q)* , 
k(pi,P2)P = s{pi,p2)t{p2,Pi) = 2pi.p2. (A.IO) 

These are the functions that appear in our code for the massless b quark. The massless 
case is also used when expressing the gluon polarisation vector to which we now turn. 



A. 2 Trans versality and gauge invariance 

The reference vector used for the polarisation of the gluon can be changed at will. 
Changing the reference vector from q to q' amounts essentially to a gauge transfor- 
mation. Indeed one has 



A; q') = e'^^'^' ''^^e^{p, A; q) + g^, (A-ll) 



where 



Qi<f>(<i' ,q) 



s{p,q) t{p, q') 
t{p,q) s{p,q') 



1/2 



Therefore up to the phase factor, the difference is contained in the momentum vector 
of the gluon. QCD gauge invariance for our process leads to the important identity 

|^(Ai, A2; As, A4;gi,g2)P = |^(Ai, A2; A3, A4; g'l, ga)^. (A.13) 

as long as q'12 satisfy condition (]A.4I) . We have carefully checked that the numerical 
result for the norm of each helicity amplitude at various point in phase space is 
independent of the reference vectors gi,2 up to 12 digits using double precision. By 
default, our numerical evaluation is based on the use of gi_2 = {p2,Pi)- For the checks 
in the case of massive b quarks the result with qi^2 = {P2,Pi) is compared with the one 
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using any gi 2 such as conditions (lA.4p are obeyed. In the case of massless b quarks 
it is simplest to take gi,2 = (^35^4)- 

This check is a an important check on many ingredients that enter the calculation: 
the Dirac spinors, the gluon polarization vectors, the propagators, the Lorentz indices, 
the tensorial loop integrals. It has been used extensively in our numerical calculation. 



Appendix B 
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B.l Optimization 

Each helicity amplitude w4(Ai, A2; A3A4) = ^(A), a C-number, is calculated numeri- 
cally in the Fortran code. The price to pay is that the number of helicity amplitudes 
to be calculated can be large, 16 in our case for the electroweak loop part. Some 
optimisation is necessary. The categorisation of the full set of diagrams into three 
gauge invariant classes as shown in section 13.2.31 is a first step. We have sought to 
write each diagram as a compact product of blocks and structures containing different 
properties of the amplitude. We write the amplitude according to a colour ordering 
pattern that defines three channels. The ordering is in a one-to-one correspondance 
with the three channels or diagrams shown in Fig. 13. 1[ The T-type is the direct 
channel, the ?7-type is the crossed one obtained from the T-type by interchanging the 
two gluons and the S'-type is the one involving the triple gluon vertex. The helicity 
amplitude for each diagram can thus be represented a^l 

j[^Xf,u,s ^ CME{a, b)xCcx ^[FFE x SME{X,)], (B.l) 

where 



^The method we use here is very similar to the one described by Denner in [51] . 
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• CME{a, b) is the colour matrix element, a, b are the colour indices of the 
two initial gluonj^. The colour products can be {T"'T^), {T^T"") or [T", T^] 
corresponding to the three T, f/, S channels respectively 

• Cc contains all the common coefficients like the strong coupling constant Qs 
or factors common to all diagrams and amplitudes such as the normalisation 
factor entering the representation of the polarisation vector of the gluon, see 
Eq. ( 1A.3I) . Cc is the same for all diagrams and is included at the very end of 
the numerical evaluation stage. 

• FFE, form factor element, contains all the denominators of propagators, loop 
functions as well as various scalar products of external momenta {pi,P2,P3,P4} 
i.e. all the scalar objects which do not depend on the helicity Aj 

• SME{X), standard matrix element, is a product of the scalar spinor functions 
Ax^Xj, Bx^Xj and Cx.Xj defined in Appendix |X1 

For each channel, say ^(A)"^, the most complicated and time-consuming part is 
the FFE. That is why we want to factorise it out and put it in a common block 
so that in order to calculate all the 16 helicity configurations of .4(A)"'" we just need 
to calculate FFE once. This is done at every point in phase space. This kind of 
factorisation can be easily carried out in FORM (see section |B.2|) . 

SME{X) is also complicated because the bottom quark is massive and 75 occurs 
in the "helicity strings". Thus we have to optimize this part as well. The way we 
do it for all the 3 groups is as follows. In FORM, we have to find out all the generic 
expressions of SME{X). There are 12 of them at tree level and 68 at one-loop if we 
choose gi^2 = P2,i for the reference vectors. For instance, 



SME 



1 



[u{\'i,p^)v{\i,pi)] X \e^{\i,Pi,P2)Pi] X [e^(A2,P2,Pi)P4]> 
BMEi{M,\i) X BME2{\i) X BME^{\2), 



(B.2) 



^ Other colour indices of the bottom quarks are omitted here for simphcity 
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can be expressed in terms of 3 basic matrix elements [BME). Each BME occurs 
several times when calculating all the SME[X). The number of BME is 31. Each 
BME is written in terms of scalar spinor functions Ax^Xj, Bx^Xj, Cx^Xj- All the SME 
or BME can be found and abbreviated in FORM. As an alternative, we can use Perl 
for such an operation. The FORM output is converted directly into a Fortran code 
for numerical evaluation. Needless to say, all the abbreviations of SME or BME 
must be put in common blocks. 

To get the final result, we have to sum over all the channels. The grouping can 
be re-arranged in terms of an Abelian part and a non-Abelian part according to 

^(A) = A{Xf + A{Xf + A{Xf, 

= {T", r^}^(A)^^^' + [T", ^^]^(A)^^'"^^ (B.3) 



where 

^(A)^"^' = ^(^(A)^ + ^(A)^), 
^^XfAbei ^ J^^xf ^^{A{Xf - A{Xf) (B.4) 

corresponding to the Abelian and Non-Abelian parts respectively. Amplitudes with- 
out color factor are denoted by a tilde. The amplitude squared then contains no 
interference term between the Abelian and Non-Abelian parts: 

I ^(A) |2= (^^ I ^(A)^^^' |2 +12 I ^(A)^^^^' A (B.5) 



256 V 3 

where 2i6~i^l^l^^ spin- and colour- averaging factor. 

The task of our FORM code is to calculate J^r-diagrams A>^f, T.U -diagrams ^(A)^ 

and TliS-diagrams -^(A)*^ ^s fuuctious of EFEs and SMEs. These algebraic expressions 
are fed into a Fortran code which uses Eqs. (lB.4p and (1B.5P to calculate the total 
amplitude squared. 
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B.2 Technical details 



What we actually do in our FORM code to find all FFE and SME{\) is the follow- 
ing. Consider, for example, a pentagon diagram in the class (b) of Fig. 13. 4[ In the 
FORM code, we have to write a subroutine to calculate the loop integral given by 



(g2 - ml) [{q + /ci)^ - mj] [{q + - ml] [{q + k^Y - m|] [{q + k^Y - ml] ' 
= (l + A75)(mi + ^ + ^i)(m2 + ^ + ^2) 

X 7'^(^3 + ^ + ^3)7'^(m4 + ^ + ^4)(l-^75), (B.6) 

where A is just some coupling constant; ki = p^, k2 = P5 + Ps, k^ = k2 — pi, 
^4 = /^s — P2, mQ = Mi4/, mi = m2 = m^ = m^ = m^. In order to write E^^" in 
terms of Passarino-Veltman functions Eijk..., Dij,,, and Cj... we have to expand the 
numerator N'^'^. We then use Dirac algebra to re-organise N^'^ such that each term 
has the form 75(. . .y x or 75(. . .)^'^ x ^. For terms with g^, or q.ki, we use 

q-ki = ^ {[{q + k^f - ml] - {q^ - ml) + {ml - ml - kl)} , 

= [{q + k^)' - ml] + {ml - kl) - 2k'lq,. (B.7) 

D and C functions appear when terms in the right-hand-side cancel with denominator 
factors. The advantage of using Eq. flB.7p is that the rank of tensorial functions is 
reduced. In the present example, Eq. ( ]B.7p helps us avoid tensorial five-point functions 
with rank 4, which are very complicated. Obviously, numerical evaluation is much 
faster this way. However, there may be a problem when using Eq. (IB. 71) if the library 
of scalar loop integrals is not complete. This is the situation of section 15.51 In that 
calculation, we have to deal with complex internal masses and our library for scalar 
four-point functions includes only special cases with 2 lightlike external momenta. 
If we use Eq. (IB. 71) . it will create four-point functions with 3 or 4 massive external 
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particles since the momenta in the denominator of Eq. (IB. 60 are shifted: ki ki — ki, 
i = 2,3,4 and {ki — fci)^ are not necessarily zero. We then have a problem with these 
four-point functions. In this situation, one should not use Eq. flB.7p but simply write 
q'^ = g-qpfffl^) q-ki = k^qri instead. The first aim is to write E!^'^ as the following 
generic expression 

Er{k,m,) = J2frh^,h)Ts{klm^), (B.8) 

s 

where Tg are Passarino-Veltman loop functions, /^''(75, ki) are complicated functions 
of g'^'^, 7^, fcf and 75. 

In order to write helicity amplitudes in simple forms, we have to choose a good 
momentum basis. In the calculation of pp — > bbH our basis is {pi,P2,P3,P4}- We 
have to replace ki by pi to get E^'^{pi,mi). One can then use the on-shell condition 
to simplify i?^*^. In practice, one knows that each Ts{pj,mf) appears several times 
in the calculation and has a very lengthy expression. For optimization and having 
compact expression of E!f^{pi,mi) we introduce abbreviations for Ts{kf,mf). Those 
abbreviations serve as the library of Passarino-Veltman loop functions for the present 
calculation. The way to do this in FORM is as follows. 

local F = Y^Ei^^{pi,mi)] bracket T^; .sort 

collect cfl; .sort 

poly fun cfl; .sort 

poly fun; id cfl(x?)=l; .sort 

print +s F; .end 

The result of this simple FORM script is that F is a sum of independent Passarino- 
Veltman loop functions Ts{k'f,mf). One can use the command "#write" to produce 
an output file if one wishes. Now one can use Perl p]^3j to introduce an abbreviation 
for each term: Ts{kf,mf) fTg. We have 

= $^/r(75,P.) X /T,. (B.9) 
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Apart from CME and Cc, a generic T-channel helicity amplitude is calculated in 
FORM as 

^(A)^ = K(A3,P3)E^r)eM(Pi,?i)e.(P2,g2)^(A4,P4). (B.IO) 

r 

We can now simplify this expression by using the transversality condition for the 
gluon polarization vectors and Dirac equation. For the latter, we have to re-organise 
fr(l5,Pi) to have the form: |/375|^i7''7V2|^4. 

U{^3,P3)^3 = "^;,fZ(A3,P3), ^4v{X4, Pa) ^ -mbV{X4, P4) . (B.ll) 

Notice that if one chooses gi = p2 and q2 — Pi then ^(A) is further simplified by 
using 

P^e^{pi:P2)^0: pre.(p2,Pi) = 0. (B.12) 

We are now in the position to factorizc each term of A{\)^ as a product of FFE and 
SME{\i). The trick to find all SME{Xi) is the same as above: using the combination 
(bracket, collect, polyfun). 

local FT = ^(A)^; local FU = ^(A)^; local FS = ^(A)^; 
bracket fZ, v, e, 7; .sort 
collect cfO; .sort 

local FTUS ^FT + FU + FS; .sort 
polyfun cfO; .sort 

polyfun; id cfO(x?) = l; .sort 
print +s FTUS; .end 

The result of this simple FORM script is that FTUS is a sum of independent helicity 
structures u{\2„pz)f^''{l5',Pi)^ti{pi, <li)^u{p2, q2)v{X4,P4)- Now one can use Perl to give 
each term a name SME{Xi). For FFE, it is just slightly more complicated 
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local FT = ^(A)^; local FU = ^(A)^; local FS = ^(A)^; 

bracket u,v,e,^; .sort 

collect cfO; normalize cfO; .sort 

local FTUS ^FT + FU + FS; bracket cfO; .sort 

collect cfl; .sort 

polyfun cfl; .sort 

polyfun; id cfl(x?)=l; .sort 

print +s FTUS; .end 

where the command "normalize cfO;" is very important. The result of this FORM 
script is that FTUS is a sum of cfO[^{X s{pi.pj) x fTg)] where Xs{pi.pj) are just 
simple algebra expressions, fTg are loop functions. We can then use Perl to give the 
argument of each function cfO a name FFE. Thus FFE = J2{^s{Pi-Pj) x fTs). 

To sum up, the working stream of our FORM code is the following. Input: all 
expressions of helicity amplitudes for all Feynman diagrams. This is just simply 
applying the Feynman rules. Output: Y^T-diagrams A>^f , Eu -diagrams A^'f and 
"Ys-diagrams fuuctious of FFEs and SMEs. The major source of bugs is at 

the beginning when we type in the Feynman-rule-amplitude-expressions. The rest is 
almost automatic. One might wonder about the connection between FORM and Perl. 
This is semi-automatic in our code, i.e. we have to run FORM and Perl separately. 
The three FORM scripts described above generate their output files. We write three 
very simple Perl scripts to read those files and introduce abbreviations. Those Perl 
scripts also prepare three output files to be read by FORM. Indeed if one does not hke 
using Perl and wants to do everything automatically within FORM, this is possible. 



B.3 Automation with FORM 



In this section, we would like to show that the working stream described in the 
previous section can be automatized in FORM without invoking Perl. We have not 
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done this in the pp bbH calculation. However, the implementation is straight 
forward. The only difficulty is "How to introduce abbreviations?". The answer is in 
the following FORM examplq^- 



Symbol a,b,c,d,e,x,y,n; CFunction fl,f2,f3; 
Local F = x*(l+a+b+c)"3+y*(l+a+b+c)"2; 
AntiBracket x,y; .sort 

Collect fl; Makeinteger fl; Bracket fl; .sort 

*** which terms to be abbreviated? *** 

Keep Brackets; 

id fl(x?) = fl(nterms_(x),x); 

id fl(l,x?) = x; 

id fl(n?,x?) = fl(-termsinbracket_(0),x); 
id fl(-l,x?) = x; 
Bracket fl; .sort 

*** give it a name and store it by using $ variable *** 
Keep Brackets; 
#$cou = 0; 

if ( count (f 1,1) !=0 ); $cou = $cou + 1; id fl(n?,x?) = f2($cou)*f3($cou,x); 
endif; 

Bracket f2,f3; .sort 

#do i = l,'$cou' id f3('i',x?$t'i') = 1; #enddo .sort 

*** using temporary expressions for writing output file *** 
#do i = l,'$cou' local XX'i'=$t'i'; #enddo .sort 

*** for fortran output files *** 
format doublefortran; 

^write <abbrf2.F> " Subroutine abbreviation(x,y)" 

#write <abbrf2.F> " IMPLICIT DOUBLE PRECISION (A-H, 0-Z)" 



^We have learnt these tricks from a private communication with Vermaseren. 
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#write <abbrf2.F> " DOUBLE PRECISION f2('$cou')" 
#write <abbrf2.F> " Common/abbr/f2" 
#do i = l,'$cou' 

* #write <abbrf2.F> " f2('i') = '$t'i"" 
#write <abbrf2.F> " f2('i') = %e",XX'i' 
#enddo 

#write <abbrf2.F> " End" 

T^write <funct.F> " Function fun(x,y,a,b,c)" 

#write <funct.F> " IMPLICIT DOUBLE PRECISION (A-H, 0-Z)" 

#write <funct.F> " DOUBLE PRECISION f2('$cou')" 

#write <funct.F> " Common/abbr/f2" 

#write <funct.F> " fun=%e",F(fun) 

#write <funct.F> " Return" 

#write <funct.F> " End" 

.end 

The working stream of this example is the following. Input: an algebraic expression 
named F. Output: two Fortran files to calculate F: "funct.F" and "abbrf2.F". The 
latter computes all abbreviations which are "complicated" functions of (x,y) and 
appear several times in the final result. This is nothing but the idea of optimization. 
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Phase space integral 



c.l 2^3 phase space integral 





Figure C.l: A typical 2^3 Feynman diagram. The arrow gives the momentum 
direction. 



The phase space integral is given by 



d^P3 d^Pi d^P5 c4 



S {P1+P2-P3-P4-P5) 



(C.l] 



2e3 2e4 2e5 

There are 9 integration variables with 4 constraints from the Dirac delta function. 
All the interactions we consider are spherically symmetric, it means that there is one 
trivial variable corresponding to rotation around the z-axis. Integration over 
gives a factor of 27r. Thus, there are 4 essential final state variable^]. We define the 



^For a, 2 —> n process, the number of essential final state variables is 3n — 5. 
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kinematical function 



X{x,y,z) = + + — 2xy — 2xz — 2yz. 



(C.2) 



i?3 can be factorized into 2-^2 and 1 2 processes (see Fig. IC.lb ) 



i?3 



/ 



(iS45-R2(s, S45, m^)/?2(s45, "^fe, m^) , 



(C.3) 



with 



I 



(C.4) 



where 6{(f))l are the angles determined in the rest frame of {i + j), sij = {pi + Pj^ 
with i,j = 3,4,5. Clearly, formula (IC.Sp is just one way of factorizing the phase 
space integral. We can replace S45 by S34 or S35. In practice, choosing a good set of 
integration variables makes the integral convergent much faster. For a complicated 
calculation it is very difficult to know which choice is the best. In that case, we should 
always start with the tree level and try all the possibilities of phase space factorization 
with the same number of Monte Carlo points (if one uses the Monte Carlo method) 
and compare the integration errors to judge the best choice. That is what we did to 
find out that Eq. fIC.Sp is the best way to parameterise the phase space in the case 
of gg bbH calculation. 

We choose Oz \\ p2. In the center-of-mass system (CMS) of (4 + 5), we call this 
the CMS45 hereafter, one has P4 + P5 = Pi + P2 — P3 = 0. Thus p^, P2 and pg define 
a plane chosen to be Oxz (see Fig. IC.2p . P4 is defined by two angles 6*24 and 024- 




(C.5) 



In the CMS of (1 + 2), we call this the CMSgg hereafter, one gets 




(C.6) 



C.l. 2^3 phase space integral 



151 



P3 = Pi + P2 




Figure C.2: Coordinate system in the frame P4 + P5 = 0. 



where the factor 2n comes from the trivial integration over 023- The range of S2 is 
given by 

s> S2 >{mb + Mnf, 
So>s = X1X2S0 > (2mfe + MhY, (C.7) 

where sq is the invariant mass of the initial protons, xi^2 G [0, 1] are the momentum 
fractions carried by the initial gluons. The second equation in flC.7p implies that 

So SqXi 

The integration formula we actually use in the Fortran code, taking into account the 
convolution of the gluon structure functions, reads 



/»1 /•! /»1 /•! /»1 r2n 

Rs = 2tC dx[ / dx'2 / ds'2 / (icOS^^23 / (icOS^^24 / '^024, 

Jo Jo Jo J~i J~i Jo 



,S2,ml) X^^^{s2,ml,Mjj) 
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where we have changed the integration variables as follows 



Xi 
S2 



X2 



1 - 



(2mfe + MhY 

So 

(2mfe + Mg)^ 

SoXi 



+ 



+ 



(2m6 + MhY 
So 

(2mb + Mnf 

SqXi 



and J is the Jacobian 

(2mfe + Mnf 



J 



1 - 



So 



1 - 



(2mb + Mh) 

SqXi 



[s - {rrib + Mh) 



(C.IO) 



(C.ll) 



We should stress again that 624 and 024 are defined in the CMS45 while cos ^23 is 
defined in the CMSgg. In order to calculate the hehcity amplitudes, which are Lorentz 

invariant, one has to reconstruct from {s2, cos ^^23; cos ^^24, ^24} all the components of 
5 external momenta in some reference frame. The way we do this for the CMSgg is 
as follows. First, the components of p4 and can be easily calculated in the CMS45 

Xy\s2,mlMj,) 



P5 

P4z 
Piy 

In the CMSgg we get 



-P4> I P5 



P4 1 = 



2^ 



P4 I cos 6124, 

P4 I sin 6*24 COS 024, 

P4 I sin 6*24 sin 024- 



(C.12) 



Pi = (a/s, 0, 0, , p2 = (Vs, 0,0, y/s), 

\y\s,S2,ml) 



P45 



P4 + P5 = -Pa, I Pa 



2^s 



, 645 = VS2+ I Pa P, 



Pa^ = I P3 I cos ^23, P3^ =1 P3 I sin^23, Pa2/ = 0. 



(C.13) 



In the CMSgg one sees that the CMS45 is moving with the 4-component momentum 
P45 = (645,— Ps)- We then boost calculated above in the CMS45 to be in the 
CMSgg by using the following general Lorentz transformation 



7oe - 7oVo-P, 

,7oVo.p 



P = P + 7oVo( 



7o + 1 



(C.14) 
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with 

70 = —, mo = -y/eg -pg,vo = — (po^, Poj,, PoJ, (C.15) 
mo ^ Co 

where po is the 4-component momentum of a reference frame K' observed in K, p is 
any momentum observed in K and p' is the same momentum observed in K'. The 
inverse of equations (lC.14p and (1C.15P are obtained by changing vo to — vo and by 
interchanging primed and unprimed variables. Now we have all the components of 
Pi) P2, P3, Pi, P5 = Pi + P2 — P3 — Pi in one reference frame, CMSgg. 

In the CMS of two initial protons (CMSPP), the proton and gluon momenta are 

= (v/i^,0,0,v/i^) ,Pf = (v/i^,0,0,-v^), 
Pi = xiPi ,p2 = X2P2, (C.16) 

where we have neglected the proton masj^. All the kinematical cuts are defined in 
the CMSPP. One can move from the CMSgg to CMSPP by using the following boost 
matrix along the z axis 





7 








-Pi 


\ 







1 


















1 







v 


-Pi 








7 


/ 



(C.17) 



with 7 = ^i+fiL S = ^1 ^2 _ This boost matrix can be easily found from the Lorentz 
transformation of the four vector (pi + p-2) from the CMSgg to CMSPP. The helicity 
amplitudes can be calculated in the CMSgg or CMSPP as one wishes. 

For experimental purpose, one has to impose cuts on the transverse momenta and 
pseudorapidities of the bottom and anti-bottom in CMSPP 

IPstI ) IP4tI — |PT|min 

and |r73,4| < ?7max (C.18) 

^The proton mass — 0.9383GeV. 
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where the values of Ip^.] 
pseudorapidity is 



and r]jnax depend on the experiment, the definition of 



V 



In 



, e 

tan - 



(C.19) 



Those kinematical cuts also help to avoid some possible zero poles associated with 
some bottom-quark propagators in the massless limit and with some Gram determi- 
nants related to the tensorial reduction of loop integrals. 



C.2 Numerical integration with BASES 

BASES is a Monte Carlo integrator which functions by means of the importance and 
stratified sampling method [37J. Executions of BASES consists of the grid optimiza- 
tion and integration steps. Both steps are done by performing a number of iterations. 
An iteration is the process of computing the estimate of an integral and its variance. 
Each iteration is a Monte Carlo integration with Ncaii sample points and is realised 
as follows. The full multi-dimensional integral volume is covered by a grid of Ncube, 
the number of hypercubes (each hypercube is divided into many subregions). In each 
hypercube, the integral and its variance are evaluated with Ntriai = Ncaii/ Ncube sam- 
ple points. The results of each iteration are obtained by summing up results of all 
hypercubes. Ncube is calculated as folows 

(AT \ '^/^wild 
< 25, (C.20) 

with Nuiiid is the number of wild variables on which the integrand depends strongly 
or exhibits singular behavior. Ncaii and Nu,iid are BASES input parameters. Clearly 
Nwiid < Ndimi Ndim the number of integration variables, and the maximum number 
of wild variables is 15. As input information for BASES, one has to decide which 
are the "wild" variables and place them at the beginning of the integration variable 
array. 
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• Grid optimization step: at the first iteration the grid is uniformly defined for 
each variable axis. After each iteration the grid is adjusted so as to make the 
size of the subregions narrower at the parts with larger function value and wider 
at the parts with the smaller one. In this way a suited grid to the integrand is 
obtained. The number of iterations for this step is denoted ITMXl (default 
15), a BASES input parameter. 

• Integration step: the probability to select each hypercube and the maximum 
value of the function in it are calculated as well as the estimate of integral with 
the frozen grid determined in the former step. The number of iterations for this 
step is denoted ITMX2 (default 100), a BASES input parameter. 

The typical BASES input parameters for calculating (IC.9P are: Ndim = 6, N^ud = 2, 
Ncaii = 10^ ITMXl = 20 and ITMX2 = 130. The two wild variables are cos ^23 
and cos 62A placed at the beginning of the integration variable array. With those input 
parameters, the typical error we obtained for the pp bbH calculation is 0.08%. 
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D.l Logarithms and Powers 

The natural logarithm In (2;) is defined as 

ln(z)=ln(|z|)+iarg(z), (D.l) 

with — TT < arg(2;) < tt. The logarithm In 2; has a branch cut along the negative real 
axis. The general power w — z"' {a is a, complex constant) is defined with the aid of 
the exponential function 

With those definitions, one has the following rules 

ln{ziZ2) = \n{zi) +\n{z2) +r]{zi,Z2), (D.3) 
77(^1, Z2) = 2'jTi[9{— Im zi)9{— Im 2;2)^^(Im -212:2) — ^^(Im 2;i)6'(Im Z2)9{— Im Z1Z2)], 
{Z1Z2)" — e"''^*-^^^^^ = eO^[in(^i)+iii(^i)+'?(2i'^2)] _ e"''(^i'^2)^a^a^ (D-4) 
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\n{-z) 



which have important consequences 

\n{ziZ2) = l^{zi) + ln(2;2) if Im Zi and Im Z2 have different sign 
In — = In(zi) — ln(z2) if Ini Zi and Im Z2 have the same sign 
(-21-22)" = z'^Z2 if Im Zi and Im Z2 have different sign. (D-5) 
For —z = a — ip with a real and p — 0^ we have 

In \a\ if a > 

In \a\ — iir if a < 

aYg[—z) = aig{z) — n, 

\n{—z) = ln(|z|) + i arg(— z) = ln(|z|) + i arg(z) — in = ln(z) — in, 
(-^)" = e-^^e"'''^^) = e-*™z". (D.6) 

If A and B are real then 

\n{AB - ip) = \n{A - ip') + \n{B - ip/A), (D.7) 

where p' is infinitesimal and has the same sign as p. From this we get 

{AB - ip)'' = e'^''^(^^-^^) = e''[HA-ip')+HB-ip/A)] = (^A- ipT{B - ip/A)''. (D.8) 



D.2 Dilogarithms 



The dilogarithm or Spence function is defined by |119l 1124^ 1125] 

where z may be complex. The logarithm has a branch cut along the negative real 
axis, implying for the Spence function a cut along the positive real axis from 1 to 
+00. When one is in a problematic situation, the following transformation formulae 
may be helpful 

Sp(^) = -Spi-)-]:n'-l\n\-z), (D.IO) 
z b Z 

Sp(z) = -Sp(l-z) + -n^ -lnn-z)\n(z). (D.ll) 

6 
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More transformation formulae can be found in |124l I125j . 



D.3 Gamma and Beta functions 



The gamma function T{z) is a function of the complex variable z. For Re 2 > it is 
defined by 



r(z) 



dte'H 



t,z-l 



(D.12) 



where the principal value of ^ is to be taken. For Re 2; < the alternative definition 
reads |126] 

1 



where the path of integration C starts at —00 on the real axis, circles the origin once 
in the positive direction, and returns to —00; the initial and final arguments of t are 
to be — TT and vr, respectively. The latter defines an analytic function for all z other 
then 0, ±1, ±2, .... For positive integers n definition (ID. 121) gives 



r(n + 1) = nT[n) = n\. 



(D.14) 



T(z) is analytic everywhere, except at the points z = 0, —1, —2, .... The following 
properties of the gamma function are very useful 



T{z + 1) = zT{z) 
vr 



r ^ F 1-;^ 



Sm TTZ 

r(e) = l-7^ + 0(e), e 
e 



(D.15) 



where 7_e is Euler constant. 

The beta function is defined by 



B{p,q) = / dttP-\l-t) 
'0 



(D.16) 
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27]-A72 

dti--- dtK = I r^-^drdQK-i, I dQx-i = tttttt^- (D.18) 



where Rep > and Reg > 0; the principal values of the various powers are to be 
taken. The analytic continuation of B[p, q) onto the left halves of the p and q planes 
is achieved by using 

D.4 Integrals 

Formulae to move to spherical coordinates: 

T{K/2) 

The following integral formula is very useful in many cases 

{z + sr r(/5) ' ^ ^ 

where z can be complex. 

The following integral usually appears in loop calculation |119] 

^3 = r dy ^'^^^'^^ + by + c)- \n{ay^ + byp + c) ^g) 

Jo y~yo 

where a is real, while b, c and yo may be complex, with the restriction that lm{ay^ + 
by + c) has the same sign for < ?/ < 1. 

Let e and p be infinitesimally real quantities having the opposite sign to Im(a?/^ + 
by + c) and Im(a|/o + byo + c) respectively. Using Eq. (ID. 71) we get 

ln(a?/2 + by + c) = ln(a - ie) + \n[{y - yi){y ~ ^2)] , 
\n{ayl + byo + c) = ln(a - zp) + ln[(?/o - 2/i)(yo - 2/2)], (D.21) 

where 7/1^2 are two roots of equation 

y^ + -y + - = 0. (D.22) 

a a 



D.4. Integrals 
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We then use Eq. (1D.3P to get 

^3 = 



^ dy ^^^^^ ~ ^^^'^^ ~ ^^^^ ~ ~ ^^"^^^^ ~ ^^^^ 

y-yo 



— rj{a — ie, 



a — ip 



In 



(D.23) 



We write 



ln[(2/o - yi){yQ - 2/2)] 
ln[(z/-l/i)(z/-l/2)] 



ln(?/o - Z/i) + ln(?/o - y2) + ?7(yo - Z/i, Z/o - 1/2) 

in(?/ - yi) + in(?/ - 1/2) + viy -yuy - ^2) (D.24) 



with 



^(y - 1/1, y - 1/2) = ^(-2/1, -2/2) 



(D.25) 



since y is real and Im[(y — yi){y — 2/2)] = + = -'■^(a) ~ Ii^(z/i?/2) as we 

assumed at the beginning. We have 



^ ^y ^y - yi) - Myo - yi) + My - y^) - in(yo - ^2) 



2/ - I/O 



+ 



-2/2) - v{yo - yi, yo - 2/2) - vi^ - ie, — -) 



a — ip 



For this, we have to calculate 

R{yi,yo) = 



- \n{yo-yi) 
2/ - 2/0 



In^^ -. (D.26) 

2/0 



(D.27) 



We change the integration variable y = y' + yi to get 

R = 



(D.28) 



-s/i 2/ - 2/0 + 2/1 

Since the residue of the pole is zero and the logarithmic cut along the negative real 
axis is outside the triangle [0, — yi, 1 — we can write 

/-i-j/i r-yi 

yi Jo Jo 



(D.29) 
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We then make the substitutions ?/ = (1 — yi)y' and y = —yiy' to get 

i-yi 



R 







fo ^^dy 


In ( 






fo ^^dy 


In ( 



yi -vo 
yi 



{ln[(l -yi)y\ -ln(?/o-yi)} 
{\n{-yiy)-\n{yo-yi)]. (D.30) 



yi - 2/0, 

Since y is real and positive the logarithmic arguments never cross the cut along the 
negative real axis. We do partial integration to get 



1-2/0 



R = in( [ln(l-yi)-M2/o-2/i)]-ln ^^ [ln(-|/i)-ln(yo-2/i)] 

,2/1-2/0/ V2/1-2/0/ 



-2/0 



Sp 



2/1 



Sp 



1-2/1 

.2/1 -2/0/ ' V2/0 -2/1 
One uses Eq. (1D.11I) to obtain 

2/0 



:D.31) 



^(2/1,2/0) = Sp 
- In 



2/0 - '2/1 
1-2/0 
2/1 - 2/0 



Sp(^^l+ln 

2/0-2/1 

1 



2/0 



2/0 - 2/1 



v[-yu 



^(1 - 2/1 



2/0 - 2/1 ■ 



2/0 - 2/1 

(D.32) 



The result for S^i reads 



5'3 = i?(2/i, 2/0) + ^(2/2,2/0) 



+ 



^(-'2/1, -2/2) - v{yo - 2/1, yo - 2/2) - - ^e, 



1 



a — ip 



— ^ (D.33) 
2/0 



which contains 4 Spence functions. 



Appendix E 

Scalar box integrals with complex 
masses 

The calculation of the scalar one-loop function for the box (N = 4) with imaginary 
internal masses in the most general case with no restriction on the external invariants 
is not tractable. The standard technique of 't Hooft and Veltman [119J (see also [127] ) 
has some restriction on the values of external momenta. In particular, the method 
works if at least one of the external momenta is lightlike. In our present calculation, 
there are at least 2 lightlike external momenta in all boxes. We explain here our 
derivation based on the method given in [119] for this special case. 

With N = D = A, from Eq. fH35ll we get 



Do ^ i^^m 



dx 



dy 



y 



dz 
1 



X 



{ax"^ + hy"^ + gz^ + cxy + hxz + jyz + dx + ey + kz + /)^ ' 



(E.l) 



where we have changed the integration variables as t = J2i=i^ij ^ = J2i=i^i^ V = 
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2Q34)=pI, b=-{Q22 + Q33-2Q23)=pl, 

2Q12) =pI, C = Q23 + Q3A - Q33 - <524 = 2p2-P3, 
h = Ql3 + Q24 - Ql4 - Q23 = 2pi.p3, j = Q12 + Q23 - Q22 - Ql3 = 2pi.p2, 

d = Q34- Q44 = ml - ml - pI, e = Q24 - Q3A = ml - ml - pi- 2p2.p3, 

k = Qu - Q24 = ml-ml+pj + 2pip4, f = -ie = ml- ie, (E.2) 

with Qij is defined in Eq. fl4.25l) . d, e, k, f are complex while other parameters are 
real. There are two cases corresponding to the fact that the positions of two lightlike 
momenta are opposite or adjacent. 



Xi + X2, z = Xi, and 
« = ^{Q33 + Qm- 

9 = ^(Qll + <522- 



E.l Integral with two opposite lightlike external 
momenta 




Figure E.l: A box diagram with two opposite lightlike external momenta pi and ps. 
Double line means massless. 



For the box shown in Fig. IE. II with Pi=pl = one gets a = g = and writes 

p1 px py 1 

D^^'^^ = dx dy dz—- : — . (E.3) 

Jo io Jo [by + cxy + hxz + jyz + dx + ey + kz + f )^ 



E.l. Integral with two opposite lightlike external momenta 
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Integrating over z to get 



Jo Jo 
with 

A = cy + d, B = by^ + ey + f, 

C = (c + % + d, L' = (6 + jV + (e + % + /. (E.5) 



One changes the integration order as 



1 fX pi pi 

dx dy = dy dx. (E.6) 
Jo Jo Jy 



We get 



= / dyy [ dx ^^^ , , (E.7) 

Jo Jy 



1 

' (Ax + B){Cx + D) ' 



where A, B, C, D are complex. Integrating over x as follows 

"1 1 1 /■! 



(dx = J- / 

jy iAx + B){Cx + D) ACjy (x + f)(x+g) 



dx 



1 /-W 1 





BC 


1 




AD- 


BC 


1 




AD- 


BC 



dx 



, x + - x + - 

1 + B 
In ^ - In g 

+ S Cy + Dj ' ^ ' 

where wc have made sure that the arguments of the logarithms never cross the cut 
along the negative real axis. One easily gets 

Df' = fdy 
Jo 

^ {cj - bh)y^ + {dj + ck- eh)y + dk - fh \ ^ Ay + B ~ ^ Cy + D J ^ 
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where the discriminant of the quadratic denominator in the prefactor is nothing but 
the Landau determinant 

det(Q4) = {dj + ck- ehf - 4(cj - bh){dk - fh). (E.IO) 

We write 

,(13) 1 /"Vl I \ A + B ^ C + D 



icj-bh)iy2-yi)Jo [y-y2 y-yJ^'^Ay + B Cy + D ' ^^'^^^ 

with 



-{dj + ck - eh) T ^/det{Q4} 

y^'^ = 2(^1^) ' ^^•'2) 

where the indices 1, 2 correspond to — and + signs respectively. 

Now we have to look at the imaginary parts of the arguments of the logarithms 
in Eq. (lE.lip . We write them explicitly 



A + B = by^ + {c + e)y + d + f, 

Ay + B = {b + c)y^ + {e + d)y + f, 

C + D = {b + j)y'^ + {e + k + c + h)y + d + f, 

Cy + D = {b + j + c + h)y^ + {e + k + d)y + f. (E.13) 

Imaginary parts read 

lm{A + B) = lm{ey + d + f) = lm[yml + {1 - y)ml- ie] < 0, 

lm{Ay + B) = lm{ey + dy + f) = lm[yml + (1 — y)m1 — ie] < 0, 

lm{C + D) = Im[{e + k)y + d + f] = lm[yml + {l-y)ml-ie]<0, 

lm{Cy + D) = lm[{e + k)y + dy + f] = lm[yml + {1 - y)ml- ie] < 0. (E.U) 

Using formula ln(a/6) = In a — Info for Im(a) Im(6) > 0, we rewrite Eq. (IE. lip as 

-,2 4 /•! 1 

dJ"^ = , y Yi-iy^' / dy^- HA.y^ + B.y + C.) (E.15) 



i=i j=i 



E.l. Integral with two opposite lightlike external momenta 
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with 



A, 
A2 
A3 
A^ 



b + c, Bi = e + d, Ci = f, 

b, B2 = c + e, C2 = d + f, 

b + j, B3 = e + k + c + h, C3 = d + /, 

b + j + c + h, B4 = e + k + d, C4 = /• 



(E.16) 



We would like to make an important remark here. From Eq. (]E.14I) we can re-write 
Eq. (jE.lip in the form 



(13) 



1 



1 



1 



In 



A + B 



° (cj - bh){y2 - yi) Jq \y - y2 y - Vi ) K" C + D 
We notice that if y = yi.2 then AD = BC which means 



In 



Ay + B 
Cy + D 



A + B 



Thus, we get 



1 



C + D 



1 



Ay + B 



y=vi.2 



Cy + D 



y=yi,2 



B 
D 



y=yi,2 



In- 



A + B 



In 



Ay + B 



Cy + D 



y=yi,2 



^y-y2 y-yij \ C + Dy=yr,2 

Subtracting this zero contribution from Eq. flE.lSp we get another form 

2 4 



0. 



(E.17) 



(E.18) 



(E.19) 



D. 



(13) 



1 



i+j 



/det(g4) ,=1 ,=1 
^ j^y ^^^y' + B,y + Q) - HA.yl + B,y, + C,) 

lo y-yi 



+ 



1 



i/det((54 



X^(-l)^+^r^.,ln 



yi 



Vi 



\y=yi 



(E.20) 
with 



where T]a = v{A + B, 1/{C + D))\y=y^ and T]i2 = T]{Ay + B, l/{Cy + D))\, 
i = 1,2. This representation is more convenient for the evaluation in terms of Spence 
functions. 



Each integral in Eq. (iKlSj) or (|K20|) can be written in terms of 4 Spence functions 
as given m Eq. flD33l) . Thus D^^^^ can be written in terms of 32 Spence functions. 
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E.2 Integral with two adjacent light like external 
momenta 





Figure E.2: A box diagram with two adjacent lightlike external momenta pi and p2- 
Double line means massless. 



For the box shown in Fig. IE. 21 with Pi = P2 = ^ one gets b = g = and writes 



D 



q"^^-* = I dx I dy I dz . ^ , . , , „,„ 

Jo Jo {(^^ + '^^y + '^^^ + Jy^ + dx + ey + kz + J Y 

)(13) 



. (E.21: 



As in the case oi Dq , integrating over z gives 



with 



(12) 



X 



1 



- I dx dy—-+Sk / dx / dy — — , 

Jo Jo Jo -Skai{aiy + bi) 

V ' V ' 

h I2 



1 



(E.22) 



Sk = sign{lm.{k)) , —Skai = —Sk{hx + jy + k)—ie', 
bi = ax^ + cxy + dx + ey + f, 
o-iy + bi = ax'^ + jy"^ + {c + h)xy + dx + {e + k)y + f — ie, (E.23) 

where we have used the fact that lm{aiy + 61) = Im[(ix + (e + k)y + /] = Im[(x — 
y)ml + (1 — x)m1 + yml — ie] < because < y < x < 1. e and e' are infinitesmally 



E.2. Integral with two adjacent lightlike external momenta 
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positive and carry the sign of the imaginary parts of — s^ai and aiy + bi. For /i, we 
integrate over y, similar to Eq. (lE.SI) . to get 



h 



dy 



{ja — hc)y'^ + {jd — he — kc)y + jf — ke 



(j + h)y + k- te' _ (a + c)y^ + (d + e)y + f 



hy + k — ie' 



ay"^ + dy + f 



(E.24) 



Consider the prefactor 



det{Qi) = {jd — he — kcY — A{ja — hc){j f — ke), 
{he + kc — jd) =F ^/det{Q4) 

2/11(12) = 



(E.25) 



2{ja — he) 

where the indices 11, 12 correspond to — and + signs respectively. We rewrite Ji as 

2 



dy- 



(j + h)y + k- le' _^^{a + c)y^ + {d + e)y + f 



y-yii[ hy + k-ie' 



2 4 



Vd^t{Q^~t~t Jo y-yu 



dy- 



ay"^ + dy + f 



ln{A,0^ + Byy + Cy) (E.26) 



with 



^11 

^12 
^13 
Ai4 



0, Bu = h, Cn = K 
0, B^2=J + h, Ci2 = k, 
a + c, Bi3 = d + e, C13 = /, 



a, B 



14 



d, Ci4 — /. 



(E.27) 



Thus Ji can be written in terms of 24 Spence functions. For I2 we shift y = y + ax, 
a such that 



ja^ + (c + h)a + a = 0. 



(E.28) 



There are, in general, two values of a. The final result does not depend on which value 
of a we take. We have used this freedom to find bugs in the numerical calculation and 
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it turns out to be a very powerful method to check the correctness of the imaginary 
part which can be very tricky for the case of equal masses. One gets 

72 = / dx / dy— — -, (E.29) 

Jo J -ax {Gx + H - ie'){Ex + F - iey ^ ' 

with 

G = -Skh - Skja, H = -Skjy - Skk, 

E = {2ja + c + h)y + d + a{e + k), F ^ jy^ + {e + k)y + f. (E.30) 
For real a we have 



rl /•(l—a)x i-l /•(l—a)x i-l n—ax 

dx dy — dx dy — dx dy 

Jo J-ax Jo Jo Jo Jo 

1*1— a 1*1 I* — a j*l 

= / dy j dx- j dy j dx. (E.31) 

Jo Jy/(l-a) Jo J-y/a 



ly/{l-a) Jo J-y/c 

We write 

1 1 / G E 



{Gx + H - ie'){Ex + F - ie) GF - HE \Gx + H - ie' Ex + F - ie 
Integrating over x, we get 



. (E.32) 



GF-HE E + F Jo GF-HE ^ + F 
Jo GF-HE ^ + F ^ ^ 



The prefactor 
GF-HE 

— + c)?/^ + (2ajk + jd — he + kc)y + a{ke + k"^ — jf) + kd — hf 

= i(ia + c)(y-Z/2i)(y-2/22), (E.34) 



with 



-{2ajk + jd-he + kc) ^ y/detiQ^) 

2/21(22) = ov ^ ^ ' ^-2^ 

2j{ja + c) 
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where the indices 21, 22 correspond to — and + signs respectively. We rewrite I2 as 



h = — -i — 

Sk^/det{Q4) 

J-a y-y2i E + F io y-y2^ ^ + f 

—^\n^ . (E.36) 

y-y2^ ^ + F ^ 



We make the substitutions y = y — a for the first integral, y = [1 — a)y for the second 
integral and y = —ay for the third integral to get 



dy , -Skjy - Skh - Skk - ie' 



= / In 



y -a-y2i jy'^ + {c + h + e + k)y + a + d + f - ie 

^ {l-a)dy -Sk{j + h)y - Skk - ie' 

In ■ 



{l-a)y-y2i {a + c + j + h)y'^ + {d + e + k)y + f - ie 

—ady , —Skhy — Skk — ie' ,^ 

+ / — 2Z, If ■ ■ (E-37) 

Jo -ay - y2i ay^ + dy + f - le 



Consider the arguments of the three logarithms, as demonstrated in Eq. ( lE.14p . it 
is easy to see that the sign of the imaginary parts of the denominators is negative 
as indicated by —ie. The derivation is for real a. However, this result can be easily 
generalized to cover the case of complex a as shown below. We can now rewrite I2 as 



/i^TTTTt EE^-^)^ f'^y T HA2,y' + B2,y + C2,) (E.38) 
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with 





= 1, Oi = 1, 61 = a, 




Co 


= — fl — ot). flo = 1 — a. 6o = 0. 






= — a, 03 = — tt, 63 = 0, 




Ca 


= — 1, (34 = 1, b/^ = ot, 






= 1 — a, 05 = 1 — a, 65 = 0, 






= tt, (36 ~ — Ct) ^6 ~ 0, 






= _Boi = — 7 (^oi = — ,Hi.k — 


Huh, 


A22 


= 0, B22 = -Sk{j + h), C22 = - 


7 

Skit, 


^23 


= 0, B23 = -Skh, C23 = -Skk, 




^24 


= j, B2A = c + h + e + k, C24 = 


a + d + f. 


A25 


= a + c + j + h, B25 = d + e + k. 


C25 = /; 


^26 


= a, B2Q = d, C26 = /• 





I2 can be written in terms of 36 Spence functions. Thus 

4''^ = h + Skh (E.40) 

contains 60 Spence functions. For the evaluation of D^^^ in terms of Spence functions 
and to generahze Eq. ( IE. 381) for complex a, we have to do the following replacement 
for each logarithm in Ji 2: 

ln(v4ij?/2 ^ Q^,y ^ ln(Aij?/2 ^ Q^^y ^ Q^^^ _ \n{Aijyf. + B^^yu + Cy), 

ln{A2jy^ + B2jy + ^2^) ^ ln(A2,i/' + 52,1/ + C2j) - HA2jyli + B2jy2i + C2,), 

with ^2j = {y2i + bj)/aj and add the corresponding extra terms related to the eta 
functions. The argument for this is similar to that explained in the previous section, 
see Eq. flKM . 

For the boxes with one lightlike external momentum, the result is written in terms 
of 72 Spence functions by using exactly the same method. 
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